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Abstract 

Introduction. Vibration interaction control is timely in production processes related to liquid and bulk media, systems of 
solids experiencing kinematic or force disturbances. At the same time, there is no single methodological basis for the 
formation of vibrational interactions. The issues of constructing optimal vibration fields of technical objects have not 


been addressed. The objective of the study is to develop a structural approach to the development of mathematical models 
in the problems of formation, evaluation, and correction of vibration fields of technical objects under conditions of intense 
force and kinematic loads. The task is to build vibration fields that are optimal in terms of the set of requirements, with 
the possibility of selecting the criterion of optimality of the vibration field of a technical object. 

Materials and Methods. A structural approach was used as the basic methodology. It was based on a comparison of 
mechanical vibratory systems used as computational schemes of technical objects, and structural schemes of automatic 
control systems, which are equivalent in dynamic terms. Lagrange formalism, elements of operational calculus based on 
Laplace integral transformations, sections of vibration theories, algebraic methods, and the theory of spline functions 
were used for structural mathematical modeling. 

Results. An approach to the selection of criteria for the optimality of vibration fields based on minimizing the residual of 
vibration fields for various required conditions was proposed. The problem was considered within the framework of a 
mechanical vibratory system formed by solids. It was shown that the optimal vibration field was determined by an external 
disturbance and was to satisfy condition Ay = b. There, A — matrix mapping the operator of conditions to the shape of 
the vibration field at control points; 5 — vector of values of vibration field characteristics; “—” above y meant the vibration 
amplitude of the steady-state component of the coordinate. To evaluate the field with account for noisy or unreliable 
requirements for dynamic characteristics, the smoothing parameter was used, indicating the priority of the criterion of 
optimality of the vibration field shape. The construction of a field for a mechanical vibratory system showed that the 
value of the vibration amplitudes of generalized coordinates remained constant when the frequency of external kinematic 
disturbances changed. Two approaches to the correction of the field optimality criteria were considered: equalization of 
the vibration amplitudes of the coordinates of a technical object and the selection of an energy operator. 

Discussion and Conclusion. The development of the applied theory of optimal vibration fields involved, firstly, the 
correlation of the energy operator and the operator of the requirements for the shape of the vibration field in the theory of 
abstract splines. The second pair of comparable elements was the criterion of optimality of the vibration field and a system 
of requirements for the characteristics of the field at control points. The structural theory of optimal vibration fields 
improved in this way will find application in various industries. Accurate calculations in the formation, assessment, and 
correction of the states of systems under vibration loading are required in the tasks of increasing the durability of 
structures, improving measurements in complex vibratory systems, and developing new technologies and materials. 


Keywords: structural mathematical modeling, mechanical vibratory systems, optimal vibration field, minimizing the 
residual of vibration fields 
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AHHOTalna 


Beedenue. Yupapyenve BUOpalMOHHbIMH B3aMMOJeHCTBUAMH AKTYAJIbHO B TIPOH3BOJCTBCHHBIX TIPOeccax, CBA3AHHBIX C 
2KMKMMH HM CbIITyYYHMH cpeyaMu, CHCTeMaMH TBeEpAIX TeJI, HCIBbITbIBAIOWIMX KHHEMaTH4eCKHe WIM CHJIOBbIeC BOSMYLICHHA. 
IIpu 9TOM HeT eAHHOM MeTOAOOrM4eCKOM OCHOBbI JJIt (POPMHpOBaHHA BHOpaHOHHBIX B3aHMoOselcTBHH. He pemieHsi 
BOMPOCbI MOCTPOeHHA ONTHMAJIbHBIX BHOpalMOHHBIX WouleHi TexHHuecKHX OOBeKToB. Llenb uccueqoBaHHa — pa3BHTHe 
CTpyKTypHoro Moyxoya K pa3paOoTKe MaTeMaTHYecKHX MOJeIeH B 3afayax (OpMHpoBaHHA, OIeCHKH MW KOppeKIHH 
BHOpalMOHHBIX Moe TeXHHYeCKHX OOBCKTOB B YCJIOBHAX HMHTCHCHBHBIX CHJIOBbIX H KHHeMaTHYeCKHX Harpy2KeHHit. 
CraButca 3ayaua MOCTpoeHHA BHOpalMOHHBIX Mose, OMTHMaJIbHBIX TIO COBOKYHHOCTH TpeOoBaHHi, C BO3MO%KHOCTBIO 
BbIOOpa KPHTepHA ONTHMAJIbHOCTH BHOpal{MOHHOTO MOJIA TEXHHYECKOTO OObEKTA. 

Mamepuaanei u memooovi. B xayectBe Oa30B0H MeTOAOIOrHH MIpHMeHAeTCA CTPYKTYPHbIM TozxoA. OH ocHoBaH Ha 
CONOCTABJICHHH MCXaHH4eCKHX KOJeOaTeIbHBIX CHCTEM, HCIOJIb3YCMBIX KaK PaCUeCTHbIe CXEMbI TEXHHYCCKHX OOBEKTOB, 
WM CTPyKTypHbIX CX€M CHCTeM aBTOMATHYeCKOTO yIpaBJICHHA, IKBUBAJICHTHBIX B JMHaMM4ecKOM oOTHOWeHHH. Jia 
CTpyKTypHOrO MaTeMaTHYeCKOrO MOA eMpOBaHHA HCHONb30BaIM PopMaiH3M JlarpaH>Ka, 9JIEMeHTbI OMepalHOHHoro 
MCUHCIICHHA Ha OCHOBe HHTerpasibHbIxX MpeoOpa3soBanui Jlanmaca, pa3qerbl Teopul KoNeOaHNi, anreOpanueckne 
MeTOJbI, TeOpHro criaitH-PyHKUHi. 

Pe3ynemamoi uccnedosanua. Upeanoxen Moxxoy K BbIOOpy KpUTepHeB ONTHMaJIbBHOCTH BHOpalMOHHBIX Mose Ha 
OCHOBe MHHMMH3alHH HeBA3KH BHOpallMOHHBIX ose WIA pa3svIM4HbIX HeoOxoAMMBIX ycoBui. I[poOmema 
paccmMaTpuBaeTCA B paMKax MexXaHH4eckKOH KOJIeOaTeIbHOM CHCTeMBI, OOpa30BaHHO TBEpALIMH TesaMu. loKa3aHo, 4TO 
ONTHMasIbHOe BHOpalHOHHOe Mose ONpeeAeTCA BHCLUIHHM BOSMYLICHHeM HU JOJNKHO YOBJICTBOpATS ycuoBH1o0 Ay = b. 
3yecb A — MaTpula, oToOpaxarollad olepaTop ycoBuit Ha (bopMy BHOpalMOHHOTO MONA B KOHTPOJIBHBIX TOUKAX; 
b — pBextop 3HayeHHii xXapakTepHcTHK BHOpalMOHHOrO MOA; «—» Ha y O3Ha4aeT aMIMIMTyy KoeOaHna 
YcCTaHOBUBIIelHcad KOMIOHeHTbI KOOpAMHatel. JI OWeCHKH TOA C Y4eTOM 3alLIyYMJICHHbIX WJM HeOCTOBepHBIX 
TpeOoBaHHit K WHHAaMMYeCKHM XapaKTepHCTHKaM HCMONb3yeTCA MapaMeTp criaKHBaHHA, OOOZHAYaIOLWIHH MpHoOpuTer 
KpHTepHaA ONTHMaIbHOCTH (opMbI BuOparHMoHHoro mond. TlocrpoeHue moma Id MexaHMuecKoH KoeOaTebHOH 
CHCTeMBI NloKa3aI0, YTO 3HaYeHHe aMIVIMTy KoeOaHHuaA OOOOMICHHBIX KOOPAMHaT COXpaHAeTCA MOCTOAHHBIM IIpv 
VW3MCHEHHM YaCTOTbI BHCINHHX KMHe€MaTHYeCKHX BO3MYINeHHH. PaccMoTpeHb!l WBa MOAXOZa K KOppeKUMH KpuTepueB 
ONTHMAJIbHOCTH O0JId: ypaBHHBaHWe aMIVIMTy{ KoNeOaHHit KOOpaMHaT TexHWYecKoro OOBeKTa HU BbIOOp 
9HepreTHyecKoro OllepaTopa. 

O6cystcoenue u 3akmo4enue. PasBuTve IpUKIaHON TeOpHu ONTHMAJIBHBIX BHOpalMOHHBIX NWoOsel Wpeslouaraert, Bo- 
TI€pBbIX, COMOCTaBJIeHHe OMepaTopa 39HeprHH HM OnepaTopa TpeOoBaHHii K opMe BuOpaltHOHHOrO MONA B TeOpHHu 
aOcTpakTHBIx ciulaiHoB. Bropas apa COMOCTaBJIAeMBbIX 9JIEMCHTOB — KpHTepHit ONTHMaJIbHOCTH BHOpal{MOHHO"O MOA 
WM CHcTeMa TpeOoBaHHit K XapaKTepHCTHKaM T1OJIA B KOHTPOJIBHBIX TOUKaX. YCOBepIeCHCTBOBaHHad TaKHM OOpa30M 
CTpyKTypHaai TEOPHA OMTHMAJIbHBIX BUOpalHOHHBIX Mosel HalieT IPHMeHeHHe B pa3HbIxX OTpacsAx. TOUHbIe pacuerTEl B 
(POpMUpoOBaHHH, OWeHKe HU KOPpeKUMH COCTOAHHM CHCTeM IPH BHOpal{MOHHOM Harpy2KeHuH HeOOxOAMMBI B 3aqayax 
MOBBILNCHHA JOUTOBCYHOCTH KOHCTpyKUMi, yIYUNeHHA H3MepeHHi B CIIOKHBIX KOeOAaTeIbHbIX CHCTeMax, pa3paooTKe 
HOBBIX TEXHOJIO“HH HW MaTepHasIOB. 


Korrouesbie cJI0Ba: CTpyKTypHoe MaTeMaTHYecKoOe MOICJIMpOBaHHe, MexaHHdeckKHe KoeOaTeIbHBIe CHCTCMBI, 
ONTHMaJIbBHOe BHOpallHwoHHoe TlOI€, MMHUMU3allHA HEBASKH BHOpallHOHHBIX nose 


BaarogxapHocrn. ABTOpBI BbIparKaroT OaroapHOCTE petakWHv U peweH3eHTaM 3a BHHMAaTeCJIbBHOe OTHOIMICHHE K CTaTbe 
VW yYKa3aHHble 3aMCuaHHA, KOTOPbIe MO3BOJIMJIM NOBbICHTb Ce KadecTBO. 
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Introduction. Features of vibration interactions are taken into account when creating production technologies [1]. 
The use of vibration methods for the intensification of production processes requires the development of mathematical 
modeling of the dynamic states of technical objects under conditions of vibration loading. The scientific and applied 
literature studies issues of the theory of vibrations and automatic control, theoretical mechanics, and dynamics of 
machines [2]. In particular, in the dynamics of machines, features of vibration interactions are taken into account in two 
types of tasks. The first one is the tasks of vibration isolation and vibroprotection [3]. The second includes tasks related 
to the assessment, formation, and correction of the states of technical objects under conditions of intense force or 
kinematic loads [4]. In the latter case, attention is drawn to the effects of vibrational interactions, which are used in the 
following production processes: 

— transportation of rocks; 

— fractionation of bulk materials; 

— hardening of surfaces of long elements of aircraft structures; 

— spatial orientation of parts during assembly. 

Tasks of vibroprotection are related to solving the problems of railway transport safety. It is referred, specifically, to the 
issues of the dynamic state of the suspension elements [5], noise and vibration [6]. The problems of estimating dynamic states 
caused by significant loads arise when evaluating the vibration characteristics of a rail track with composite sleepers [7]. Harmful 
vibrations propagating in the urban environment from a moving train require control, and in this case, an assessment of complex 
dynamic states is also needed [8]. Similar examples are related to tracking the interactions of a wagon trolley with rails [9]. The 
dynamic state is formed when the sliding contact of the pantograph head and the contact wire is violated [10]. In this context, it 
is worth mentioning the tasks of analyzing the reaction of the rail to vertical impacts from a moving vehicle [11]. 

Another extensive group of tasks requiring an assessment of the dynamic states of elements is the vibration processes of 
production machines [12]. Thus, the dynamic state of vibration screening machines depends on the formation of dynamic states 
of several housings, which can be considered as solids, taking into account the nature of elastic linkage between them [13]. An 
essential factor determining the dynamic states of production facilities is the vibration frequency [14]. Separately, it is necessary 
to mention vibrating transport machines, whose dynamic state determines characteristics of movements [15]. 

It is worth noting the result of generalization of the tasks of assessment, formation, and correction of dynamic states 
of technical objects under conditions of intense loads. This approach has opened up the possibility to clarify the idea of 
how the distributions of the amplitudes of vibrations in the coordinates of control points depend on the connectivity of 
external disturbances. 

Parameters of external vibration disturbances determine the possibilities of implementing dynamic effects that can act 
in different directions: 

— constitute a danger to the performance of technical facilities; 

— increase the process performance. 

An example of a negative effect is exceeding the permissible values of the vibration amplitudes. Useful effects include 
a significant decrease in the amplitudes of coordinate vibrations and the establishment of a certain amplitude distribution. 

In the tasks of vibration hardening of long-length parts, dynamic effects can be understood as the coincidence of the 
vibration patterns of the control points of the working bodies of vibrating production machines. This maintains the one- 
dimensionality and uniformity of the vibration field, which is consistent with the requirements for the dynamic quality of 
the interactions of the working medium and the surface of the hardened part. The required motion mode of the working 
body can be expressed using the conditions of equality of the vibration amplitudes of the control points. To study the 
concept of a vibrational field, the concept of a dynamic state at a point is generalized — at the same time, a set of points 
distributed over a certain variety of points of a solid body or a system of solids is considered [16]. 

In problems of vibrational interactions, a number of dynamic effects can be displayed using lever, elastic, and 
interpartial unilateral constraints. All of them are factors of specific dynamic effects that manifest themselves under the 
influence of external disturbances of a special configuration. To determine such configurations of external disturbances 
providing dynamic effects, it is proposed to develop a methodological framework that allows building mathematical 
models of vibration fields, taking into account the system of requirements for the features of interactions of system 
elements. The proposed methodological framework should have: 

— an effective mathematical apparatus; 

—a tool for visual display of structural features; 

— capabilities to build dynamic analogies to generalize the developed methods to a broad class of physical processes. 
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Materials and Methods. A structural approach based on a comparison of mechanical vibratory systems and automatic 
control systems can serve as a methodological framework [17]. In this case, the mechanical vibratory system is considered 
as a design scheme of a technical object. It is compared to the structural scheme of a dynamically equivalent automatic 
control system. The established correspondence makes it possible to analyze and evaluate the dynamic characteristics of 
technical objects, taking into account the ideas about the steady-state forms of vibrations of a mechanical system based 
on the apparatus of the theory of automatic control [18]. The technology of structural mathematical modeling has features 
related to the assessment of dynamic states. Transfer functions of structural schemes are used for it. The input signals for 
these functions are external disturbances, the output signals are vibrations in the coordinates of the object being evaluated. 
Each coordinate of the system can be correlated with an indicator of the amplitude of the steady-state vibration. This 
allows us to consider the distribution of amplitudes over the points of the system as a vibration field of a technical object. 
The construction of a vibration field, taking into account a number of requirements, belongs to the class of inverse 
problems of dynamics, which in general may be incorrect. One of the methods for solving ill-posed problems is related 
to the regularization of the initial problem by introducing optimality criteria. 

The listed features of the assessment of the dynamic states of technical objects indicate the need to develop a common 
approach in the methodology of constructing vibration fields, taking into account the criteria of optimality. At the same 
time, the methodology for the formation, evaluation, and correction of vibration fields with specified characteristics is 
not sufficiently detailed. Specifically, it is not possible to build vibration fields according to optimality criteria. 

The presented research is devoted to the development of the techniques of structural mathematical modeling. For this 
purpose, we solve a class of specific problems of constructing vibration fields of technical objects with conditions for 
vibration amplitudes that reflect the requirements for dynamic effects and take into account the criteria of optimality. 

Basic provisions. Formulation of the problem. A model mechanical vibratory system formed by a set of solids 
under conditions of vibrational loading of a kinematic origin is considered. The model system (Fig. 1) has 5 degrees of 
freedom, consists of 4 links in the form of solid-state rods with masses M;, moments of inertia J;, centers of gravity at 
point O;, spaced from extreme points A; and A;+; at distances /\; and /2;, where i = 1...4. 

Solid bodies in the form of rods are articulated at the extreme points — these are the attachment points of elastic 
elements mounted on support surfaces that perform in-phase harmonic vibrations (Fig. 1). 
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Fig. 1. Mechanical vibratory system formed by solids A1A2, A243, A3 A4, A4.As: A — articulation points of solids; 
O;— centers of gravity; zi — support surfaces; Ji — moments of inertia 


Under the influence of external kinematic disturbances z; from the support surfaces, the system performs small 
steady-state vibrations relative to the static equilibrium position. External kinematic disturbances are connected 
in-phase vibrations: 

Zi =ZegVi.- (1) 

Here, zg — harmonic disturbance of type Asin(wé), where @ — frequency and y;— connection coefficients of external 
disturbances; i = 1...5. In general, the connection coefficients (1) are external disturbances of a certain structure. 

The required dynamic effects are provided by a set of model conditions, namely: the vibration amplitudes of the 
coordinates of the control points of a technical object take fixed values regardless of the frequency of external disturbances. 

The task is to determine for a mechanical vibratory system such external kinematic disturbances that the corresponding 
vibration field at the control points satisfies a given set of conditions in the form of fixed vibration amplitudes. 

Research Results 

Mathematical model. To compile a system of differential equations, the coordinates {)1, v2, y3, ya, ys} of the 
displacements of point A;, i= 1...5 relative to the positions of static equilibrium are studied. The following coordinates 
are considered along with coordinates {y;}: 

— coordinates {Zo1, Zo2, Zo3, Zoa, Zos}, displaying the displacements of the centers of gravity O; of rods relative to the 
positions of static equilibrium; 

— coordinates {@1, 2, 3, @4, Ms} of small rotation angles around the centers of gravity. 
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Coordinates Zo;, i, yi are related by expressions: 
Zoi =i Vi + bi Via 7 
Qi =Ci( Visi — Yi)’ ( ) 
lin li qe 
lath’ = datla? da thin 


Here, aj = 


A mathematical model in the form of Lagrange differential equations of the 2nd kind is based on expressions for 
potential and kinematic energy. They can be represented in a matrix form using the scalar product (,) vectors in R”: 


1=5(K(3-2)( 3-2), T = >(Mbo,20)+>(J6.8), o 


here, Zo,@ — coordinate derivatives Zo, @. 


ky 0 Ji 0 M, 0 
k J M 
K= 2 Fe 2 Gyr 2 
k3 J3 M;3 
0 ka 0 J4 0 M, 
v1 Z| 
M1 ZO01 
y2 Z2 
= @2 = Z02 = S 
Q= »Z0= »>V=} V3 1,2 =] 23 
3 Z03 
v4 24 
4 Z04 
V5 Z5 


Vectors 7) ¢ R*, GE R*, PER? taking into account (2), are interconnected through the following relations: 
g g g 


Zo =Uy, @=Vy. (4) 
a bh 0 —c1 CC] 0 
b = 
Here, U = be ‘ V= ae 
a3 b3 —C3 C3 
0 a4 by 0 —C4 C4 


In coordinates {y;}, kinetic energy (3) takes the form: 
al T T = 5 
T=5((u MU+4+V IV)3,9). (5) 
Taking into account expressions (2)-(5), Lagrange equations of the 2nd kind take the form: 
(UT MU +V" JV) 9+ Ky = Kz. (6) 
It is assumed that the initial conditions are zero. In this case, the integral Laplace transformations lead (6) to a system 


of algebraic equations: 
(Mop? +K)y = Kz, (7) 


where My =U?MU+V" JV —matrix of the system, sign “~” above vectors y and z denotes the Laplace transformation [19]. 
Using the methods of [18], algebraic system (7) can be represented as a block diagram of a dynamically equivalent 
automatic control system (Fig. 2). 
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Fig. 2. Block diagram of mechanical vibratory system (Fig. 1): p=/®, j=v-1 


Transfer functions can be constructed based on a block diagram for a fixed set of connection coefficients y;. They are 
determined by the selection of an object, whose dynamic state is estimated by the expressions: 


Wil p) = 2, Wal p) = 22, Wal p) = 22, Walp) = 22, Ws pp HZ. (8) 


Zg Zg Zg Zg Zg 

Physically, we consider the amplitude-frequency characteristics of transfer functions (8) as lever connections 
representing the ratio of the amplitudes of the coordinate vibrations {y;} to the vibration amplitude of the external 
kinematic disturbance {zg}. 

Within the framework of the lever interpretation, the positive branches of the amplitude-frequency response display 
lever connections of the 1st kind (implemented by double-arm levers). Negative branches represent lever connections of the 
2nd kind (implemented by double-arm levers). The modes of resonance and dynamic damping of vibrations display degenerate 
lever connections, which can be represented by conditional virtual levers with zero or infinite arms. 

The selection of the structure of external kinematic disturbances through setting the connection coefficients {y;} 
uniquely determines the amplitudes of system coordinate vibrations. At the same time, the requirements imposed on the 
vibration amplitudes may: 

— either imply ambiguity of the solution; 

— either directly contradict each other; 

— or, due to large errors, allow only partial compliance. 

In these cases, criteria for the selection of external disturbances, which uniquely determine the vibration field 
according to the set of requirements, are needed. The logical approach for this task is based on the fact that external 
kinematic disturbances providing the required conditions can be determined by the criterion of optimality of vibration 
fields. The principle of this criterion is as follows: the vibration amplitudes of a technical object are selected in such a 
way that they deliver a minimum to a certain energy functional [20]. 

Formation of vibration fields of mechanical vibratory systems based on the optimality criterion. Within the framework 
of the model problem, it is assumed that a certain condition is set at points Ai, A... Aix of vibration field {y;}, i= 1...n of the 
mechanical vibratory system. It consists in the fact that amplitudes {yx}, A = 1...m take fixed values {yo1...Vom}: 


Vik = Yor. (9) 
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Here, 4 = 1...m, where m displays the number of points at which conditions for amplitudes are set; in this case, and 
further in the text, the sign “—” above y or z means the vibration amplitude of the steady component of the coordinate, 
unless otherwise specified. At the remaining points {y;, i= 1...n} / {yiz, k = 1...m}, the requirements for the values of the 
vibration field are not set. 

Needed requirements (9) can be presented in the form: 

Ay =b, (10) 
where A — matrix displaying the condition operator on the shape of the vibration field at the control points; b — vector 
of vibration field characteristics values [20]. 

In general, the conditions for vibration field (10) may be uncertain or incompatible. To rationally account for 
incompatible or uncertain requirements for the vibration field, a criterion of its optimality is introduced, as well as the 
assumption that the set of requirements can be fulfilled with some approximation. 

Values ¥1, ¥2...¥»n mean the vibration amplitudes of the steady-state components; therefore, for each fixed frequency 


@ of external disturbances z, it is possible to set a variational task of constructing an optimal vibration field, which: 
— delivers a minimum to the energy functional; 
— displays representations of vibration patterns or values of potential or kinetic energy. 
In general, the energy functional determines features of the distribution of the amplitudes of dynamic characteristics 
over the points of a mechanical vibratory system. 
Within the framework of the considered model problem, the energy functional specified by operator T., displays the 
nonuniformity of the vibration field. Its indicator is the sum of the squares of the differences in the deviations of the 
vibration amplitudes of nearly points: 


nel 
Te(F)= > (Fin -Hi}?. (11) 
i=l 
Zero value of energy functional (11) determines a uniform vibration field, whose vibration amplitudes of the points 
are equal to each other. Minimizing the energy functional, taking into account the specified conditions for the vibration 
amplitudes at the control points (10), determines the shape of the vibration field, extremely close to uniform. 
We make allowance for a number of conditions: 
— fixed frequency of external disturbances o; 
— vibration field is optimal in the terms of the criterion given by the energy functional; 
— amplitudes of the vibration field at the control points satisfy conditions (10). 
In this case, the formal formulation of the problem of constructing a vibration field can be represented as: 


Vo(@) = arg min TY), (12) 
Ay=b,¥=(—Mow? +K) Kz 


where y — vibration field defined by external kinematic disturbances Zz, which are considered as independent 


variational parameters. 
For fixed frequency o, vibration field y,.(@), which exactly satisfies the set of requirements, is formally found as 
solution (12) to the conditional minimization of the quadratic function (11) taking into account linear conditions (10): 


L(¥,A) = (Tey, ¥)+( Ay —b, A) > min ; (13) 


AF=(Mop?+K) | Kz 
Here, y — vibration field determined by external kinematic disturbance z ; A = (A1, ... Am) — Lagrange coefficients. 


Minimization in problem (13) proceeds according to the parameters Zz of external disturbances. For fixed values Zz of 
amplitudes of external disturbances, the created vibrational field y is the solution to the problem: 


(—Mow’ +K)¥ = Kz, (14) 
where My =U’MU+V7JV — matrix of the system. 
Vibration field y , which exactly satisfies the set values (10), can be determined under the solution to the conditional 
minimum problem using the Lagrange function: 
iz Aye (0 (Mop?+K) Kz,(Mop?+K)- kz)+ (4(Mop? +K) Kz—b,A) > min, (15) 
where Z — vector of vibration amplitudes of external kinematic disturbances; 7.9 — matrix defining the quadratic form 
of the energy functional (11). 
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External disturbances z, , forming the required vibrational field, the Lagrange multipliers A are the solution of the system: 


JG} as 


Thus, z;— solution to the problem of conditional minimization (15), (16). The optimal vibration field y, is determined 
by external disturbance z, and has the form: 


-l =i 


K?((-Moo?+K)") To (Mao? +x) K (4(-Mow? +x) K) 


A(-Mow? +K) 'K 0 


yo =(-Mow? +K) | Kzo. (17) 


Here, y, can be called an interpolating vibration field, since under the condition of compatibility (16), it exactly 
satisfies requirements (10). Taking into account the dependence on frequency w, the compatibility and certainty of 
system (16) is characterized within the framework of the theorems of the existence and uniqueness of solutions to systems 
of equations of the theory of splines [20]. Under the condition of compatibility and certainty of the system (16) and (17), 
the optimal vibration field y, exactly satisfies the set of conditions (10). 

Necessary conditions (10) may contain noise, errors, or be contradictory. Taking into account contradictory or 
unreliable requirements, the vibration field can be constructed by solving the smoothing problem. Its result is an optimal 
vibration field that satisfies conditions different from the initial requirements, but close to them. 

Vibration field assessment with account of noisy or unreliable dynamic performance requirements. Vibration 
field yo(w) (17) satisfying conditions (10) can be approximated using vibration field y,(m) taking into account the 
unreliability or inconsistency of the necessary requirements (10): 


yo(@)=arg min [ansseS/5, -mu)'} (18) 


¥=(-Mow?+K) Kz hel 
Here, 7.y — operator displaying the criterion of optimality of the shape of the vibration field; « — smoothing 
parameter indicating the priority of the criterion of optimality of the shape of the vibration field in relation to the 
requirements of conditions (10). 
In expression (18), Z is considered as an independent variable. Vibration field ya(@) can be called smoothing, since 


it approximately satisfies the necessary requirements of the system. 
Vibration field y, is determined by a system of equations obtained through minimizing the Lagrange function: 


La(¥(Z)) = &(Teo¥.¥) + AV Bl) > min. (19) 
Here, y¥(Z) is a vibration field formed by external kinematic disturbance Z: 
(—Mow? +K)y = Kz. (20) 
A replacement based on (20) is made: 
y= (Mow? +K) | Kz. (21) 


Now, the Lagrange function (19), taking into account (20) and (21), is reduced to the form: 


-1 


La(Z)= a (To (—Mow? +K) | Kz,(-Mow? +K) Kz) A(-Moa? +K) | Keo] > min, (22) 


where | = V(.-) — length of the vector in R”, expressed in terms of the scalar product (.,-) . 
Consider the external kinematic disturbances z,. They form optimal vibration field y,, which smooths out the specified 
conditions (10) and is a solution to the minimization problem (22). Such disturbances are determined by the system: 
-1 


Ke ((-Moo® + Kk)" ) (Teo + A” A)(—Mow? +K)" Kz, =K" ((-Moo® +K) ) A‘b. (23) 


After simplifications of system (23), the desired amplitudes of kinematic disturbances z, are determined by the 
expression: 

za =K~'(Mop?+K)(aTo+A™A) | Ab. (24) 

As shown above, external disturbances z, (24) are a solution to problem (23), taking into account the selection of the 


smoothing parameter a. These disturbances provide the construction of vibration field ya, which approximately satisfies 
the set of requirements (10), taking into account possible incompatibility or unreliability: 


Vo =(Mop?+K) | Kza. (25) 
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The type of vibration field y, (25) is clearly independent of the frequencies of external disturbances a: 
Vu =(aTeo + ATA) A. (26) 


External disturbance z,, which provides the construction of an optimal smoothing vibration field ya (24), can be used 
to construct an interpolation vibration field y, as a0. The selection of smoothing parameter o is determined from the 
conditions of smallness of the discrepancy: 

2 


(a) =|Ava -9| 


(27) 


where yo = (—Mow? + kK)" Kz, Za — kinematic disturbance (27). 


It should be noted that vibration field y, depends on frequency o implicitly, through the smoothing parameter a = a(@), 
whose selection is determined by the frequency of external disturbance @ and the value of discrepancy (27). 

Construction of a vibration field for a mechanical oscillatory system. Computational experiment. Let the 
parameters of the mechanical vibratory system (Fig. 1) be given by stiffness k;, masses Mi, moments of inertia J; and 
geometric characteristics /\;, /2; (Table 1). 


Table 1 
Parameters of the model mechanical vibratory system 
No. M,kg J, kg-m? k, N/m l,m b,m 
1 100 10 1,000 1 1 
2 200 20 2,000 1 1 
3 300 30 3,000 1 1 
4 400 40 4,000 1 1 
As requirements for the amplitudes of the vibration field, conditions are set for the values at points A; (Table 2). 
Table 2 
Vibration amplitudes at the points of the vibration field (y, m) 
Values A Ad A3 Ag As 
Required 0.002 ~ — - 0.005 
Smoothed out 0.002 0.0027 0.0035 0.0043 0.005 


Expression (11) is used as an energy functional, which can be represented using operator T-: 


|) (y2-H) (1-1 ee 

y2 a y2 

- ¥3—Y2 L. =f _ 

Te:| y3 || _  _ = ¥3 

_ y4-ys 1 -l _ 

4 aa aa 4 

: ys—ya 0 1 -l > 

Ms Vs 

Expression (T. eV le ey) can be represented as: 

1 -l 0\ yn yi 
= a y2 || ye 
(TeV, TeV) = (Te" TeV,V) = (TeV, ¥) = -1 2 -l Vs |.) Vs 
at 2 <1llgiil le 
0 -1 1 J ys) (ys 


Here, 7.9 — matrix that defines a quadratic shape. The value of operator 7.y is zero on a uniform vibration field 
y =const. The vibration amplitudes of the points of this field are equal to each other for an external kinematic disturbance 


at an arbitrary frequency. 
We set conditions (10) for the vibration field. For this purpose, data operator Ai; can be used, which maps the 
amplitude values at points A), As to vibration field y: 


Vi v1 
ale loft © 9 % Ws [_(% 
Pee. gy Oe ae hae 

ya ya 

Vs V5 
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The vibration field interpolating the set values at points 41, As (Table 2), can be approximated using the smoothing 
vibration field y, for smoothing parameter a = 0.01. We find the difference between the values of the vibration field and 
the required values for the selected smoothing coefficient a. It is determined by the residual function (27), whose value 
is approximately 0.00001 m, where b = (0.002, 0.005)’. The constructed vibration field, smoothing the interpolation 
conditions, has a linear shape (Table 2). For fixed smoothing parameter a = 0.01, the set of vibration amplitudes of the 
points of vibration field y, which provides proximity to the interpolated data with residual @(a) = 0.00001 m (Table 2), 
has the form of constant functions of the frequency of external disturbance @ (Fig. 3 a). 


ym (4) (5) z,m (5) 


10 «@, rad/s 
ie) 


—0.015 
0 2 4 6 @, rad/s 
a) b) 
ym 9(@), m 
(5) 
(4) 0.0003 
0.006 (3) 
Gd) \) 
0 1 2 3 4 5 0 2 6 @, rad/s 
n, generalized coordinate number 
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Fig. 3. Vibration field characteristics: a — values of vibration amplitudes of smoothing vibration field for frequencies of external 
disturbances in the range @ € (0.100) rad/s; 6 — amplitudes of external kinematic disturbances providing the required vibration 
amplitudes depending on the frequency of disturbances w; c — shape of the vibration field for frequencies in the range w€(0...100); 
(1), (5) — required values of the vibration field; (2), (3), (4) — values based on the energy approach; d— deviation @ (@) of the 
vibration field from the required amplitude values at fixed points of the technical object for constant smoothing coefficient a = 0.1 


In this case, the value of the vibration amplitudes of the generalized coordinates remains constant when the 
frequency of external kinematic disturbances z, changes. These disturbances, which are needed to provide the required 
amplitudes at fixed points of the vibration field, depend significantly on frequency of external disturbances (Fig. 3 b). 

We describe a vibration field for a fixed frequency or a frequency domain in which it remains unchanged. It can be 
represented by a graph of a function depending on the number of the generalized coordinate or the coordinate of the 
point of the system in which the vibration field is estimated. The shape of the vibration field is determined by the 
optimality criterion. Specifically, energy functional 7.0, which determines the optimality criterion, forms a linear shape 
of the vibration field (Fig. 3 c). It should be noted that the “interpolation” vibration field is constructed using a 
“smoothing” vibration field corresponding to smoothing parameter a. The deviation of the “smoothing” vibration field 
should be negligible. In the considered model example, for a = 0.1, the deviation of the smoothing vibration field is 
represented by constant @ (a) = 0.0001 m (Fig. 3 d). It should be noted that in the general case of the energy functional, 
the value of the smoothing coefficient a, which provides a fixed level of residual @(a) = 9", depends on frequency o 
and the system of requirements. 

Thus, the options for the formation of dynamic states of mechanical vibratory systems in the form of a vibration field 
satisfying a system of conditions depend on: 

— energy functional reflecting the characteristics of the vibration field; 

— operator displaying a set of requirements for the vibration field. 

Other options for the formation of a vibration field based on the selection of an energy operator are also possible. 
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Correction of optimality criteria for shaping vibration fields using the weight coefficients of energy operators. 
The shape of the vibration field with the possibility of equalizing the vibration amplitudes of the coordinates of a technical 
object is considered above. There are other approaches to minimizing the amplitudes of generalized coordinates, e.g., by 
selecting energy operator 7; in the form of an identity map: 

Tr: y > Ey. 3 
where £ — identity operator. 

The formation of a vibrational field based on the energy operator is associated with minimizing or zeroing the 
quadratic function of the amplitudes of in-phase harmonic vibrations of a fixed frequency. As a result of solving the 
smoothing problem, a vibration field is formed for the energy functional corresponding to unit operator Tz. The 
amplitudes of its vibrations at points A2, A3, Aq are zero (Fig. 4 a), i.e., at points A2, A3 and Ay, a dynamic vibration 
damping mode is implemented. 
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0.001 2 8) @ 
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Fig. 4. Possibilities of vibration field correction: a — shape of vibration field for energy functional displays minimum possible 
vibration amplitudes, taking into account the necessary requirements; 5 — amplitudes of external kinematic disturbances zi, 
providing the necessary requirements for dynamic features of vibration field; c —vibration amplitudes of vibration field of technical 
object depending on frequency of external disturbance; d — options of vibration fields depending 
on weight coefficient of energy operator 7p 


The totality of external kinematic disturbances is a solution to the problem of minimizing the energy functional, 
generally defined on the dynamic states of mechanical vibratory systems, taking into account the requirements for 
vibration amplitudes at fixed points of the vibration field (Fig. 4 b). 

The characteristic property of forced kinematic disturbances, which provide the required shape of the vibration 
field, is an increase in the vibration amplitude when the frequency of disturbances grows (Fig. 4 c). External kinematic 
disturbances z provide a fixed vibration field of a technical object at various frequencies w (Fig. 4d). Smoothing 
parameter a = 0.01 provides an absolute deviation from the required values at level ~(a) = 0.000053. This shows that 
the totality of the vibration amplitudes at the remaining points of the vibration field determines the characteristics of 
the energy operator based on the required values of the vibration field amplitudes at fixed points of the mechanical 
vibratory system. Specifically, operator 7. is aimed at forming a uniform vibration field, and operator Tz is aimed at 
zeroing the vibration field. 

Weight coefficients that display a combination of optimality criteria can be used. This allows combining the following: 

— requirements for convergence of vibration amplitudes of points of vibration fields; 

— requirements aimed at zeroing the vibration amplitudes of the vibration field points. 

Varying the weight coefficients that determine the combination of optimal criteria for the formation of a vibration 
field is actually a way to correct the vibration field. 
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We build a family of forms of vibration fields of a technical object that meet certain conditions. For this purpose, we use a 

combined optimality criterion implemented by a weighted sum of energy functionals, taking into account the weight factor: 
m5)" =|r-7]° +B|rex)’. (28) 

Here, 6 — weight factor that determines the combination of the optimality criterion based on the smallness of the 
amplitudes and the criterion based on the uniformity of the amplitudes of the vibration fields. The variation of the weight 
coefficient B leads to the construction of a family of vibration fields. The criterion of their optimality has a complex 
structure that takes into account the differences and absolute values of the vibration amplitudes of the points of the 
vibration field. The selection of a set of weight coefficients determines the appropriate forms of vibration fields, which 
are realized under the influence of external disturbances, which are the solution to the variational problem [20]. 

Based on the set of weight coefficients B, it is possible to construct optimal vibration fields that meet the requirements 
for vibration amplitudes at fixed points with sufficient accuracy (Table 3). 


Tey 


Table 3 
Forms of optimal vibration fields for weight coefficients of the energy functional 
1 1 0.01 | 0.0000841 10 0.0019704 | 0.0009850 | 0.0009845 | 0.0019686 | 0.0049213 
2 0.5 0.01 | 0.0000503 10 0.0019847 | 0.0014467 | 0.0016322 | 0.0026337 | 0.0049521 
3} 0.25 0.01 | 0.0000305 10 0.0019940 | 0.0018969 | 0.0022740 | 0.0032196 | 0.0049701 
4] 0.125 | 0.001 | 0.0000195 10 0.0020000 | 0.0022503 | 0.0027818 | 0.0036611 | 0.0049980 
5 | 0.0625 | 0.001 | 0.0000014 10 0.0020004 | 0.0024754 | 0.0031053 | 0.0039292 | 0.0049986 


Combining optimal criteria for the formation of vibration fields is a way to correct the dynamic states of a technical 
object (Fig. 4 d). Note that the found kinematic disturbances z;, considered as functions of frequency w, have one zero 


each (Fig. 5 /f). 
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Fig. 5. Optimal vibration field characteristics for weight coefficient B = 0.25 and smoothing parameter a = 0.01: 
a — ratio of the vibration amplitudes of points of optimal vibration field (va)1 to amplitude of kinematic disturbances (Za)1; 
b — ratio of amplitudes (ya)2 to (Za)2; c — ratio of amplitudes (ya)3 to (Za)3; d — ratio of amplitudes (va)4 to (Za)4; 
e — ratio of amplitudes (va)s to (Za)s; f— value of amplitudes of external disturbances (za) depending on frequency @ 
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Amplitude ratios y/z;, for which y; = y;(@) and z; = z;(@), are found on the basis of solution to the problem of 
determining optimal vibration fields, have one discontinuity of the 2nd kind (Fig. 5a—d). 

The specified vibration amplitudes of fixed points of the vibration field of the mechanical vibration system should be 
provided. For this purpose, combined energy functionals with account for the weight coefficients can be selected, and, 
based on these indicators, the required external kinematic effects with account for the variational principles can be 
determined. 

If the requirements for the vibration amplitudes of the vibration field are highly noisy or incompatible, it is possible 
to construct a vibration field that meets the requirements with a certain accuracy and is optimal for the combined criterion, 
but for other conditions. 

Determination of the optimal criterion for shaping vibration fields based on the residual function. In general, 
the requirements for a vibration field can be a complex set of joint or incompatible conditions reflecting limitations on 
the amplitudes of displacements, velocities, and dynamic reactions, including non-retentive conditions. 

Taking into account the criterion of optimality of the vibration field of a mechanical vibratory system, the set of joint 
conditions allows for a single solution in the form of a certain structure of external disturbances. In the case of 
incompatibility of conditions on the vibration field, the optimality criterion actually acts as a method of regularization of 
an incorrect inverse problem. 

The energy functional determines the specifics of the formation of vibration fields based on the values of the vibration 
amplitudes of displacements, velocities, force disturbances, or reactions occurring between the elements of the system. 
Energy functionals can convey ideas about the maximum potential or kinetic energy, or serve as a combination of them 
with account for weight coefficients that act as methods for correcting vibration fields. 

It is possible to add conditions according to the nature of the external kinematic or force disturbances themselves to 
the conditions for the values of the amplitudes of the vibration fields. They are expressed in the fact that the amplitudes 
of external disturbances take fixed absolute or relative values. The need to set requirements for input disturbances is due 
to the fact that to provide the specified forms of vibration fields, the amplitudes of external disturbances can increase 
indefinitely with growing frequency. 

Thus, here is a complex task of selecting: 

— an energy functional; 

— an operator of requirements for a vibration field; 

— an additional operator of requirements for external force or kinematic disturbances. 

Within the framework of the development of the structural theory of the optimal vibration field, it is possible to set a 
generalized task of its construction: 

Vo(@) = arg min Te(¥ ), 


Ay=b,3=(-Mow?+K) | Kz,Cz=f 


where C-Z = f — additional connection limiting the dimension of external kinematic disturbances. 


The technique for selecting an effective operator can be based on minimizing the function of additional residual, 


depending on the parameters of the energy operator: 
2 


Bo =arg min||Byp —c| 


where yg — solution to the problem of determining the vibration field with a combined optimality 
criterion 7, ((28) and Fig. 4) and conditions A (10); B — operator for determining values in a set of additional points; 
c — vector of values at additional points. 

This problem can be considered as the problem of selecting a “natural” optimality criterion [21]. Thus, the presented set of 
tasks and solution techniques can be the fundamental basis for the structural theory of the construction of vibration fields. 

Discussion and Conclusion. Elements of the structural theory of optimal vibration fields of mechanical vibratory 
systems have been developed. The main components of this theory are methods for constructing vibration fields of 
mechanical vibratory systems based on variational principles. The basic concept is the vibration field of a mechanical 
vibratory system. It represents the distribution of vibration amplitudes of generalized coordinates and other dynamic 
characteristics of a mechanical vibratory system over the points of an object with an estimated dynamic state. 

It is assumed that a system of conditions is set at a number of points of the mechanical vibratory system. It is expressed 
in terms of absolute or relative values of displacements, velocities, and reactions. Conditions are requirements for the 
characteristics of the vibration field of a mechanical vibratory system. There is also a system of conditions related to the 
features of external force or kinematic disturbances, and it also reflects the requirements for a set of absolute and relative 
values of the amplitudes of kinematic or force disturbances. 

For the development of the applied theory of optimal vibration fields within the framework of structural 
representations, an analogy can be established between the two pairs. The first one is formed from energy operator T and 
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condition operator A in the theory of abstract splines. The second one is formed by the criterion of optimality of the 
vibration field and a system of requirements for the characteristics of the vibration field at control points. 

Interpolation, smoothing and mixed splines are distinguished within the framework of the theory of variational spline 
approximation. Within the framework of the concepts of vibration fields of mechanical vibratory systems based on 
structural representations, the role of the spline is performed by a set of amplitudes of external force and kinematic 
disturbances, and the role of interpolation conditions is performed by a set of requirements for the amplitudes of input or 
output signals. 

It can be assumed that there are differences between the vibration fields of technical objects under conditions of 
external forced or kinematic disturbances, which are considered as a key factor in the formation of a dynamic state. At 
the same time, when constructing a vibration field, it is considered natural to estimate energy ratios depending on 
generalized coordinates. Their values are determined by an external force or kinematic disturbance. In fact, we need to 
find two vectors — generalized coordinates and external disturbances. 

We take into account the variational approach, as well as the set of requirements for the vibration field and external 
disturbances. We present the problem of determining external disturbances that form a set of dynamic states in the form 
of a vibration field with account for a system of conditions. In this interpretation, we are talking about the optimization 
problem in different formulations within the framework of variational principles. Specifically, it is possible to set tasks 
similar to the tasks of constructing a smoothing or mixed spline. It depends on the compatibility of the conditions for the 
input and output signals. We note the following an essential feature of the problems of constructing vibration fields: the 
frequency of an external disturbance is a parameter of the optimization problem, and in the case of systems with 
concentrated parameters, amplitude vectors serve as elements of an n-dimensional Cartesian space. 

Analogies between the tasks of approximation and construction of vibration fields provide developing techniques for 
the formation, assessment, and correction of dynamic states of technical objects. 

The improved structural theory of optimal vibration fields will be in demand in the aviation and space industries. Its 
apparatus can be used to maintain flight safety and durability of structures. In addition, it is logical to use this approach 
in mechanical engineering, electronics, and the development of medical equipment. The theory will allow solving the 
problems of forming, evaluating, and correcting the states of technical systems and devices under conditions of vibration 
loading. Another area of its practical application is the creation and improvement of measuring methods and means of 
complex vibratory systems. We also believe that mathematical models of optimal vibration fields are of interest for the 
development of new technologies and materials, in particular, materials with high damping, which can reduce the level 
of vibrations and noise in technical systems. 

As the main limitation of the technique, it is necessary to indicate that vibrations in the coordinates of an object, which 
are steady-state forms of harmonic vibrations or vibrations of a fixed amplitude, are considered. 
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Abstract 

Introduction. To describe the operation of temperature piezoceramic structures, the theory of thermoelectroelasticity is 
used, in which the mathematical model is formulated as a system of nonself-adjoint differential equations. The complexity 
of its integration in general leads to the study of problems in an unrelated formulation. This does not allow us to evaluate 
the effect of electroelastic fields on temperature. The literature does not present studies on these problems in a three- 
dimensional coupled formulation in which closed solutions would be constructed. At the same time, conducting such 
studies allows us to understand the interaction picture of mechanical, thermal and electric fields in a structure. To solve 
this problem, a new closed solution of a coupled problem for a piezoceramic round rigidly fixed plate has been constructed 
in this research. It provides for qualitative assessment of the cross impact of thermoelectroelastic fields in this 
electroelastic system. 

Materials and Methods. The object of the study is a piezoceramic plate. The case of unsteady temperature change on its 
upper front surface is considered, taking into account the convection heat exchange of the lower plane with the 
environment (boundary conditions of the Ist and 3rd kind). The electric field induced as a result of the thermal strain 
generation is fixed by connecting the electrodated surfaces to the measuring device. The thermoelectroelasticity problem 
includes the equations of equilibrium, electrostatics, and the unsteady hyperbolic heat equation. It is solved by the 
generalized method of finite biorthogonal transformation, which makes it possible to construct a closed solution of a 
nonself-adjoint system of equations. 

Results. A new closed solution of the coupled axisymmetric thermoelectroelasticity problem for a round plate made of 
piezoceramic material was constructed. 

Discussion and Conclusion. The obtained solution to the initial boundary value problem made it possible to determine 
the temperature, electric and elastic fields induced in a piezoceramic element under arbitrary temperature axisymmetric 
external action. The calculations performed provided determining the dimensions of solid electrodes, which made it 
possible to increase the functionality of piezoceramic transducers. Numerical analysis of the results enabled us to identify 
new connections between the nature of external temperature action, the deformation process, and the value of the electric 
field in a piezoceramic structure. This can validate a proper program of experiments under their designing and 
significantly reduce the volume of field studies. 


Keywords: thermoelectroelasticity problem, coupled problem, round piezoceramic rigidly fixed plate, biorthogonal finite 
integral transformations 
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CBxa3aHHad OCeCHMMeTpHYHAA 3a/1ada TEPMOIIEKTPOYNPYLocirH AIA KpyrJion *KeCTKO 
3akpellIeHHou M1acTHHbI 


JIA. Masxnn®, E.B. Casunora@X 
Camapcxkuii rocyapcTBeHHbIli TeXHHYeCKHH yHuBepcnuter, r. Camapa, Poccuiickaa Deyepayna 
slenax@yandex.ru 


AHHOTalna 

Beedenue. J[na onucanua paOoTbI TeMMepaTypHbIX Ibe3OKepaMH4eCKHX KOHCTPyKUM UCMONb3yeTcA TeOpHA 
TEPMOIICKTPOYMpyrocTH, B KOTOPOM MaTemMaTHYecKad MOeIb CPOPMYIHpPOBaHa B BUC CHCTEMbI HECaMOCONIPs.KEHHBIX 
TuddepenuMarbHblx ypaBHeHHit. Ciox%KHOCTb ee MHTErpHpOBaHHA B OOMIEM BUe IPHBOAUT K HCCeOBAaHHIO 3aay B 
HeCBA3aHHOM MOCTaHOBKe. ITO He MO3BOIACT OLCHHTh 3PPeKT BIIMAHUA ZICKTPOYMpyrux Moe Ha TeMMepaTypHoe. 
B mmtepatype He lIpeycTaBIeHbI HCCIeqOBaHHA JaHHBIX 3aa4 B TPeXMepHOM CBA3AHHOM MOCTaHOBKe, B KOTOPBIX OBLIH 
ObI MOCTpOeHbI 3AaMKHyTHIe pemtenua. [pu 3TOM MpoBeyeHHe MMCHHO TaKHX MCCIeAOBaHH MO3BONACT MOHAT KapTHHy 
B3aMMOJCHCTBUA MCXaHHYeCKHX, TeMJIOBbIX H 3JICKTPHYeCKHX MOJel B KOHCTpyKuMH. JIA pewleHHA TaHHOH MpoOseMbBI 
B HacTOAMelH padoTe MOcTpoeHO HOBOe 3aMKHYTOe PelIeHHe CBA3AaHHOM 3aa4H JIA Mbe30KepaMHMyeckOH Kpyriol 
%KECTKO 3aKpelJIeCHHOH TWIaCTHHbI, MO3BOJIMIOWIee KAaYeCTBCHHO OLWCHHTb B3aHMHOe BIIMAHHe TepMO3IeKTPOypyrux 
nose B TaHHOM dIeKTpoylpyroi cucTeme. 

Mamepuaaoi u memoovi. OObexToM UccieqOBaHMA ABJIAeTCA Mbe3OKepaMuyeckasd MlacTuHa. PaccMatpuBaetca Cryyali 
HeCTalMOHapHOTO H3MeHeHHA TeMIlepaTypbI Ha ee BepxHel JIMeBOM MOBepXHOCTH MpH yueTe KOHBeEKIMOHHOrO 
TeMOOOMeHa HWKHeM IWJIOCKOCTH Cc OKpy2KaloleH cpeyow (rpaHuyHbie ycuosua | u 3 poga). Unayuupyemoe B 
pe3ybTaTe oOpa30BaHHaA TeMMepaTypHbIX AeopMalHi seKTpuyeckoe Moe PuKcHpyeTcaA MyTeM NOAKTOUeHHA 
3JIEKTPOAMpOBaHHBIX TOBepXHOCTeli K H3MepHTeIbHOMy UpHOopy. 3aqaya TepMO3IeKTpOyMpyrocTH BKIIOUAeT 
ypaBHeHHA PaBHOBeCHA, IICKTPOCTaTHKM MH HecTalMoHapHoe rulepOomMuecKoe ypaBHeHHe TerIonpoBoyHocTu. Oua 
pewiaetca oO0ON[eHHBIM MeTOJOM KOHeYHOrO OMOpTOrOHabHOrO peoOpa30BaHHA, TMOSBOJAIOWerO MOCTPOHTb 
3AMKHYTOe pellleHHe HeCaMOCONps.KeHHOM CHCTeEMBI ypaBHeHHi. 

Pesynbmamet ucciedoeanua. IloctpoeHo HOBOe 3aMKHyTOe pellieHve CBaA3aHHOM OCeCHMMeTpPHYHOM 3aa4qH 
TEPMOIIEKTPOYNPYTOCTH AIA KPyTJIOH MWAaCTHHBI, BbITIOJIHCHHOM U3 Tbe30KepaMH4ecKOrO MaTepualia. 

Ooécyocdenue u_ 3aksjo4uenue. TlomyayeHHoe pellieHve HavasIbHO-KpaeBOH 3aauH NO3BONIMeT OlpeeuMTb 
TeMUlepaTypHoe, dIeKTpH4eckoe UM yipyroe NOsA, HHAYUMpYyeMBle B IIbe30KepaMHYECKOM IJICEMEHTE Tp MpOW3BOJIBHOM 
TeMUepaTypHOM OCeECHMMeTPH4HOM BHELIHeM BO3JelcTBHH. ITpoBeqeHHble pacdeTbI MO3BOJIAIOT Ope eIMTb pa3sMepsl 
CIJIOUIHBIX 3JIEKTPOAOB, KOTOPbIe alOT BO3MOXKHOCTL MOBbICHTL (PYHKIMOHAJIbHbIC BO3MOXKHOCTH Ibe30OKepaMH4eCcKHX 
TpeoOpasoBatenet. UnceHHbIM aHasIM3 Pe3yIbTATOB MO3BOJIACT BbIABHTb HOBbIC CBA3H Me%Ky XAapaKTepOM BHELIHero 
TeMIIepaTypPHOro BO3ACHCTBHA, IIPOeccoM AesOpMUpoOBaHHA M BEIMYHMHON IICKTPHYECKOTO MOA B Ibe30KepaMuyeckoni 
KOHCTpyKUMH. SOTO aeT BO3MO%KHOCTb OOOCHOBaTb pallMOHaIbHy!O MporpaMMy 9KCIIepHMeHTOB pH Ux 
IIpOCKTHPOBaHHH HM 3HAYHTEIIbHO COKpaTHTb OOBEM HaTypHBIX UCCIeLOBaHHH. 


Korouesble cJ10Ba: 3aqada TeCPMOSJICKTpOYHUpyrocTH, CBA3aHHaA OCCCHMMeCTpH4HaA 3ajtada, 2HeCCTKO 3akpeIvIeHHaAd 
TJlacTHHa, OvMOpTOrOHasIBHBIe KOHCUHbIC HHTCIpaJIBHbIle IIpeoOpa30BaHH4a 


Buaarogapuoctu. ABTopbl BbIPaxKaloT OaroqapHOCTb pelleH3eHTaM 3a TIpoBeseHHy!IO paboTy, KOTOpad MO3BOMIAa 
MOBbICHTb KAYeCTBCHHBIM YPOBCHb CTATbH. 


Aaa uutTupospanna. Wsaxuu J].A., Casnunosa E.B. Caa3anHaa ocecHMMeTpH4uHad 3aada TepMOIIEKTPOYIIpyrocTH WA 
KpyrIoH %KeCTKO 3aKpellieHHow macTuHbt. Advanced Engineering Research (Rostov-on-Don). 2024;24(1):23-35. 
https://doi.org/10.23947/2687-1653-2024-24- 1-23-35 


Introduction. Various mathematical models are used to improve the functionality of piezoceramic sensors [1-3] 
based on the interdependence of thermoelectroelastic fields. To more accurately account for the effect of coupling of these 
fields, it is needed to construct closed solutions. Some simplifications are used to solve systems of initial nonself-adjoint 
differential equations. Thus, the problems can be considered in an uncoupled formulation, or the problems consider and 
analyze elements that have a degenerate geometry. An uncoupled stationary problem for a long electroelastic cylinder is 
considered in [4, 5], and article [6] is devoted to the analysis of thermal stresses in a hollow sphere. Papers [7, 8] are 
related to the determination of the temperature field in a piezoceramic shell and a round plate in solving uncoupled 
problems. Coupled dynamic problems for a homogeneous piezoceramic layer, as well as dynamic problems in a coupled 
formulation for a gradient-inhomogeneous piezoceramic layer, were considered in [9 10]. In [11, 12], fields in an 
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unbounded medium were analyzed. In [13, 14], a long hollow cylinder was considered, and thermoelectroelastic fields 
were analyzed. 

Currently, the literature does not describe the results of constructing closed solutions to the mentioned non-stationary 
problems in a three-dimensional coupled formulation. Therefore, in this paper, we consider a round plate made of 
piezoceramic composition and having a rigid fixation, for which a new closed solution to the problem of 
thermoelectroelasticity is obtained. The use of a limit on the rate of temperature change on its front surface [10] makes it 
possible not to include the inertial characteristics of the system under study and apply the equilibrium equations in the 
calculated ratios. 

Materials and Methods. In the process of solving, a generalized finite biorthogonal transformation was used, which 
provided the reduction of the dimension of a nonself-adjoint system of equations and the construction of a closed solution 
through significant simplifying research in the image space. 


Mathematical model. Consider certain area Q: {0<7r, <b,0<0<22,0<z, <h"}, which is occupied by a 


piezoceramic solid circular plate in the cylindrical coordinate system (r,,0,z..) Arbitrary temperature boundary conditions 
can be used for the problem under study. However, for the certainty of the solution, on the upper (z, = 0) front surface, 
the temperature change j(7,,t,) at a given ambient temperature 9° on the lower (z,=h") 


plane (t. — time) is considered. The cylindrical thermally insulated surface is rigidly fixed: there is no radial component 


of the displacement vector and the angle of rotation, and its lower part is fixed in the vertical plane. The lower plane of 
the round plate in question is grounded. The front electrodated planes of the plate are connected to the measuring device. 
The design scheme of the plate is shown in Figure 1. 


eZ 4 


2x 


Fig. 1. Plate design diagram 


The mathematical formulation of the problem under consideration in a dimensionless form for an axially polarized 
piezoceramic material with a hexagonal crystal lattice of 6 mm composition has the form: 
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6O 
o)2-0 = 0, oi--1 a 0, ©),-0 =@), [+10 | = 1795 
Oz |z=h 
0VU,W, . 
es -0,2 -6,; (4) 
Ot \'=0 Ot \=0 


t=0,{U,W,, @}=0, 


where 
* * €31 * hk e31 ek 
{U,W,r,z} ={U iW teste) Tb, {do} = {6 5}, {6.60} =<—H{o oo}, 
és x GC C3 4+C. es te 
{0,0,,9}=t{o (of -7),(9 -T))\,a, = 8 pa es UTS pcg MI Ea, 
Cy Cy Cy e31 
_ Cx _ eis _ 33 _ 133 _ €33 _ 231 (15 +e31) _ &15€31 
, ’ é L x 2 ) 
a4 as a6 ay ag ay a9 
Ci e31 e31 Yu e11 Cy E11 Cie 
€33 811831 Y11¥ 33 83311 
a1; = 49 »Qj2 = , a4 =Ty ,a5 =Th » 416 =a, —a), 
C15 Yu81 Ck e3k 
a, =a-b/A, O*(r,,2.,t.) — temperature increment in dimensional form; U(r, 2.,t.),W "(te Ze te) 


o* (7, Ze, ts ) — components of the displacement vector, electric field potential; o-, (7, z, f), Oz (”, z, ) — components 

of the mechanical stress tensor; D,(r,z,#)— radial component of the electric field induction vector; A, k, a;— 

coefficients of thermal conductivity, volumetric heat capacity, and linear thermal expansion; 49 tr t) — electric 

potential induced on the upper front surface; yi, gi — components of the tensor of temperature stresses and pyroelectric 

coefficients (i= 1,3, yu = Cia); e15, €31, €33, €11, €33 — piezoelectric modules and permittivity coefficients; ©” = T— T; 

T, To) — current temperature and temperature of the original state of the body; B,.; — relaxation time; a — heat transfer 
1 


coefficient, ©)» — known rate of temperature change; V = ap =; 
ror 


To determine the potential of the electric field induced under deformation on the upper front surface 9 (r, ae in the 


case of connecting electrodes to a measuring device with a large input resistance, an additional boundary condition is used: 


0 

— | D,,-odS =0 

= [ Poe (5) 
(S) 


where D. (r, z, t) — axial component of the induction vector; S — surface area. 
Construction of a general solution. To fulfill the condition of fixing the cylindrical surface of the plate in the vertical 
plane, new functions w (7, z, t), Wi(t) are introduced: 
W(r, Zs t)=Wi(t)+w(r, Z, t), (6) 
this makes it possible to form a boundary value problem with respect to functions U, w, o, ©, which is investigated by the 


method of finite Fourier-Bessel transformations: 
1 


Uy (1,201) = [U(r.208) 0 (Gyr dar (7) 


0 


{wi (n,z,t),dn (n,z,t),Nu (n,z,t)} = [ 6-(0.2.1).0(r.2.1),.(r.2.0) rJ (inr)dr, 


ee ay ee Ga (8) 
{w(r,z,t),0(r,2,t),O(r,z,t)} = 25° {wi (n,z,t) dy (n,z,t),Nu (m2.t)} (in), 


n=0 Jo GP y 


where /, — positive zeros of the function J; (in ) (n = 0,0; Jo= 0), J,(...) — Bessel functions. 

It should be noted here that to satisfy the last boundary condition (2), it is necessary to assume that the original 
temperature of plate Tp is equal to the ambient temperature 9”, and the temperature increment function on the upper front 
surface (1, ¢) = 0. These assumptions, without much error, allow assuming that on the cylindrical surface of the plate, 
OC, z, ) =0. 
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As a result of using the transformation algorithm in the image area, the following initial boundary value problem is 


obtained: 
Ou Ow ob 
—j7uy + i i, —t +a3j, —-+i,Ny =0; 9 
Inn Tay az? A2Jn az 43) n az InN (9) 
Ow Ou eo ON 
=f Weta, Fee a be td #_Q, 
A JnWy + a4 an? a2J a Asindy +A az? ay a 
; . Ou . ow ; ON 
indy — 4s ae as ae AiodnWa +4 ap? ada + a3 a 
0?N 0 ig Ow oo 
i> Ny + us + Nyt i Uy +— |-a H | — 0; 
JnNu ap2 at baa HT 4} JU az ae 
Z=0,h; aig jy + a4 Oni + d¢6 On —a,Ny = 0, Mn JnWu = 9, (10) 
Oz Iz Oz 
ON 
O r\2-0 = 90H» Vx zn =9, Naz-0 = on Sts anNy | = 4179 15 
|z=h 
OWn> 
t= 0,{u7.04.Na} =0, Wy = Won (nds 9 (11) 
Ot \t=0 
Owy ; ON . 
= Won» = Non; 
Gy GE py 
I 
. dW, (t ‘ 
where [or aun Non} = |, 8. (0), 2) 18 joa 
t= 


0 
At the next stage of the solution, the introduction of functions U;(n, z, t), Wx(n, z, t), x(n, Z, t), Ox(n, Z, t) using the 


following relations: 
UH (n,z,t) =H (n,z,t)+Un (n,z,t), WH (n,z,t) = A, (n,z,t)+Wu (n,z,t), (12) 


Ou (n,z,t) =H; (n,z,t)+@u (n,z,t), Nu (n, z,t) Ay (n,z,t)+ On (n,2,t), 
allows the reduction of conditions (9) to homogeneous. 
Here, {H)...H4}={Hi...Hi\+{ fo(z)...fi2(z)} bou(t), {Hi...Ha} ={fi(z)...fa(z)pour(t)+{ fs (z).-.fa(z)} 9a, 
fi(z)...fi2(z) — twice differentiable functions. 
Substitution (12) in (9) — (11) when the conditions are satisfied: 
OH 0H3 


oH 
z=0,h3 ais jnHi +a4 +6 a7H4 =0; jn»H. -—+ =0, (13) 
Oz Oz Oz 


OH. 
F13)2=0 = on, F3)2=h =0, 42-0 =O1H; ES ants | =a79n, 
IZ 


|z=h 


provides the formulation of the following task: 


U OW, . O : 
= j2U 1 +) a jn + 3 jn + jn = Fi, (14) 
Oz Oz Oz 
. OWn . Uy j 0°OH OOu 
2 2 = 
—ajnaWy +a4 Ae? + 42jn As jn Ont + as cae Fy, 


. fo _ OU ; OW, ; 
jaOu —ag ue +daojn = ao jaWn tan _ + a2 jnOu +a3 
Oz Oz Oz 


fol a) o? : OW, 0 
—jrOn + Qu -( +B || av jn + wi) a oH | ra 
z z 


OQxu 


=F, 


Oz? ot Ot? 
Z> 0,h; a6 JnU 1 +a4 OW + a6 Opn a7Ou = 0; ea TiWay = 0, (15) 
Oz Oz Oz 


fa] 
On|z=0 = Ou|z=n = 9, Onz-0 = 9, (2H anQu | = 0; 
IZ 


|z=h 
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t=0;UH =Uon; WH =Won; Ou = Gon; Qu =Qon, 


OU rn OWn c On : 0On 


=Uon, =Won, =Qou, =Oon; 
Ot \=0 Ot |r=0 


Ot \=0 Ot \r=0 


where 


7M, OH 0H; 
Fi =j7H\-a +42 Jn a3Jn int 4, 
1= Jnl, -a 352 2J es 3] a Jnfl4 


07H: OH 4 
Oz? Oz 


O°H oH 
Fy =a jrHr a4 e a2 jn at asin a6 
IZ 


Oz? 


. o? _ OH 0?H oH 
FP; =—j,H3 + ag ai aon 10 fn Ha — a1 = —ay2 jn 4 — a3 a 
074 6) fol OH 0H; 
22 2 
Fy = jnHa- ae (20% AHa4t+ai4 ide — ais a |’ 
Von =—Hy-0, Won = Wow — A x-0, Pon =—H3p-0, Qon =—Hay-0, 
: OH, : : 0H . 0H; . : OH 4 
Uox =-— Won = Woe -——__— , QoH =-—___ , Oonh = Now - ‘ 
Ot |=0 Ot \=0 Ot \=0 Ot |t=0 


(16) 


Using the biorthogonal finite transformation (CMD) [15], we obtain a solution to problem (14) — (16). CMD with 


unknown components of vector functions of transformations is introduced on _ the 
KiQ\in, z). ‘ Ka(Kin, Z), N\ (Lin, Zz). : -Ni (tin, z): 


h 
G(n,hin,t) = jlo + a4 Gz ot — na we Jeane 
IZ 
0 


IZ 
{Mi (HWin,z),N2 (Win,z),.N3 (Win,z),Na (in.z)} 
[Kin 


{Un Wi. Ou,Ou}= YG n Rint) 


i=l 


Kin 


h 
; = | Ks (hineZ)Na( iz Me 
0 


2 5) 


segment [0, /] 


(17) 


where Ain, Lin — the eigenvalues of the corresponding problems with respect to the components of the vector functions of 


the CMD (k= 1...4). 
In the course of transformations, we obtain a task for determining transformants G(n, Ain, £): 


2 —— 
(poo S 98 Jo(mdat =-Fi(n,hint),(i=10 1 =0,00) 
d Kin, , 
1-0, 6Gin.n.0)= Gye, Oe) <esbigen, 
dt \t=0 


whose solution has the following form: 
G(n,hin,t) = (muin —M2in - {(Go — Gomain ) exp (mint ) - (Go — Gomiin ) exp (maint ) +, 


ri pF (1, hin, t)] exp (main (t-1))—exp(min (t —t))]dt , 
0 


in addition, two homogeneous problems with respect to components Ki(Ain, Z)...Ka(Ain, Z), 


d*K in 7 dK in . dK in : 
: a 2 ag Jn : +2 ais jnK sin =0, 
Zz 


—j7Kiin +a af 2Jn 
d*Ko; dK; d?K3; dK 4; 
22 2in . lin +2 3in 2 4in 
a jnK rin + a4 ae +42Jn a a0 jnK3in +a FE: Mina Zz 0, 


d°*K3; dK; d*K >i dK 4j 
A) in ‘ in 2 in 2 in 
Jn K3in as dz a3Jn re as jnK2In + d6 ie +h5,a15 dz 0, 
a°K in ° dK in ° dK in 
(13, — ja) Kain + F = + jinKiin +7 = 0G {Rag tig = 0: 
IZ a 


(18) 


(19) 


(20) 


(21) 
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dK in dK in 
Z= 0, h, a6 JnK iin + a4 2 +a) 3 V2 aisKain = = 0, K3in\z=0 — K3in\z=h = 0, (22) 
dz dz 
dK in : . dK in 
dz ay dz = 


and Ni (tin, Z)...-Na( Lin, Z): 


d?Niin a 5 dN2in | ‘ AN 3in 
ge ee ge 


d °Noin dNiin 2 d7 Nain AN ain 
+ a2 jn ——— a in N3in +a a 
dz? ag dz eo ° dz? : dz 


s d? N in dN in ‘. d?N in " dN. in 
TaN sin ag = + dojn : aro jnNoin + ay = + a2 jnNain +413 <n = 0, 
dz” dz Iz dz 


—f2 Min + ay + jnNain = 9, (23) 


=0, 


2 
—a jn N2in + a4 


d’ in dN in dN in 
—j2Nain + al +ui,| Nain + ai4jnNiin + aia 2 =e )- 0; 
dz dz 
dN in dN in ON in : 
Z=0,A, aie jnNiin +44 zZ + a6 z a7 Nain = 90, — JaN in =9, (24) 
dz dz fe) 


ON sin + aNen| _ 0; 
Oz 


|z=h 


N3injz=0 = N3injz=n = 0, Nainiz-0 = 0, ( 
where 


h 
(n, Met) =[( FiKiin + F2K2in + F3K3in + F4K 4in dz, 
0 


h Peel a5 d {pou ou} IK 4indz, 


(G01. on} =| ae {Uon Uo} + i iE 
Min, M2in — roots of the characteristic equation: Bm?, + mj, +7, = 0. 

Constructed homogeneous problem (23), (24) with respect i functions Nj (tin, Z)...Na4([in, Z) is invariant to the initial 
calculated relations (14), (15). 

Systems (21), (23) are reduced to the following equations with respect to Ko(Ain, z), N2(bin, Z): 


a’ 6 d 4 d 2 
(S + @lin dz + €2in dz4 + €3in iE + €4in oa 7 Noin } = 0. (25) 
In the paper, coefficients e1j,...e4in are not given due to the limitation of its volume. 

In equation (25), the left part is decomposed into commutative factors, presented below: 


2 2 4 
(5-4 (S45, [Soom < + | an Na) = 0, (26) 
Zz ZL ZL 


e4i _ : 
where Alin = V Buin, Ain = ~VStin » M5, = Ctin + Bin + Shins Main = a a » Bin, S1in — real positive roots of the following 
2in 


ee . lin 
characteristic equations: 
4 3 2 
Bin at e1inBin + €2in Bj, + €3in Bin + Cain = 0, 


C4in 
So = (tin ag Bin Si te (e1inBuin ot he + €2in )Sin = B. a 0. 


lin 
When examining a round rigidly fixed piezoceramic plate, the general integral of equations (26) has the following form: 
{Kain , Noin } = {Din Ein \ exp ( Atinz) a {Doin , Edin } exp (—Atinz) + {Dsin EL 3in \ sin (Aainz) at (27) 


+{ Dain Lain} cos (Arinz) ot {Dsin, Esin \ sin(A3inz) ate { Dein, Eoin} cos ( A3inZ) tr 
+{Drin, Erin} sin (Asinz) + {Dsin, Esin } cos (AginZ), 


where 
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0.5 0.5 
2 ly 4 2 2 4 2 
M3in + M3 in —4m4,, A _ M3 in — YN 3in — 4mj;, 
Asin = > A4in . 


2 


It should be noted here that the condition of the actual positive values of coefficients Byin, Siin, Aiin..-A4in iS fulfilled 
for most structures made of piezoceramic material. Otherwise, the formula structure (26), (27) simply changes. 

Considering that the connections were previously obtained as a result of reducing (21), (23) to (25), we get expressions 
for functions K\(Ain, z), K3(Ain, Z), Ka(Ain, Z), Ni(Ains Z)s N3(Ains Z), Na(Ain, Z)- 

Substituting Ki(Ain, z)...Ka(Ain, Z), Ni(tin, Z)...Ni([in, Z) In conditions (22), (24) provides determining constants 
Diin...Dsin, Ein ...E'in and eigenvalues Ain, Hin. 

The final expressions of functions U (n, z, f), W (n, z, f), 0 (n, z, 1), © (n, z, f) are obtained by applying the inversion 
formulas (17), (8). Then, taking into account (6), (12), we have: 


(%, z,t) ay Jil (ine an n, z,t) + Grn, Din t )N1 ( Hin, Z \ Kin || 
n=l Jo( (jn) i=l 


| (28) 


W (r,z,t) =Wi(t my) #0 


ue 1 (n,z,t) +L Omran Nol bin,z)||Kin | | 
jn) = 

) 
2 


(7, Z, t) )=2)) an fn n 12,1) +) G( 2, dinst JN3 (meal | 


(r,z,t) ye 0 (Jnr | (n,z,t) + Grn, Din, t )Na ( Hin»Z )|Kin| | 


+ Jo( (jn) i=l 


Functions /i(z)... fi2(z) are calculated from the simplification condition F\... F4 when conditions (13) are satisfied: 


7H, ee) 0H; 
29) , , : _ 
Jnl a az? + a2Jn ae a3Jn Be JnHl4 =0, 


2H 


; . 0H F 
aijnHr-a4 3 a2 jn + as jn Hs — ao 7 ta Or =0, 
oH OH : OH . OH. 
— jn Hs + ag —— an? : — ao Jn jt ao jn He — a a2 Aina — an =O 
0°H, 
23) _ 
jJnHa - az2 = 0. 


Function W;(f) is determined from condition W(1, h, t) = 0: 


t)=-2))] Ha (ht) + > G(n, dint )No (Wink) Kill Yo (in) 


n=0 i=l 
For a qualitative assessment of the induced electric pulse on its upper front surface, it is required to form two electrodes 
with a radius of separation R and connect them to a measuring device. In this case, potential do(r, , induced on two 
equipotential surfaces is represented as: 


b0 (1,1) = 01 (t)H (R-r) +600 (t)H (rR), (29) 


where H(...) — the Heaviside step function. 
Substituting (29) into (5) makes it possible to define expressions for determining potentials o(), o2(): 


R 1 
| Dap-otdr = | Dap-otdr =0. (30) 
0 R 


As a result of solution (30), functions do(f), bo2(¢) are defined as follows: 


_ p2 
on (0) = O51 [01(0)+02(+ Os (0) 4 (0)=2e!| -0.()+ 0. 0)+04(0), 
where 
Ji (jnR) { aro ai OH2*(n,z,t) as OH3"(n,z,t) 
=R) —— >, —A4A1'(n,0, pa Ls ae 
ae n=l oe a ar Oz je=0 Jn Oz \2=0 
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yet 2,0) ans) Jae) } 
0 


i=l ||Kin||° as Jn dz iz=0 Jn dz |z= 
0H 2" (0,z,t) 0H3"(0,z,t) ies G(Ai0,0,t) dK» (Ai0,z) dK3(Xi0,z) 
jak . 
2 ( ) a Oz \z=0 a Oz \z=0 2 Kio |’ . dz |z=0 ae dz \z=0 


Q3(t)= anf or (r,t) rdr, O4(t) = anf or (r,t)rdr, 
0 R 


on [En et oF) 20] AU] | 0 pepe tele ati) I 


dz |z=0 Jo (i, ) a5 Jn dz \z=0 Jn dz 


n=1 


Qo2 2 Jo (hr ) as Jn dz \z=0 Jn dz 


In this case, the potential difference V(#) is determined by the equality: 
V(t) = b01 (¢)— 02 (7). (31) 
Research Results. Numerical results are presented for a plate made of piezoceramics of the composition 
PZT-4 [4, 11, 16]: 
{Cu, Cra, Cis, C33, Css} = {13.9, 7.78, 7.3, 11.5, 2.26}<10"° Pa, {e11, 633} = {6.46, 5.62} x 10° F/m, 
{é15, 31, 33} = £1257, —5.2, 15.1} Ci’, A=1,6 W/(m-:K), Ou = 0.4«10> KI, 
k= 3x10° J/ (m*-K), gi = 933 =-0.6x104 C/(m?-K), Bret = 10+ 5, a = 5.6 W/(m?-K). 


The following case of temperature change j (7+,t«) is investigated: 


0} (T,t«) -(1-4 |r sn( tH (ta —te)+H (ts — thax ) ; 


where Tinax, tmax —- Maximum temperature value and the corresponding time (Gas =100°C, To = 20°C). 


dz z=0 a z= n=l 


Figure 2 shows graphs reflecting at various points in time (thax =01 s) the change in temperature ©*(0, z, f) in the 
thickness of the plate (b = 14x10 m, h* = 1x10 m). 

According to the calculation result, it is observed that due to the high coefficient of thermal conductivity and the small 
thickness of the piezoceramic plate, the steady-state temperature regime is formed quite quickly (tna =10 s) when it 
reaches ©*(0, z, #) on the lower front surface (z = h) 78°C (Fig. 2). 

Figure 3 shows the change ©* (0, Up ; t) in time (toe =3x10° s) taking into account (represented by a solid line) 


and without account for (represented by dotted line, B = 0) the relaxation of the heat flux (b = 14x10> m, h” = 1x10 m). 
It should be emphasized that the application of the hyperbolic Lord-Shulman heat conduction equation is needed only in 


the study of a piezoceramic micro-dimensional structure with a very rapid change @/ ( Vs,te ). 


©" (0, z,t),°C 


60 


40 


20 


0 0.02 0.04 0.06 Z 


(rN Hol 7, Hnl?) Jo2ey[ 2 Yo fo(z)+ diol) _as afin(z) I 
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Fig. 2. Diagrams ©*(0,z,t)—z:1—t=tmax, 2—t =10tnax, 3—t =100thar 


0 2x10+ 4x10+ 6x104 t 
Fig. 3. Diagrams o°(0 WA t)-t 1-44 : 
g. 3. Diag We pit | 
solid line — B =10~°(s), dotted line — B =0 


The numerical results of determining function ©°(r, z, t) show that when conducting a study of a structure made of 
piezoceramic material, it is possible to neglect the impact of the rate of change of body volume and tension on the 
temperature field, i.e., to use only the equation of thermal conductivity in calculations. 

Figure 4 shows a diagram of movements W*(0, z, f) over time ¢, and Figure 5 shows the dependence of the change in 


the radial component of normal stresses o,(r,z,f) along coordinate r at different points in time: 


1-t =tmax, 2—t =10t max (thax =1 s), solid line — z= 0, dotted line — z= A. 


W*(0,z,t)x107,m 


—10 


—15 
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Fig. 4. Diagrams W*(0, z, t) —¢ 


Gr(r, Z, 2)x10-7, Pa 


0.2 0.4 0.6 0.8 


N 


~ 


Fig. 5. Diagrams o,,(r, z, 2) — 1: 1 —t = tary 2 — t= 10tnax 


It should be noted that under heating, the plate bends with increasing thickness; due to fixation, compressive normal 
stresses o,,(r, z, f) are formed at all points. In the case of complete heating of the structure (¢= 10tnax), the value of normal 
stresses 67(r, z, 10tmax) in the height of the section practically coincide (Fig. 5, Diagram 2, solid and dotted lines). At this, 
or(r, 0, t) remains constant over the entire time interval t > tnax (Fig. 5, solid line), and on the lower plane at the initial 
moment of time o,(r, h, f), it is significantly less (Fig. 5, Diagram 1, dotted line). 

For a qualitative assessment of the induced electric pulse in the form of a potential difference V(f) (31), two electrodes 
with a radius of separation R = 0.7 and connected to a measuring device (Fig. 6, solid line) must be formed on the upper 
front surface of the element in question. At this, determination of V(¢) by connecting the upper and lower (grounded) solid 
electrodated surfaces of the plate to the voltmeter (Fig. 6, dotted line) is ineffective. 


V(t)x10° 


0.5 


0 2x 104 4x10 6x 104 t 


Fig. 6. Diagrams V(t) — ¢ 
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Discussion and Conclusion. The developed closed solution of the coupled axisymmetric thermoelectroelasticity 
problem for a round plate made of piezoceramic material is more accurate than the solution that was developed when 
solving problems in an uncoupled formulation. This is due to the fact that the calculated ratios obtained make it possible 
to determine how the non-stationary temperature field affects the stress-strain state and the electric field of the element 
in question, which makes it possible to describe the behavior of a round piezoceramic plate under the influence of 
thermal and electrical loads with greater accuracy. In addition, it becomes possible to scientifically establish the 
dimensions of two uncoupled electrodes, which provides measuring the induced electric pulse most effectively. 
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Abstract 

Introduction. All polymer materials and composites based on them are characterized by pronounced rheological 
properties, the prediction of which is one of the most critical tasks of polymer mechanics. Machine learning methods open 
up great opportunities in predicting the rheological parameters of polymers. Previously, studies were conducted on the 
construction of predictive models using artificial neural networks and the CatBoost algorithm. Along with these methods, 
due to the capability to process data with highly nonlinear dependences between features, machine learning methods such 
as the k-nearest neighbor method, and the support vector machine (SVM) method, are widely used in related areas. 
However, these methods have not been applied to the problem discussed in this article before. The objective of the 
research was to develop a predictive model for evaluating the rheological parameters of polymers using artificial 
intelligence methods by the example of polyvinyl chloride. 

Materials and Methods. This paper used k-nearest neighbor method and the support vector machine to determine the 
rheological parameters of polymers based on stress relaxation curves. The models were trained on synthetic data generated 
from theoretical relaxation curves constructed using the nonlinear Maxwell-Gurevich equation. The input parameters of 
the models were the amount of deformation at which the experiment was performed, the initial stress, the stress at the end 
of the relaxation process, the relaxation time, and the conditional end time of the process. The output parameters included 
velocity modulus and initial relaxation viscosity coefficient. The models were developed in the Jupyter Notebook 
environment in Python. 

Results. New predictive models were built to determine the rheological parameters of polymers based on artificial 
intelligence methods. The proposed models provided high quality prediction. The model quality metrics in the SVR 
algorithm were: MAE — 1.67 and 0.72; MSE — 5.75 and 1.21; RMSE — 1.67 and 1.1; MAPE — 8.92 and 7.3 for the 
parameters of the initial relaxation viscosity and velocity modulus, respectively, with the coefficient of determination 
R? — 0.98. The developed models showed an average absolute percentage error in the range of 5.9-8.9%. In addition to 
synthetic data, the developed models were also tested on real experimental data for polyvinyl chloride in the temperature 
range from 20° to 60°C. 

Discussion and Conclusion. The approbation of the developed models on real experimental curves showed a high quality 
of their approximation, comparable to other methods. Thus, the k-nearest neighbor algorithm and SVM can be used to 
predict the rheological parameters of polymers as an alternative to artificial neural networks and the CatBoost algorithm, 
requiring less effort to preset adjustment. At the same time, in this research, the SVM method turned out to be the most 
preferred method of machine learning, since it is more effective in processing a large number of features. 
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AHHOTalna 

Beedenue. J\na Bcex NMOMMepHbIX MaTepHasIOB HM KOMMO3HTOB Ha HX OCHOBe XapaKTepHbI ABHO BbIpaxKeHHBIe 
peoslorMyeckue CBOMCTBAa, IPOTHO3MpOBaHHe KOTOPHIX ABIIACTCA ONHOM 3 BaxKHeMIMX 3aad MeXaHHKM MOIMMepos. 
Bonpuive BO3MO2%KHOCTH JIA UPOrHO3HPpOBaHHA peOOrM4ecKHX TapaMeTPpOB TOJIMMepOB OTKPBIBAaIOT MeTOJBI 
MallMHHoro oOy4eHHa. PaHee poBOAMIMCch UcclIeqoBaHua Ha WpeqMeT MOCTpoOeHuA MpOrHO3HBIX MOJeleH c 
MCHOJIb30BaHHeM HCKYCCTBCHHBIX HelpOHHBIX ceTeli HM amroputTMa CatBoost. Hapayy c 9TuMu MeToZaMH, Onaroyaps 
BO3MOXKHOCTH OOpabaTHIBaTh JaHHbIe C CHJIbHO HeJIMHeHMHbIMM 3aBHCMMOCTAMH Me@2xKTY MpH3HakaMu, IWIMpoKoe 
IIPHMeHeHHe B CMEXKHBIX OOACTAX HAXOJAT MCTObI MANIMHHOTO oOOyyeHuaA — MeTor k-OnmKaliIHXx coceyel H MeTOT 
ONOPHBIX BeKTOpoB (SVM). OqHako paHee K IpoOsieMe, paccMOTpeHHOH B JaHHOl CTaTbe, ITH MCTOJBI He IPHMeHAIIMCh. 
Llemb1o paOoTbl aBuIacb pa3pa0oTKa IIporHo3HOM MOJeIH WIA OICHKH peoNOrM4yecKHX TWapaMeTpoB MosMMepoB 
MeTOJJaMH HCKYCCTBCHHOTO HHTeIIeKTa Ha IipHMepe MOJMBHHWIXJIOpHa. 

Mamepuanoi u memoodoi. B padote mpuMeHeHbI MeToy k-OnmKaMINIMX coceyqeH WH MeTO ONOPHbIX BeCKTOPOB JIA 
onpeyeseHua peoworMyecKHx MapaMeTpoB TOJMMepoB Ha OCHOBe KPHBBIX pellakcallwH HallpsKeHHM. OdOyaenne 
MOJesel BEIMOJHATIOCh Ha CHHTCTHUCCKUX JJAHHBIX, CrCHePHPOBAHHBIX Ha OCHOBe TeOpeTH4YeECKNX KPHBBIX pesakcallHn, 
TIOCTPOCHHBIX C HCIOIb30BaHHeM HesIMHeMHOrO ypaBHeHua Maxcpema-l ypepuua. BxoqHbIMM MapaMeTpaMu Moelle 
BBICTYMaIW BeIMYNHa TedopMalluu, Ip KOTOPOM NpOv3BOTMJICA IKCIepHMeHT, HavasIbHOe HallpwKeHHe, HallpsKeHNe 
B KOHIe Mpolecca pewlakcalluu, BPeMA pellakcallMM U yCOBHOe BpeMaA OKOHYAHHA Mporecca. BEIxoMHbIe WapaMerTps!: 
MOJyJIb CKOpocTH HM KoopduIMeHT HavabHOW peslakcaltwoHHOH Ba3KocTu. Moyen pa3padoTaHbBI B cpeye 
Jupyter Notebook Ha a35rKe Python. 

Pe3zyibmamoti ucciedoeanua. TloctpoeHbl HOBbIe IIpOrHO3Hble MOJeIM IIA OlIpeyeueHuA peouOrM4ecKuXx MapaMeTpoB 
TIOIMMepoB Ha OCHOBe MeTONOB MCKyccCTBeHHOrO MHTeeKtTa. IIpeqnoxKeHHbIe MOeuM OOecIedHBaIOT BbICOKOe 
KayecTBO MporHo3HpoBaHna. Metpuku kayecTBa Moyen B alroputrme SVR coctapiaiot: MAE — 1,67 u 0,72; 
MSE — 5,75 u 1,21; RMSE — 1,67 u 1,1; MAPE — 8,92 u 7,3 aya MapamMeTpoB HayasIbHOH peslakcalMOHHOH BASKOCTU 
HM MOJlyJIA CKOPOCTH COOTBETCTBEHHO C KOIPPMIMeEHTOM JeTepMuHauMU R? — 0,98. PaspaSoraHHble Moje NOKa3aIH 
CpeHIOIO aOCOuIOTHYIO MPOIeHTHyIO omMOKy B uama3oHe 5,9-8,9%. Tlomumo cuHHTeTHYeCKHX aHHBIX, 
pa3paOoTaHHBle MOJesIM TaloKe alpoOupoBaslach Ha peasIbHBIX IKCIePHMCHTAJIBHBIX JAHHBIX JIA WOTMBUHWIXJIOpua B 
qMata30He Temilepatyp oT 20 Wo 60 °C. 

OG6cystcoenue u 3aKirouenue. AtipoOata pa3paboTaHHbIxX MOjeei Ha peasIbHbIX 3KCIICPHMCHTaJIbHBIX KPHBbIX 
ToKa3ajIa BbICOKOe Ka4YeCTBO HX alMPOKCHMAalIIHH, COMOCTaBHMOe C UpyrHMu MeToaMu. TakuM o0pa30M, asITOpHTMBI 
k-Onmxaimux cocezeti u SVM MoryT Hcnonb30BaTbCA [IA IPOTHOSHPOBaHHA PeOMOTH4YeCKNX MapaMeTPoB MOJMMepoB 
Kak ajJIbTepHaTHBa HCKYCCTBCHHBIM HelpOHHbIM ceTaM HM alropuTMy CatBoost, TpeOyrollad MeHBINMXx ycuIMi 0 
IIpeBapuTesbHOH Hactpolike. II[pu 9TOM B WaHHOM UccieqOBaHuU HanOoliece IpeMOUTHTeIbHbIM MCTOJOM MallIMHHOro 
oOyyeHna OKa3asica MeTo SVM, Tak Kak OH Ooslee a(p(PeKTHBEH B OOpadoTKe OOOO YMCA MIPpH3HAKOB. 


KosroueBble CJIOBAa: peoIOorHaA, NOJIMMepsi, MCKYCCTBeHHBIM HMHTCIVICKT, MalLIMHHOe ooyueHue, k-OnuKaiimme cocean, 
OMOpHadA BCKTOPHadA perpeccna 


BaarogqapHocru. ABTODBI BbIpaxKaroT OaroapHOcTB petakWHu U peweH3eHTaM 3a BHHMAaTeCJIBHOeC OTHOIMMCHHE K CTaTbe 
MW yKa3aHHble 3aMCvaHHA, KOTOPble MO3BOJIMJIM HOBbICHTb Ce KayecTBO. 


Asia WaTupoBanusa. Konypatbesa T.H., Yerrypuenxo A.C. IIporHo3uposanue peouorMuueckux WapaMeTpoB MOIMMepoB 
MeTOqaMH MalMHHoro obyyeHua. Advanced Engineering Research (Rostov-on-Don). 2024;24(1):36-47. 
https://doi.org/10.23947/2687-1653-2024-24-1-36-47 


Introduction. Polymers are used in various industries, including the production of plastics, textiles, packaging 
materials, and more. Accurate prediction of the rheological parameters of polymers is a complex task that is important 
for optimizing production processes and creating products with desired properties. 

Today, machine learning methods have gained great popularity in various fields, including chemistry and materials 
science, due to their ability to efficiently process and analyze large amounts of data. These methods make it possible to 
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predict the properties of materials. In [1], a platform based on machine learning was described, and the integration of 
metrological support in the context of digital transformation was proposed. In [2], the local distribution of deformation, 
the development of plastic anisotropy, and fracture in additively manufactured alloys were predicted. The problems of 
developing measuring control regulators on digital platforms were formulated in [3]. An intelligent model for controlling 
the parameters of overlap joint welding was built in [4]. However, the issues of using machine learning methods to predict 
the rheological properties of polymers remain insufficiently investigated. This is caused by both technical and 
methodological difficulties, such as the heterogeneity of the polymer structure, their sensitivity to external conditions, 
and complex interactions between molecules during deformation. 

Research in the field of rheological properties of polymers and composites using machine learning methods has great 
prospects in the construction industry [5]. For numerous polymers, the experimental data are well described by the 
generalized nonlinear Maxwell-Gurevich equation [6], which has the form for a uniaxial stress state [7]: 


os" _ ft” 
Ot ia 
f =0-f£8, (1) 
11 rn" 
ye * > 
NN NNo m 


where ¢* — creep deformation, f~ — stress function, o — stress, E. — module of high elasticity, 19 — initial 
relaxation viscosity, m“ — velocity module. 

Various intelligent machine learning models can be used to determine the rheological parameters of polymers, such 
as the initial relaxation viscosity (hereinafter just “viscosity”) and the velocity module [8, 9]. For example, one such 
model is a neural network that can be trained on generated datasets to determine optimal polymer parameters [10]. 

Prediction based on synthesized data is a fairly common practice, including for nonlinear optimization 
methods [11, 12]. One of the ways to generate data is the use of Rosenbrock, Himmelblau, and Booth functions [13], 
which are effectively applied to test optimization methods such as gradient descent methods, genetic algorithms, and the 
Newton method. This approach was applied in [14], where a data set based on theoretical stress relaxation curves using 
the nonlinear Maxwell-Gurevich equation was generated to test the efficiency of various optimization methods. 

In [15], several machine learning approaches were given to predict the durability of a reinforced concrete beam, such 
as a neural network of back propagation, linear and ridge regression, a decision tree, and a random forest. The input 
parameters of the study were both various characteristics of the material and their properties, depending on the 
environment (temperature, humidity). Finally, according to the results of the study, the back propagation model 
determined a more accurate forecast (85%), the average values (MAE) and MAPE were 1.13% and 14.5%, respectively. 

Another approach to solving inverse problems of creep theory using the neural network method is based on training a 
model on large amounts of experimental data. In [16], a neural network model was developed, which was trained on data 
obtained as a result of long-term experiments on polymer materials, and successfully predicted the viscoelastic behavior 
of these materials. The data obtained from experiments on samples of various materials were used for the study. 

Unlike the above-mentioned papers, the presented research is intended to promote the development of more accurate 
and reliable methods for predicting polymer properties, such as the k-nearest neighbor method and the support vector 
machine, which is important for various industries and science. 

The research objective was to develop a predictive model based on artificial intelligence methods for analyzing the 
rheological properties of polymers. Previously, the authors had already used a machine learning algorithm based on 
gradient boosting CatBoost to process stress relaxation curves [17, 18]. CatBoost is one of the most powerful machine 
learning algorithms applicable to solving not only regression problems, but also classification and ranking problems [19]. 

The CatBoost method can be useful for solving some tasks, but it also has its limitations and disadvantages. In this 
regard, there is an interest in using other algorithms mentioned earlier [20] to solve the problem. 

Materials and Methods. The generated data array is partially presented in Table 1. This array was formed on the 
basis of theoretical stress relaxation curves described by the Maxwell-Gurevich equation, according to the technique 
presented in [14]. The variation ranges of the velocity modulus and the initial relaxation viscosity in the generated array 
correspond to the real ranges for polyvinyl chloride in the temperature interval from 20° to 60°C. The total number of 
numerical experiments (7) was 30,000. 
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Table 1 
Table of initial data for training the model 
. as ie Oe | Bae aie Relaxation | Conditional end | Velocity | Viscosity 
Deformation, beginning of end of the : . r 
No time time of the module no, 10°, 
% the process oo, | process ox, ; 
MPa MPa th, h process fos, h m, MPa MPa's 
1 1.000 10.000 0.909 0.277 1.484 6.000 3.000 
2 2.000 20.000 1.818 0.109 1.003 6.000 3.000 
3 3.000 30.000 2.727 0.046 0.820 6.000 3.000 
4 1.000 10.000 0.909 0.861 4.615 6.000 9.333 
5 2.000 20.000 1.818 0.339 3.122 6.000 9.333 
6 3.000 30.000 2.727 0.142 2.552 6.000 9.333 
7 1.000 10.000 0.909 1.445 7.747 6.000 15.667 
29997 3 45 37.5 0.285 2.476 15 53.666 
29998 1 15 12.5 1.003 4.255 15 60 
29999 2 30 25 0.558 3.371 15 60 
30000 3 45 37.5 0.319 2.769 15 60 


The data set consisted of five input variables and two output variables. Input variables (unit measure): 
deformation — ¢ (%); stress at the beginning of the process — oo (MPa); stress at the end of the process — o~ (MPa); 
relaxation time — ¢, (h); conditional end time of the process — fo; (h). Output variables (unit measure): velocity 


module — m* (MPa); initial relaxation viscosity — 1o (in Table 1 and further, simply “viscosity”) (10° MPa:s). 


Values Oo,, Ow,, tn, and fos are schematically shown on the typical stress relaxation curve (Fig. 1). 


Fig. 1. Typical stress relaxation curve 


The k-nearest neighbor (A-NN) algorithm is based on the similarity analysis of nearby objects. The K-NN method is in 
great demand for solving various types of machine learning tasks. 

Formula (2) represents the general form of the algorithm, where w(i, x) — weight function evaluating the importance 
of the i-th neighbor. 
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F (x)=argmax yey aE = y|w(i,x). (2) 


The maximum total weight can be achieved for several objects at the same time. The entropy of this process can be 
adjusted using nonlinear sequence w(i, x) = [i < k]q' (exponentially weighted k-nearest neighbor method) provided that 
0<q<05. 

Representing a fairly simple machine learning algorithm, k-NN is well applicable to solving classification and 
regression problems. The advantages of this method are ease of implementation, no need for pre-training of the model. It 
is used for all types of data, including categorical and numeric. Disadvantages: a tendency to over-training (provided that 
kis too small), poor performance with large amounts of data, it is not possible to take into account the relationship between 
the signs. 

The support vector algorithm — support vector regression (SVR) — solves the problem of minimizing the sum of the 
mean absolute error. SVR is more resistant to outliers, unlike the least squares method, due to the regularization 
coefficient (C) and the “epsilon-insensitive tube” (¢). In this case, ¢ determines the width of the tube in which errors are 
ignored. Stochastic gradient descent is used to find the minimum of the function. 

The support vector machine learning algorithm is function F(x) of approximation and regularization of empirical risk, 
which converts training and test samples into output data for each object of the corresponding sample. Formula (3) 
represents the general form of the algorithm, (4) is a linearly separable sample, (5) is a linearly inseparable sample, where 


C — regularization coefficient, Mw, wo) — scalar product of vectors (feature and support vector), wi— weight 
coefficients. 
I 1 2 ; 
F(x)=CY),_,(I-Mi (a0) +5 poo? > min F (x) 6) 
tae tc . 
gel > min F(x); ‘s 


M;(wwo)>Li={1:]}. 


1) 42 : ; : 
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S047). 

Function K (ea) is a function of a pair of objects (eae) , 1 representable as a scalar product in some space H, for 


which transformation y:X —> AH takes place. Function K:X¥xX >R — kernel if K (x,x") =(w(x),w(x’)),» 


provided that K is symmetric: K(x,x’) = K(x,x’) and nonnegative definite: [fee (x) g(x')dxdx', Vg :X > R. 


The regularization coefficient is determined by the sliding mode control method. 

Advantages of the SVM method are as follows: high accuracy in classification problems in nonlinear spaces; ability 
to work with a large number of features (including categorical and numerical), generalize data (which provides 
applying the model to new data), work with data that are not linearly separable due to the use of kernel functions. 

Disadvantages of the SVM method include inefficiency of working with large amounts of data; low interpretability 
of the model; the requirement to configure numerous parameters, such as the type of kernel (its parameters, 
regularization parameters), etc. 

In this research, algorithms are developed in the Jupyter Notebook intelligent computing environment using 
machine learning methods. 


As a learning algorithm, function F(x) is considered. It transforms training sample {xi} 5 eX” and test 
sample ee se eX’ into output data when training {yds eX” and testing {yi ‘ 7 e X' for each object of the 
corresponding selection. The training of the vector of parameters w; €W is embedded inside the algorithm. 

Under the conditions of the presented problem: { vite — actual values of viscosity nj (at the beginning of the 
i=l 


relaxation process) and velocity module m’; { Vi " _, — predicted values of viscosity ng (at the beginning of the 


relaxation process) and velocity module m”. 
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The selection of such a parameter as the number of neighbors affects the generalizing ability of the developed 
model and is important for its correct operation. The most suitable algorithm for calculating distance based on data is 
Distance, in which the weights of objects are inversely proportional to their distance. Accordingly, in the case of closer 
neighbors of the query object, they have more influence than their neighbors located at a greater distance from the 
object. 

The data set was divided into training and test samples in a ratio of 75/25. In turn, 20% of the training sample 
became validation. The sample size was: training — Xjrain = 20,400; test — xresr = 6,000; validation — xe.q = 3,600. For 
variables Virain, Viest; Veval, the data were distributed in a similar way. 

To build the é-nearest neighbor model, the following parameters were selected: number of neighbors, sheet size, 
interval, and weight function. The range and functionality of the values for the configurable parameters are shown in 
Table 2. 


Table 2 
Parameter table for kK-NN model 
No Parameter Value Functional 
1 Number of neighbors (4) 3,5, 7,9 Determines optimal number of neighbors for query 
> SiseeeizeiG) 15, 20, 30 Determines speed of querying and required memory for 
storing the tree 
3 Interval (p) 1 (11), 2 (12) Defines power parameter (Minkowski metric) 
4 Weight function (w) ‘uniform’, ‘distance! Predicting weights 


To build the SVR model, the following parameters were selected: kernel type, kernel order, regularization 
coefficient (quadratic regularizer), . The range and functional values for the adjustable parameters are presented in 
Table 3. 


Table 3 
Parameter table for SVR model 

No Parameter Value Functional 

1 Kernel type ‘linear’; 'poly'; 'rbf'; 'sigmoid' Defines type of hyperplane (linear/nonlinear) 

2 Kernel order 1,2, 3,4, 5,7 Defines degree of polynomial function of kernel 

: : Solves problems of vector multicollinearity 
3 Quadratic regularizer (C) ;3;4;5; 7; 10 snd madelietaininip 
4 ‘ 0.1; 0.2: 0.5: 1: 1.5: 2:3 Determines deviation of the object (proximity 
measure) 


Research Results 

Figure 2 shows the correlations between the variables. 

The following types of linear correlations between individual input and output variables of the model can be noted: 

— strong enough — between the variables “Deformation” and “Stress at the beginning” p.)< = 0.93; “Relaxation 
time” and “End time of the process” (;,1; = 0.93; 

— average — between the variables “Deformation” and “Stress at the end” pg... = 0.71; “Stress at the beginning” 
and “Stress at the end” poos., = 0.75; 

- weak — between the variables “End time of the process” and “Viscosity” Py»; =0.58; “Viscosity” and 
“Relaxation time” py; = 0.46. 

The presence of a moderate correlation between variables or its absence indicates only the absence of a linear 
relationship; therefore, it is possible to have a nonlinear relationship between variables. 
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Fig. 2. Correlation matrix 
Table 4 shows the statistical characteristics of the original data set. 
Table 4 
Statistical characteristics of the original data set 
Parameter € 60 Ox tn tos m No 
Unit measure % MPa MPa h h MPa 10° MPa:s 
count 30,000.00 | 30,000.00 | 30,000.00 | 30,000.00 | 30,000.00 | 30,000.00 30,000.00 
mean 2.00 25.00 15.78 0.75 4.4] 10.50 31.50 
std 0.82 10.77 9.10 0.94 4.40 2.87 18.19 
min 1.00 10.00 0.91 0.00 0.07 6.00 3.00 
max 3.00 45.00 37.50 10.04 38.02 15.00 60.00 


The best parameters for the 4-nearest neighbor model were determined as a result of 5-block cross-validation (Table 5). 


Best k-NN model parameters 


Table 5 


Parameter Number of neighbors (4) Sheet size (n) Interval (p) Weight function (w) 
No 3 15 2 ‘distance’ 
m 5 15 2 ‘distance’ 
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The best parameters of the SVR model for viscosity parameters 1% (at the beginning of the relaxation process) and 


velocity module m* were obtained empirically (Table 6). 


Table 6 
Best parameters of SVR model 
Parameter Kernel type Kernel order Quadratic regularizer € 
Nb 'rbf 2 5 0.3 
m 'rbf 3 6 0.3 


The ratio between the real and predicted values for the ANN model in terms of the parameters “Viscosity” and 
“Velocity modulus” is shown in Figures 3, 4. 
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Fig. 3. Diagrams of prediction errors of k-NN, “Viscosity” 
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Fig. 4. Diagrams of prediction errors of k-NN, “Velocity module” 
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The ratio between the real and predicted values for the SVR model according to the parameters “Viscosity” and 
“Velocity module” is shown in Figures 5, 6. 
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Fig. 6. Diagrams of prediction errors of k-NN, “Velocity module” 
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The metrics of the developed models of k-nearest neighbors and support vectors are presented in Tables 7 and 8, 


respectively. 
Table 7 
Metrics of the developed k-NN models 
Parameter RMSE MAPE (%) R? train R’ test 
1.00 0.98 
0.99 0.98 
Table 8 
Metrics of the developed SVR models 
Parameter MAE MSE RMSE MAPE (%) R* train R’ test 
N16 1.67 5.75 1.67 8.92 0.98 0.97 
m 0.72 1.21 1.1 7.3 0.89 0.87 


In addition to synthetic data, the developed models were also tested on real experimental data presented in [13]. 
Experimental relaxation curves of polyvinyl chloride were used for various temperatures in the range from 20° to 60°C. 
In Figure 7, the experimental stress values at different temperatures at different points in time are marked with felt-tip 
pens, and solid lines show stress relaxation curves based on values m* and nj predicted by the models. 
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Fig. 7. Results of testing the model on experimental data 


Discussion and Conclusion. Figure 5 shows that the quality of prediction based on experimental data is quite high, 
specifically, for temperatures of 30°C, 50°C and 60°C. For other temperatures, the prediction quality is somewhat lower, 
which is due to the quality of the experimental curves themselves. It was necessary to extend the experiment time and 
wait for the curves to reach the horizontal asymptote. 

In this research, the most preferred method is the support vector machine (SVM). This is due to the fact that SVM can 
process data with a large number of features, which is important for the analysis of rheological parameters of materials. 
In addition, SVM works with nonlinear dependences between features, it is applicable to solve the regression problem, 
which is required to determine the rheological parameters of materials. 

However, the CatBoost method can also be effective in this task, especially, if there are categorical features in the data. In 
addition, CatBoost can process missing data, which can be important for analyzing rheological parameters of materials. 

The k-nearest neighbor method is less preferable in this task due to its low efficiency in processing a large number of 
features, as well as the presence of problems with high data dimensionality. 

In the course of the investigation, it has been shown that the use of machine learning methods makes it possible to 
effectively analyze and process large amounts of data, including information about the characteristics of polymers and 
their rheological properties. The model developed on the basis of such an analysis maintains high accuracy in predicting 
the rheological parameters of polyvinyl chloride, which is confirmed by the results of cross-validation and comparison to 


experimental data. 
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One of the key advantages of this approach is the ability to automate the process of predicting the rheological 
parameters of polymers, which reduces the time and cost of research and development of new materials. In addition, the 
model can be easily adapted to analyze other types of polymers and predict their properties. 

As a result of this research, a predictive model has been developed to evaluate the rheological parameters of polyvinyl 
chloride using artificial intelligence methods based on data of its characteristics and rheological properties. The model 
demonstrates high prediction accuracy and can be used to optimize the production and development of new polymer- 
based materials. 
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Abstract 

Introduction. Improving the calculation methods of mechanical engineering facilities is an urgent and in-demand task. 
This fully applies to the techniques of calculating the strength of all-welded gastight boiler structures. Normative 
calculation techniques are based on simplified models that give limited possibilities for optimizing these structures. The 
low calculation accuracy inherent in such techniques is unacceptable under real design conditions, when an engineer is 
faced with the task of developing competitive structures in a short time, i.e., reducing metal consumption while providing 
the strength of these structures with limited development time. The use of simplified models was justified in the past, 
under conditions of insufficient development of computer technology. Application of the most advanced techniques based 
on computer modeling makes it possible to increase the accuracy of calculations, provide the optimization of such 
structures, and improve the quality of design. The objective of this study was to develop a new special procedure for 
calculating the strength of all-welded gastight structures based on computer modeling, using the most advanced methods 
of modeling the membrane wall and factors affecting it. The accompanying task was to verify the developed procedure 
based on comparing the results of calculations using the developed technique and the normative method. 

Materials and Methods. The developed technique is based on the replacement of the membrane wall with an orthotropic 
plate or shell. Computer modeling was used applying the finite element method of all-welded gastight structures, and the 
impacts to which they were subjected during operation, as well as an effective method for assessing the technical condition 
of these structures. 

Results. A new two-stage technique for calculating the strength of increased accuracy of all-welded gastight boiler 
structures was developed and patented. The calculation results were compared according to the proposed procedure and 
the normative method. It was shown that the proposed technique made it possible to increase the accuracy of modeling 
and calculation. The error in calculating all-welded gastight structures of a high-power boiler was reduced by more than 
30% for the recommended steps between stiffeners. For specially reinforced membrane walls with steps exceeding the 
permissible values, the error reduction reached 70% or higher. 

Discussion and Conclusion. The developed technique is used in the modeling and calculation of all-welded gastight 
structures. Its application makes it possible to optimize the step between the stiffeners of the structure of the support and 
connecting nodes of gastight membrane walls. Based on the results of the application of the two-stage calculation 
procedure, new designs were developed and patented. The developed technique has been used in the real design of boilers 
since 2014. 


Keywords: boiler unit, all-welded gastight structures, membrane walls, plates, shells, orthotropic plates, mathematical 
modeling, finite element method 
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AHHOTalna 

Beedenue. CopepiicHcTBoBaHHe pacuéTHbIX MeTOHK OObCKTOB MalIMHOCTpoeHHA — aKkTyaJIbHad HU BOCTpeOoBaHHaA 
3aqaua. B momHOH Mepe 3TO OTHOCHTCA HM K MeTOMKaM pacuéTa Ha IpOU4HOCTb IebHOCBapHbIX Ta30IJIOTHBIX 
KOHCTpyKUMi KoTIOarperaToB. HopMaTuBHbie pacuéTHbIe MCTOAMKH OCHOBAaHbI Ha YIPOWGHHBbIX MOJeIAX, HMCIOWIHX 
OrpaHHYeHHble BOSMOXKHOCTH JIA ONTHMM3aMM ITUX KOHCTpyKuMM. Hu3kadq TOUHOCTA pacuéta, IpHcyulad TaKHM 
MeTOJUKaM, HelIpHeMJIeMa B YCJIOBHAX PeasIbBHOTO IPOeKTHPOBAHHA, KOra Mepey MHKeHEPOM CTOAT 3aa4n paspaooTKH 
B CKATbIc CPOKH KOHKYPCHTOCHOCOOHBIX KOHCTpyKUMH, TO eCCTb CHWKeHHA MeTaJIIOEMKOCTH pH obecreyeHHu 
IIpOYHOCTH ITHX KOHCTpyKUMM HW OrTpaHH4eHHOM BpeMeHH pa3paboTKH. Ucnomb30BaHHe yiIpoueHHbIx MOJeNeH OBO 
ompaByaHo B TIpOWWIOM, B YCJIOBHAX HeJOCTaTOUHOTO pa3BHTHA KOMIIbIOTepHOM TexHuKH. IIpumMeHenHve HanOosee 
COBEPUICHHBIX Me€TOJMK, OCHOBaHHBIX Ha KOMIIbIOTEPHOM MOJCJIMPOBaHHH, MO3BOJIACT MOBbICHTb TOUHOCTb pacuéTos, 
oOecre4YHTb ONTHMH3alMI0 TaKHX KOHCTPyKUHi, YIy4UINTb KayecTBO mpoekTupoBanna. Lemp aHHoro 
YiccIeqOBaHHa — pa3paOoTKa HOBOM CielMasIbHOM MeTOAMKH pacuéTa Ha MPOYHOCTS IeIbHOCBapHBIX Ta30IJIOTHBIX 
KOHCTpyKUMi, OCHOBaHHOM Ha KOMIIbIOTepHOM MOJeIMpOBaHHH, C IpHMeHeHHeM HaMOosIee COBeEPUICHHBIX MCTOAHK 
MOJeIMpOoBaHHa MeMOpaHHoro 9KpaHa MH *akTopos, BoO3sAelCTBYIOWMX Ha Hero. ConmyTcTByIolet 3afa4ueli aBTOpoB 
CTaTbH ABJIANIACh BEPHPUuKaluA paspadOoTaHHOM MeTOAMKH Ha OCHOBEe CpaBHeHHA Pe3yJIbTATOB pacueTOB C IIpHMeHeHHeM 
pa3pa0oTaHHO MeTOAMKH H HOPMaTHBHOrO MeTOLA. 

Mamepuanoi u memoodoi. PaspadoTaHHad MeTOAHKa OCHOBaHa Ha 3aMeHe MeMOpaHHoro 9KpaHa OpTOTpOMHoH 
TWIacTHHOH vWIM OOon0UKO. Ucromb30BaHbl KOMIbIOTepHOe MOseMpoBaHHe C MIpHMeHeHHeM MeTOJa KOHCYHBIX 
QJIEMCHTOB ICJIbHOCBapHBIX Ta30IJIOTHBIX KOHCTPpyKUHMM MW BO3eHCTBHM, KOTOPbIM OHH TOABepxKeHbI B Mpotecce 
SKCIIyaTayun, a TakxKe 3*(eKTHBHBIM MeTOA OLWCHKH TEXHHYCCKOTO COCTOAHHA STHX KOHCTpyKUHH. 

Pe3yismamoet uccnedoeanua. Pa3spadoTaHa HOBad JBYX3TalHad MeTOAMKa pacuéTa Ha MpOYHOCTb LWesIbHOCBapHBbIX 
Ta30MWIOTHBIX KOHCTpyKUMM KOTIOArperarTos, MOJYYMBIat WaTeHT Ha u300petenne. IIpopeyxeno cpaBHeHHe pe3ysIbTaTOB 
pacyéTosB mo mpeszoxeHHOM MeTOAMKe HM TO HOpMaTHBHOMy MeTogy. Iloka3aHo, 4TO MpesOwKeHHat MeTOMKa 
MO3BOJIA€T MOBbICHTb TOUHOCTh MOJeMpoBanua u pacuéta. IlorpemiHocTb pacuéta WebHOCBapHBbIxX Ta30IJIOTHBIX 
KOHCTpyKUHMi KOTIAa OObUIOM MOUIHOCTH CHYKeHa Oosee 4em Ha 30 % AIA peKOMCHAOBaHHBIX WaroB MexKTYy WoAcaMu 
xecTKocTu. Jia NoAKpenéHHbIx cCielMabHbIM OOpa30M MeMOpaHHBIX 9KpaHOB Cc WaraMH, MpeBbIWalolujuMu 
JOMYCTHMBIe 3HaveHHA, CHYDKeHHe MorpeliHoctTu ZocTuraet 70 % u BBILIUe. 

O6cyotcoenue u 3akio4uenue. PaspadoTaHHad M€TOAMKa MCIOJb3yeTCA IPH MOJeIMpOBaHHH UM pacuéTe WebHOCBapHBIX 
Ta30MJIOTHBIX KOHCTpyKuMH. IIpumMenenve ee MO3BOIAeT ONTHMVM3HpOBaTb War MexKY MOACAaMH 2KeCTKOCTH 
KOHCTPyKUMH OMOPHBIX H COCMHUTEIbHBIX Y3IOB Ta30IWIOTHBIX 3KpaHos. Ilo pe3sybTaTaM MIpHMeHeHHA AByXITANHOH 
pacuéTHOH MeTOAUKU OBLIM pa3spaOoTaHbI HOBbIC KOHCTPYKIMH, MOYYMBLINe MaTeHTbI Ha H300peTenuaA. PaspadboTaHHaa 
MeTOJHKa IIpHMeHACTCA B peasIbBHOM IIPOCKTHPOBAHHH KOTIOarperatTos c 2014 roza. 


KosrroueBble cJI0Ba: KOTJOarperaT, UWCJIbHOCBapHbIe Ta30WJIOTHble KOHCTpyKUHH, MeMOpaHHEIe SKpaHbl, WIacTHHbl, 
o0o00UKH, OpTOTPOTIHbIle WlacTHHbl, MATCMaATHYUCCKOC MOJCIINPOBaHHe, MCTO], KOHCUHBIX SJICEMCHTOB 


BaarogapHocrn. ABTOpBI BbIpaxKaroT OaroqapHOcT petakWHv U peleH3eHTaM 3a BHHMAaTeCJIbBHOe OTHOMICHHE K CTaTbe 
MW yYKa3aHHble 3aMC4aHHA, KOTOPble NO3BOJIMJIM NHOBbICHTb Ce KadecTBO. 


Aaa wnTupopanua. Kypenun M.II., CepOunoscxuit M.IO. I[pumMenenne crenMambHbIxX pacuéTHBIX MeTOJMK IIpH 
IIpOCKTHPOBaHHM IeIbHOCBaPHbIX [a30NJIOTHbIX KOHCTpyKUM KoTIOarperatoB. Advanced Engineering 
Research (Rostov-on-Don). 2024;24(1):48-—57. https://doi.org/10.23947/2687-1653-2024-24- 1-48-57 
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Introduction. The normative technique of the strength calculation of all-welded gastight structures! is based on 
replacing the screen with a structurally orthotropic Kirchhoff-Love plate with considerable simplifications. In this 
model, the orthotropy of the plate properties is taken into account only when calculating bending forces. The major 
objective of the introduced simplifications is to reduce the complexity of calculation formulas and, ultimately, the 
overall complexity of calculations. The models underlying the technique include a gastight screen in the form of an 
orthotropic plate, stiffener rings in the form of constant-section beams, and fasteners in the form of connections 
between the wall and the stiffener rings. 

The following assumptions are made in the model’: 

1. The orthotropy of the gastight wall is taken into account by introducing ratio D, / D, when calculating the bending 
moment WM, for an isotropic wall with cylindrical stiffness in a section perpendicular to pipes D,, and in a section parallel 
to pipes D.. 

2. When determining stresses at the midpoint of the panel, there is a hinged connection of the gastight screens to each 
other in the screens of the all-welded box, since the fastening condition affects the nearest 6-8 pipes from the fixing point. 

3. When determining the stresses in the corner of an all-welded box, there is pinching. 

4. The action of two groups of forces is considered: active (boost and “buckling’”’) and reactive (emphasis on 
stiffener rings). The remaining impacts (from the internal pressure in the pipes, weight and temperature factors) are 
determined independently and taken into account at the stage of assessing the technical condition by the superposition 
of forces method. 

5. The calculation is based on the principle of assessing the structural strength by bearing capacity, which is 
determined by the limiting state of the transition of the most loaded section from an elastic state to a plastic one. The 
calculation is performed according to a conditionally elastic scheme. 

A significant disadvantage of this scheme is the use of simplified models that do not fully take into account the 
reinforcing, attached and other structural elements, their cross impact, as well as the effect of all influencing factors 
causing a complex stress-strain state of the membrane wall. This technique has limited modeling capabilities, specifically 
for membrane walls of complex shape, for which it is not applicable. 

There are limits to the application of normative calculation formulas, constraining them due to design or loading 
features*. For membrane walls, standard calculation formulas are not applicable for any cross-sections, except for a pipe 
with a flat spacer. They are not applicable for membrane walls of complex shape. The maximum step size between 
stiffener rings is limited. It is noted in the normative technique that it is allowed to use other methods of calculating the 
strength of all-welded membrane structures, subject to providing the normative strength reserves. The technical 
regulations of the Customs Union enable using a number of calculation methods in addition to regulatory formulas, 
including computer modeling‘. In regulatory documents, the strength calculation techniques that go beyond the 
application of regulatory formulas are called alternative, or special. In connection with the above, the objective set by the 
authors of this article — to develop a new special method for calculating all-welded gastight structures using computer 
modeling, which allows for high accuracy and reliability of modeling and calculation results — seems highly relevant. 

In the papers of Russian authors, models have already been proposed in which the membrane wall is replaced by a 
statically indeterminate frame. Special calculation methods are successfully used to assess the strength and resource of 
elements of boilers and power plants [1]. In the studies devoted to the calculation of gastight membrane walls, it is noted 
that the methods based on the use of modern achievements in the area of numerical methods using computer modeling 
are promising [2]. 

Recent studies on boiler membrane walls include the works of MiloSevic-Mitié [3], Nagiar [4], Serti¢ [5], who 
proposed a method for calculating boiler strength and determining temperature displacements and loads of supports based 
on the Kirchhoff-Love theory and computer modeling using the finite element method (FEM) [6]. However, this technique 
is characterized by an increased error due to the application of significant simplifications in the membrane wall model 
(e.g., pipes are replaced with absolutely rigid bodies) [7]. The method of reverse transition from the plate to the wall and 
assessment of its technical condition has not been proposed [8]. 


‘RD 10-249-98. Calculation Standards for the Strength of Stationary Boilers and Steam and Hot Water Pipelines. URL: 
https://docs.cntd.ru/document/1200021653 (accessed: 15.12.2023) (In Russ.). 

? RTM 24.031.06-73. Calculation and Design of Fastening Elements of the Wall System of Boiler Units with All-welded Panels: Technical Guides. Leningrad: RIO 
NPO TsKTI; 1974. 39 p. (In Russ.). 

3 GOST 34233.1-2017. Vessels and Devices. Norms and Methods of Strength Calculations. Moscow: Standartinform; 2019. 30 p. (In Russ.). 

*CUTR 032/2013. Technical Regulations of the Customs Union. On the Safety of Equipment Operating under Excessive Pressure. URL: 
http://docs.cntd.ru/document/49903 1170 (accessed: 15.12.2023) (In Russ.). 
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High calculation accuracy was obtained through modeling membrane walls based on the general mechanics of a 
deformed solid and measurements in numerical experiments with accurate solid-state models of the screen [9]. Effective 
methods for modeling loads [10] acting on all-welded gastight structures have been proposed, which make it possible to 
increase the accuracy and reliability of modeling these structures [11]. 

In general, modeling membrane walls using orthotropic plates is relevant if we take into account the complexity of 
modeling and calculating membrane walls with many pipes and spacers using traditional solid-state computer simulation 
techniques, for which there are strict limitations. These restrictions deal with the relative dimensions of the height of the 
solid-state element to its length and width, which cannot exceed ' (aspect ratio metric [12]), with a ratio of membrane 
wall thickness to its width of 1/400 or more. 

When developing the new technique, the following objectives were set: increasing accuracy, reducing the complexity 
of modeling and calculation, reducing the time of computational operations during optimization and evaluation of 
strength, rigidity, stability, and durability of structures. It is known that membrane walls operate under conditions of 
complex resistance and are simultaneously exposed to external forces, moments, pressure or rarefaction from the inside 
of the boiler, pressure in pipes, an uneven temperature field, taking into account the cyclic effect of these parameters and 
creep of materials at high temperatures. An additional, but rather important task of this development was to create a 
convenient tool for cyclic optimization of the structure, which includes visualization of the results of making changes to 
the structure in the form of stress, strain, and displacement fields. 

Materials and Methods. The technique proposed by the authors involves replacing the membrane wall with an 
orthotropic plate or shell. The membrane wall has a regular structure and a section that periodically changes in the 
direction perpendicular to the pipes. Such a membrane wall can be considered as a structurally orthotropic plate, i.e., a 
plate in a flat stressed state with elastic characteristics that differ in mutually perpendicular directions. The membrane 
wall is a plate, and the bent screen is a shell since the diameter of the pipes, which determines the thickness of the screen, 
is small compared to its other dimensions. The structural elements of a gastight boiler, which can be represented as 
structurally orthotropic plates, are, first of all, membrane walls of the furnace, the transition flue, the convective shaft of 
the boiler, and membrane superheaters. 

The authors previously developed a two-stage modeling and calculation technique for all-welded gastight structures, 
which received a patent for the invention [13]. It is based on the replacement of the membrane wall with an orthotropic 
plate or shell, reliable modeling of all-welded gastight structures and the impacts to which they are exposed during 
operation, and on an effective method for assessing the technical condition of these structures. Through the developed 
technique, there are opportunities to optimize structures and reduce their metal consumption [14], which have not existed 
before. New designs are worked out, the complexity of design is significantly reduced [15]. The problem of assessing the 
technical condition of boiler structures containing membrane walls is solved by applying the most advanced modern 
practices based on achievements in numerical modeling. 

The method of calculating the strength of all-welded gastight structures consists of two stages. At the first stage, a 
model of a box structure with screens in the form of orthotropic plates and/or shells is formed with a preliminary 
calculation of their dimensional and physico-mechanical characteristics. The attached elements, including pipe connectors 
with headers, festoons, stiffener rings, various fasteners, supports, reinforcement elements are modeled in the form of 
rods, shells, connections, and boundary conditions. Then, local zones of orthotropic plates and/or shells with increased 
displacements and deformations are determined, and the technical condition of the attached elements is assessed. At the 
second stage, models are formed for places with exceeding the specified parameters, in which local zones of screens with 
attached elements with increased displacements, stresses, and deformations are modeled by solid-state elements 
(submodeling of local zones [12]). Conditionally elastic analysis is carried out. Health assessment in critical sections is 
provided according to static strength conditions and optimization of the structure. 

Research Results. The calculation of all-welded gastight boiler structures with a steam capacity of 810 tons of steam 
per hour (TPE-360/T model) with dimensions of 14.48x14.24 m was performed. In the course of the work, the authors 
compared the calculation results using the proposed procedure and the normative method. 

Initial data for the calculation: the boiler membrane walls were welded from a 606 mm pipe of steel 20 with flat 
spacers 80 mm wide, made of sheets of steel 20 with a thickness of 6 mm. The membrane wall was exposed to heating 
from a furnace with a cooling medium in pipes with the following parameters: 16.3 MPa, 349°C (Fig. 1). The calculation 
was carried out for buckling 3,000 Pa and emergency rarefaction 5,000 Pa. 
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Fig. 1. Initial data for the calculation: heat flow perceived by the screen surface 


The calculation was carried out in accordance with the requirements and formulas of normative documents 
(“Calculation standards for the strength of stationary boilers and steam and hot water pipelines”, “Calculation and 
design of fastening elements of the screen system of boiler units with all-welded panels”, “Technical Regulations of 
the Customs Union”). 

The permissible step between the stiffener rings, according to the calculation results, is 4,012 mm at the furnace level 
and 4,050 mm in the upper part of the furnace at marks 16.8 and 50.6, respectively. 

The calculation of the boiler box using the technique presented in this paper was performed in two stages. The model 
of the first stage included gastight screens of a furnace, a convective shaft, a gas reversing chamber, an aerodynamic nose, 
and a ceiling superheater, made in the form of shells with elasticity parameters calculated using numerical methods [9]. 
The stiffener rings, festoons No. | and No. 2 with suspensions, suspensions of boiler screens, wall superheater were 
modeled with beam-rod elements (Fig. 2). When comparing the results of the two calculations, it was reported that the 
problem areas of the structures, determined by the first stage of modeling, converged with the design points of the standard 
calculation: in the middle of the screen under the stiffener rings (point No. 1, Fig. 2 C), in the corner of the furnace 
between stiffener rings (point no. 2, Fig. 2 B). At the same time, the new technique based on a more advanced model of 
the boiler box made it possible to take into account the following technical solutions applied in the design of the boiler 
box and aimed at optimizing and strengthening it: design of corner bracing in the form of lever mechanisms that allow 
the force to be transmitted strictly along the axis of the adjacent screen; intermediate fasteners reinforced with combs- 
ribs, providing free thermal expansion of the screen; longitudinal ribs under the stiffener rings in places where the standard 
values of the steps between the stiffener rings are exceeded [15]; special designs of stiffener rings [16]; angle brackets; 
seal boxes]. In addition, the model made it possible to take into account the effect of a wall-mounted radiant superheater 
(attached to the screens through welding), festoons, stiffen boxes, and other structural elements. 
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Fig. 2. Calculation of the boiler box of TPE-360/T model according to the developed technique: 
A — shell-rod model (first stage of modeling); B — solid molel of the box corner (second stage of modeling); C — solid molel of 
the midpart of the screen under the stiffener (second stage of modeling) 


The model includes structures of a ceiling superheater, an aerodynamic nose, a convective shaft, a gas reversing 
chamber, and other elements, whose technical condition is difficult to assess using the normative technique due to its 
limitations. The permissible step between the stiffener rings, according to the calculation, is 4,500 mm at the furnace 
level, and 6,000 mm at the top of the furnace, which is higher than the values obtained through the normative technique. 

Thus, it can be concluded that the impact of design solutions aimed at optimizing and strengthening the boiler box 
(without increasing the number of stiffener rings) is significant. The results of the calculation according to the normative 
technique, which does not allow taking into account these decisions, are excessively conservative. The calculation results 
using the proposed method are more reliable, as they enable us to fully take these decisions into account. 

A comparison of the calculation results obtained through different methods was carried out at two characteristic points 
of the structure specified in the normative technique: point 1, for which the influence of attached and reinforcing elements 
is minimal, at the calculated point 2, for which the influence of reinforcing elements is most significant (Table 1). 
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Table | 
Calculated stresses in screen sections according to the normative and new technique, MPa 
Design section, mark, m 16.8 50.6 
Permissible stresses in the cross section of 1—2 screen, MPa 149.8 152.6 
Permissible stresses in the cross section of 5—6 screen, MPa 158.4 161.8 
Calculated combination of loads: weight, overpressure in pipes, buckling in the furnace space 
Setting-out point according to the requirements of the normative * 7H re aH 
technique 
Equivalent stresses in the cross section of 1-2 soreen according 1193 121.7 1414 145.9 
to the requirements of the normative technique, MPa 
Equivalent stresses in the cross section of 5-6 soreen according 10.4 970 21.9 175.9 
to the requirements of the normative technique, MPa 
Equivalent st in th tion of 1-2 i 
quivalent stresses in the cross section 0 screen according 135.6 109.2 149.0 1144 
to the requirements of the developed technique, MPa 
Equivalent st in th tion of 5— i 
quivalent stresses in the cross section 0 5-6 screen according 18.1 55.4 774 92.8 
to the requirements of the developed technique, MPa 


Calculated combination of loads: weight, overpressure in pipes, emergency rarefaction in the furnace space 


Permissible stresses in the cross section of 1—2 screen, MPa 199.8 203.5 


Permissible stresses in the cross section of 5-6 screen, MPa 217.8 222.5 


Equivalent stresses in the cross section of 1—2 screen according to 


138.6 175.8 183.4 


the normative technique, MPa 


Equivalent stresses in the cross section of 5—6 screen according to 
the normative technique, MPa 


15.6 289.3 


Equivalent stresses in the cross section of 1—2 screen according to 
the developed technique, MPa 

Equivalent stresses in the cross section of 5—6 screen according to 
the developed technique, MPa 


109.0 


125.8 


* in the middle of the wall, between the stiffener rings 


** in the corner of the box 


At point 1, a good coincidence of the results was obtained — with respect to the permissible stresses, the voltage 
difference was no more than 11%. At point 2, the stress difference ranged from 33% for small and medium steps between 
the stiffener rings to 73% for large steps between the stiffener rings. High stress values in all cases were given by the 
normative technique. 

The proposed method and model increase the reliability of simulating the behavior of the structure and the accuracy 
of the calculation. One of the results of the calculations and optimization of the design is the distribution of stiffener rings 
according to the height of the boiler furnace of TPE-360/T model. Table 2 shows the results of calculating the step of 
stiffener rings: 1 — according to the normative technique, without taking into account the design features; 2 — according 
to the calculation results using the developed techniques that provide more reliable simulation of the behavior of the 
structure, taking into account the influence of all its elements, including additional reinforcing elements in the form of 
stiffener rings, etc. 

Table 2 
Calculated and accepted steps between the stiffener rings 


Design section!, m 8.3 | 16.8 | 21.3 | 25.8 | 29.1 | 32 | 36.1] 42 | 50.6 


Maximum permissible step between stiffener rings”, mm 4,012|4,012)3,961|3,974|/3,979|3,971|3,962/4,000/4,050 


The adopted step between the stiffener rings following the 
results of design optimization, mm 

Fulfillment of the strength conditions for the accepted step 
(yes/no) 


2,660/4,500}4,400/4,250)2,700)3,520)3,300)/3,500/6,000 


yes | yes | yes | yes | yes | yes | yes | yes | yes 


1 — mark for the height of the furnace, where 0 — the floor level 


2 — as per calculation according to the normative technique 


Kurepin MP, et al. Application of Special Calculation Techniques in the Design of All-Welded Gastight Structures of Boiler Units 


It can be observed that the use of special calculation methods makes it possible to avoid the installation of additional 
stiffener rings for at least two characteristic places of the boiler furnace: at the place of installation of the burners (marks 
from 16.8 to 25.8 m) and at the place of the reversing flue (50.6 m mark). This avoids overspending of metal and 
complication of the structure due to the difficulties of simultaneous placement of burner devices and stiffener rings, and 
elements of the reversing flue. 

An assessment of the metal consumption was carried out. In this case, the weight of one stiffener ring with fasteners 
is on average 3,500 kg. To meet the requirements of the regulatory calculation for the step between the stiffener rings, it 
is required to install at least three additional stiffener rings at the furnace height from 16.8 to 25.8 m and 50.6 m. The total 
metal consumption of such a constructive solution is 10,500 kg. Thus, the use of new special calculation techniques 
provides ample opportunities for optimization and reduction of metal consumption while maintaining sufficient strength 
of all-welded gastight boiler structures. 

Discussion and Conclusion. A new technique for calculating and analyzing the technical condition and optimization 
of membrane wall structures of boiler units has been developed on the basis of modern, most effective methods of 
mathematical modeling of membrane walls with equivalent orthotropic plates, calculation of their geometric and physico- 
mechanical parameters, special loads and impacts inherent in all-welded gastight structures. The use of new special 
calculation techniques provides ample opportunities for optimization and reduction of metal consumption while 
maintaining sufficient strength of all-welded gastight boiler structures. 

The results of verification of the replacement of membrane walls with orthotropic plates were compared to the results 
of field experiments on determining the angles of rotation of the vertices of the ribs under the action of a unit moment. 
They showed that the displacement error did not exceed 14% for shell models implemented through the proposed 
technique. At the same time, solid-state models give an error of up to 10%, i.e., slightly less when the number of finite 
elements of the models is an order of magnitude larger (and the dimension of the model is up to 100 times higher, 
depending on the type of element — linear or quadratic). This modeling error fully meets the requirements of modern 
design, since the errors associated with manufacturing tolerances are in the range of 10-15% of the wall thicknesses of 
pipes, sheets, and the weld leg [9]. 

The significant disadvantages of the normative technique are the scant possibilities for optimizing the design and 
strength requirements in other ways, except for increasing the number of stiffener rings, which significantly increases the 
metal consumption and causes growth of cost. This means that simplified models adopted in the normative calculation at 
the time when computer engineering was insufficiently developed, limit the possibilities of optimizing structures in 
modern design. This necessitates using special calculation techniques based on numerical computer modeling. 

The technique proposed by the authors during the actual design of boiler units has shown its high accuracy and 
efficiency, wide possibilities for optimizing the design, reducing metal consumption under maintaining sufficient strength, 
the ability to significantly reduce the complexity of model formation and machine calculation time. In addition, the 
reliability of behavioral modeling and calculation accuracy of all-welded gastight structures have been increased through 
the use of modern numerical simulation techniques. They provide most accurate accounting of the features of structures 
and their cross impact, and the effect of all influencing factors, as well as calculating screens of complex shape. Through 
the developed techniques, the calculation error has been reduced by more than 30% for the recommended steps between 
the stiffener rings. For specially reinforced membrane walls with steps exceeding the permissible values, the error 
reduction reaches 70% or higher. 

The application of the proposed technique makes it possible to decrease the number of finite elements of models of 
structures with screens by several orders of magnitude compared to similar solid-state models of the screen, reduce the 
complexity of forming models of structures with screens. The time of computational operations in calculating the stress- 
strain state of all-welded gastight structures is significantly improved. And finally, the complexity of modeling and 
calculation is significantly reduced while maintaining sufficient accuracy [13]. 

The authors note that the technique they developed for calculating the strength of all-welded gastight structures is 
generally intended for the use in the design and optimization of thin-walled structures with periodically changing cross- 
section and internal channels in the walls, which can be simultaneously exposed to external forces, moments. In this 
particular case, these are thin-walled structures made of tubular membrane walls of steam and hot water boilers, i.e., all- 
welded gastight structures with all structural elements attached to them, e.g., elements of the support and suspension 
system, fasteners, burners, blast nozzles, etc. 

This technique is implemented on medium-powered personal computers and can be put into practice of real design of 
boiler-building enterprises and design organizations of other branches of mechanical engineering, which use plates and/or 
shells with regularly changing cross-sections. 
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The proposed technique has been used in the real design of boiler units since 2014, including power boilers 
Pp-1030-25.0-570/570GM, E-220-9.8-540GM, E-500-13.8-560G (model TGE-440), PP-1900-25.8-568/568 KT, 
E-540-13.8-560GN (model TGE-225, E-540-13.8-560GN (model TGE-224/S7), E-810-13.8-560BT (model TPE-360/T). 
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Abstract 


Introduction. The combined effect on the phase transformation process, involving a combination of heat treatment and 
external action, is a major technology solution for obtaining the required properties of steel products. When hardening 
steel in a constant magnetic field with a strength of 1-2 MA/m, martensite formation is observed at higher temperatures. 
In addition, when compared to conventional hardening, there are changes in structure and properties. Such effects cannot 
be explained only in terms of thermodynamics, since the expected shift in the equilibrium temperature between austenite 
and martensite in a magnetic field of such strength does not exceed 4—-8°C. To explain the effects that occur during 
hardening in a magnetic field, it is proposed to consider the features of martensitic transformation in highspeed steel when 
exposed to an external magnetic field in the temperature range of austenite superplasticity. This research was aimed at 
identifying the features of martensitic transformation in the presence of a constant magnetic field in steel with account 
for the phenomena occurring in the premartensitic state. 

Materials and Methods. Samples made of steel R6MS5 were used. Characteristics of the martensitic transformation were 
studied using the potentiometric method of electrical resistance. The data were recorded using an L-CARD E14-440 
analog-to-digital converter with the LGraph2 software package. The sample was heated by passing current. The sample 
was placed in the interpolar space of an open-type laboratory electromagnet FL-1, which provided the creation of a 
magnetic field with a strength of 1.2 MA/m. 

Results. The obtained differentiated dependences were characterized by electrical resistance anomalies (low-temperature peaks) 
at a temperature corresponding to the appearance of a ferromagnetic phase as a result of martensitic transformation. In a 
magnetic field, the development of martensitic transformation began at a higher temperature, which could not be explained in 
terms of thermodynamics. Thus, the formation of stress martensite was observed in microvolumes of austenite with 
ferromagnetic ordering, which perceived the energy of the external field through magnetostrictive stresses. Under conditions of 
superplastic austenite, such stresses were sufficient to initiate shear transformation. The minimum possible size of lattice 
instability fluctuations (1.372 nm) was determined. 

Discussion and Conclusion. Exposure to a magnetic field under hardening intensified the processes of some magnetic 
decomposition of austenite. At temperatures close to the beginning of the martensitic transformation, the existing areas 
of magnetic inhomogeneity were superimposed on the effects of the phenomenon of instability of the crystal lattice. In 
the temperature range Mz-M,, when austenite exhibited superplasticity, the formation of stress martensite in microvolumes 
of austenite with ferromagnetic ordering was significantly facilitated. 


Keywords: stress martensite, magnetic field, superplasticity, lattice instability, steel, hardening 
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AHHOTalna 

Beedenue. KomOunuposanhoe BuIMAHHe Ha Wporecc asoBoro MpeBpall[eHuA, mpesnMosaraoujee CoueTaHve TepMM4eCcKOoH 
oOpaOoTKH C BHELIHHM BO3]{eHCTBHEM, ABJIACTCA AKTYAJIBHBIM TEXHOJIOTMYCCKHM PeIeCHHeM JVIA MOJYA4eHHA HEOOXOAMMBIX 
CBOMCTB cTaIbHOM mpoyyKuMu. IIpu 3akanKe cTasIw B MOCTOAHHOM MarHHTHOM Tose HalpsKeHHocTEIO 1-2 MA/M 
HaOjoJaeTca OOpa30BaHve MapTeHcuta pH Oosee BLICOKHX TemMIepaTypax. IlomuMo 9TOrO, MO cpaBHeHHIO c OOBMHOM 
3aKaJIKOM, MIPOHCXOJAT H3MeCHCHHA B CTpyKType u cBolicTBax. HloyqoOubIe addeKTEI He MOryT OOBACHATBCA TOJIBKO C 
TePMOJMHAMMYeCKHX MO3HUMH, Tak KaK TpequosaraeMblii CABHr TeMIepaTypbl paBHOBeCHA Mey ayCTeHHTOM HU 
MapTeHCHTOM B MarHHTHOM ToJIe TaKOM HallpmKeHHOCTH He MpeBpuuaetT 4-8 °C. JIna oObacHenua acpdexros, 


BO3HHKaIOWHX IIpH 3aKasike B MaTHUTHOM Moule, lpeylaraeTcA paCCMOTpeTb OCOOCHHOCTH MapTeHCHTHOLO ipeBpallleHHsA B 
ObICTpOpexKyLlei CTaIIM IIpH BO3AeHCTBHM BHELIHHM MarHHTHBIM IOJIeM B TeMIIepaTypHOM HHTepBasie CBepXxIIaCTHYHOCTH 
ayctennta. Llenbro WaHHoi paOoTbI cTaIO BbIABIIeHHe OCOOeHHOCTeli MapTeHCHTHOTO IIpeBpallleHHa B TIPHCyTCTBHM 
MOCTOAHHOTO MarHUTHOTLO MOJIA B CTA C YYCTOM ABJICHHM, BOSHHKaIOWMX B IIpeMapTeHCHTHOM COCTOAHHMH. 
Mamepuaaot u memoooi. VUctionb30Banucb oOpa3tibl cTasm Mapxu P6M5. UccneqoBanue ocobeHHocTeli MapTeHCHTHOTO 
TIpeBpalleHuA OCYLICCTBIAIH MOTCHIMOMeTpHYecKHM MeTOAOM 3JIeKTpoconpoTuBseHua. JlaHHble (PHKCHpOBaJIHCh C 
TOMOIMIbIO aHasloro-l[u@poporo mpeodpa3opatena L-CARD E14-440 c  ucnonb30BaHvem mporpaMMHoro 
Komiiexca LGraph2. Harpes oOpa3yja MpoBoAMsICA MpOXOWALJHM TOKOM. OOpa3ell pa3Melaicd B Me€KITOJIIOCHOM 
TIpoctpaHctTBe saOopaTopHoro 3IeKTpOMarHuta OTKpbITorO Tua DJI-1, KoToppii oOecne4nMBal Co3sqaHHe MarHuTHOrO 
TOJIA Hatips»KeHHOCTHIO 1,2 MA/M. 

Pe3yiemamoet ucciedosanua. Ha nomyyenubix uddepeHuMpoBaHHbIx 3aBHCHMOCTAX IIPHCyTCTBOBaIM aHOMaJIMH 
3JIEKTPOCOMpOTHBIIeHHA (HH3KOTeMIepaTypHble MHKH) pH TemMMepatType, COOTBeTCTBYIOWIeH  MOABJICHHIO 
(beppoMarHUTHOH a3bl B pe3syIbTaTe MapTeHCHTHOrO MpeBpalweHus. B MarHHTHOM Hose pa3BuTHe MapTeHCHTHOTO 
IipeBpalleHuaA HAYHHaeTcA Ip OosIee BLICOKOM TeMIepatype, UTO He MOXKeT HaliTH OOBACHEHHA C TePMOAMHAMMYeCKHX 
mo3sHyHuH. Takum oOpa30m, HaOs1oqamM oOpa30BaHve MapTeHCHTa Hallps»KeHHA B MHKpOOOBeMaX ayCTeHHTa Cc 
(beppOMarHUTHBIM yIOpAOUeHHeM, KOTOPbIe BOCIIPHHUMAIOT IHEPIHio BHEIHETO MOA Yepe3 MarHHTOCTPHKUMOHHbIC 
HallpsKeHua. B ycNOBHAX CBepxIWIaCTHYHOrO ayCTeHHTa Take Hallps9KeHHA OKa3bIBaIOTCA JOCTATOUHBIMH JIA 
MHHUMMpoBaHHA CBHTOBOrO UpeBpalljeHua. OmnpeyeseH MHHMMaIbHO BO3MOX%KHBIM pa3sMep  (plyKTyalHi 
HeyCTOM4MBOCTH peweTKH (1,372 HM). 

O6cystcoenue u 3akiiouenue. Bo3yevicTBue MarHUTHbIM TOJIeM Ip 3akasIke MpHBOAMT K YCHJICHHIO TIpOleccoB 
cBoeoOpa3sHoro MarHHTHOrO paccioeHua aycTeHuta. IIpu tTemnepatypax, Omu3KHx K Hayally MapTeHCHTHOrO 
TipeppalleHua, UMerloljHecd OONaCTH MarHHTHOM HeOAHOpOAHOCTH HakslafbIBaloTcA Ha 9deKTbl OT ABJICHHA 
HeYCTOMYMBOCTH KpHCTaIM4eCKON PeLIeTKH epey, MapTeHCHTHBIM IIpeBpatteHHem. B TemrepaTypHoM uHTepBasie M,- 
M,, Koryla ayCTeHUT IpoOsABIIAeT CBEPXIIaCTHYHOCTS, CYLIeECTBeEHHO OOsIeruaeTca OOpa3z0BaHHe MapTeHCHTa Halps.KeHHA 
B MHKpooObemax aycTeHHTa c (beppoMarHHTHBIM ynopsaqOueHHeM. 


Kunrouesbie CJI0Ba: MapTeCHCHT HallpsKCHHA, MarHHTHOe NOI, CBepXIWIaCTHAHOCTS, HeyCTOM4MBOCTB pewueTKH, CTaJIb, 3aKasIKa 


BuaarogxapHocrn. ABTOpBI BbIparKaroT OlaroyqapHOcTh peleH3eHTaM, 4WbA KPHTHYVeCKaA OWeCHKa IipeACTaBsICHHbIxX 
MaTepHaJIOB HW BbICKa3aHHble MpeTIO%KCHHA NO UX COBCpIICHCTBOBAHHHAO, cnocoOcTBOoBasH 3HaYHTCJIBHOMY IIOBbILICHHFO 
KauyecTBa HacToslmel CTaTbH. 


Asia unTHpopanna: Jlonrayes IO.B., ycrosoiit B.H., Beppuropos 1O.M. 3apoxyqenve MapTeHcHTa HallpKeHHA B COCTOAHHM 
TIpeyMapTeHcHTHOH HeycTowuMBoctH perlierKH. Advanced Engineering Research (Rostov-on-Don). 2024;24(1):58—65. 
https://doi.org/10.23947/2687-1653-2024-24-1-58-65 


Introduction. Before starting the martensitic transformation, instability of the crystal lattice is detected in metals and 
alloys [1] expressed in a specific premartensitic state of the initial lattice [2]. The thermodynamic analysis of the state of 
the initial phase before point M, [3] shows the possibility of the emergence of separate micro regions having their own 
short-range ordering in the arrangement of atoms. Such micro regions retain their atomic order only near the fluctuation 
nucleus. At a distance from it, there are no obvious boundaries between phases, the order gradually degrades. The 
experimentally described micro regions manifest themselves in the form of diffusion scattering effects during X-ray 
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diffraction. The presence of regions with a near atomic order explains the decrease in elastic constants near point M, and 
the convergence of structures [4] involved in the phase transformation. 

In the case where an applied external stress is present, the martensitic transformation can be initiated by deformation 
at a temperature above point M,. The maximum temperature at which the deformation of the initial phase causes the 
formation of martensite is indicated by point Mz. It has been found [5] that a maximum of plasticity is observed in the 
temperature range My — M,, which is due to the manifestation of superplasticity associated with the hardening transition. 
The imposition of a constant magnetic field under these conditions causes the appearance of forced magnetostriction and 
corresponding stresses in the austenite nanodomains having a ferromagnetic order [6], which induces the appearance of 
crystals of the hardening phase called stress martensite (by analogy with deformation martensite) [7]. The imposition of 
a constant magnetic field changes the initial magnetic state of austenite [8], increasing the number and size of 
ferromagnetically ordered clusters, which are possible sites of nucleation of the ferromagnetic a-phase [6]. Thus, in the 
temperature range M,—M,, the integration of magnetic and a number of special structural states occurs, which is 
undoubtedly of interest to investigate. 

This research was aimed at studying the features of martensitic transformation in the presence of a constant magnetic 
field in steel taking into account the phenomena occurring in the premartensitic state. 

Materials and Methods. In the course of the research, samples of steel P6M5 were used, for which there was 
information about the development of superplasticity [5] associated with phase transformation [9]. The chemical 
composition was monitored using a Q8Magellan optical emission spectrometer (Bruker). The average composition of 
samples from one melt is shown in Table 1. 

Table 1 


Average content of elements in samples, % 


To analyze the course of hardening transformation without a field and with the imposition of a magnetic field, the 
potentiometric method of resistometry was used [10], since changes in the alloy structure were reflected in the measured 
values p with a fairly high sensitivity. And although it is difficult to quantify the relationships between the transforming 
phases, the establishment of the start and finish points of the transition can be done quite accurately. The applied method 
is highly sensitive to the appearance of a ferromagnetic phase (during phase transformations of the Ist or 2nd kind), which 
is shown in abnormal behavior on curve p due to electronic interactions at “s” and “d” sublevels [11] determining the 
presence of spontaneous magnetization [12]. 

An S-type thermocouple (PTC) was used as a temperature sensor. The thermocouple was spot-welded in the central 
part of the sample to a single point to eliminate the occurrence of side stresses under separate welding due to the passage 
of current through the sample during measurement. The data was recorded using an analog-to-digital converter L-CARD 
E14—440 and the LGraph2 software package. The sample was heated by passing current from the RNO-250—5 
autotransformer. 

The sample was placed in an interpolar space (Fig. 1). An open-type electromagnet (FL-1) provided the generation of 
a magnetic field with a strength of 1.2 MA/m. When conducting the experiment without a field, the power supply of the 
electromagnet was turned off, and the poles were bridged with an ARMCO-iron plate. 


A 


Fig. 1. Sample with connected copper contacts in the interpolar space of the magnet 
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The experimental data (p(t) and 7(t) in Fig. 2), recorded with a sampling rate of 400 kHz, were numerically 
differentiated and approximated to obtain function dp/d7(7), which reflected the characteristic anomalies when a 
ferromagnetic phase appeared. 

Research Results. Primary dependences p(t) and 7(t) shown in Figure 2 already demonstrate the shift of anomalies, 
corresponding to the phase transformation, to the high-temperature region when a constant magnetic field is applied. 


T, °C p, nOhm-m Y ies ©: p, nOhm-m 
800 800 800 800 
600 600 600 600 
400 400 400 400 
200 200 200 200 
0 2 4 T, min 0 2 4 tT, min 
T T 
a en) === f— 
a) b) 


Fig. 2. Dependences p(t) and T(t) for steel P6M5 under cooling: 
a — without field; 5 — in magnetic field of 1.2 MA/m 


According to the reference data [13], for samples of steel P6MS corresponding to the chemical composition indicated 
in Table 1, temperature of 140°C corresponds to the initiation of the martensitic transformation. This is confirmed by the 
data obtained as a result of differentiation in Figure 3, which illustrate the characteristic features of the analyzed 
transformation. 

The formation of a ferromagnetic phase corresponds to the appearance of a peak at low temperatures (Fig. 3), which 
makes it possible to record the start of the martensitic reaction. In the case of treatment with the imposition of a magnetic 
field, the start of the formation of the martensitic phase is noted at a temperature of 185°C, which cannot be explained 
solely from the point of view of changes in the thermodynamic potentials of the phases. For a magnetic field with a 
strength of 1.2 MA/m, the effect of the equilibrium temperature shift is ~ 4.5°C [14]. It can be argued that the imposition 
of a magnetic field in the temperature range of superplasticity causes the formation of stress martensite higher than the 
starting point of transformation known for this steel as a result of forced magnetostriction in nanodomains with a near 
magnetic order in austenite. In this case, the stresses arising in the austenitic matrix are about 10 MPa [15], which, under 
conditions of the crystal lattice instability, contributes to the initiation of a phase reaction. 
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Fig. 3. Results of differentiation of experimental data obtained under cooling of steel P6MS: 
1 — without a field; 2 — in a magnetic field with a strength of 1.2 MA/m 


The phenomenon associated with the occurrence of instability of the crystal lattice of the initial phase, as it approaches 
the transition point, is reflected in the differentiated dependences in the form of a peak in the high-temperature area of the 
curve. The appearance of this peak under cooling indicates the approach of the superplastic state of austenite. The effect 
of an external magnetic field enhances the magnetic decomposition of austenite [8], contributing to an increase in the 
number and size of ferromagnetic nanoclusters [15], which causes a shift of this peak to higher temperatures. Anisotropic 
fields arise in local regions with unidirectional spins. They introduce elastic distortions due to the fact that the 
magnetization energy differs in different directions — this reduces the stability of y-phase lattice. The average diameter 
of a ferromagnetic nanocluster determined in [8] is about 1.7 nm for a given field strength, which is comparable to the 
wavelengths of electrons. It can deflect them and, accordingly, cause distortions in the curves. The results obtained 
indicate the presence of heterogeneity in the magnetic state of austenite, and that this heterogeneity increases when an 
external magnetic field is applied. 

The transformation under hardening is characterized by a diffusionless shear mechanism of atomic restructuring; 
therefore, external action increases the driving force of the transition. In a state where the lattice of the initial phase 
becomes unstable before point 4, spontaneous atomic displacements occur in local regions (preparing the lattice for ya 
transition), where magnetostrictive stresses can act as a trigger for the start of the transformation. 
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The conditions describing the possibility of fluctuations of atomic displacements in the premartensitic state were 
considered in [16]. The type of function, whose minimum corresponded to the transition state, was derived on the basis 
of diffuse scattering patterns and existing options for restructuring the crystal lattice. Oscillation frequencies and 
anharmonicity parameters were expressed in terms of elastic modulus, whose values determined the adjustment option. 

The analysis suggests that micro regions of atomic displacements are shaped like plates of densely packed planes that 
can experience shear. The expression for estimating the minimum possible fluctuation size H can be derived using matrix 


=), ' 
may (1) 
2c'+Ca4 


where d — interplanar spacing; c’ = 1/2(ci1 — c12) — elastic constant; c11, C12, c44 — elastic coefficients. 

Values for calculation according to formula (1) were selected for the temperature of martensitic transformation in steel 
P6MS. Lattice constant a was estimated taking into account the instructions in [17]. The interplanar spacing was 
determined for the case of a FCC lattice, when the first reflecting planes were a family of planes {111}, having the 


elasticity coefficients: 


maximum interplanar spacing d = “ =21A. Based on the data on the elastic characteristics of steel P6MS5 [9], the 


3 


elasticity coefficients and constant were determined through known ratios [18]. The calculation results based on 
dependence (1) allow us to estimate the minimum diameter of the region where instability H =1.372 nm may occur, 


which is commensurate with the average diameter of a ferromagnetic cluster in austenite. 

Discussion and Conclusion. The following script for the development of martensitic transformation can be assumed. 
Above point M, in y-phase, there are regions with ferromagnetic ordering [15]. If waves of atomic displacements are 
superimposed in or near these places, then, in the presence of an external magnetic field, under such conditions, forced 
magnetostriction is able to affect the fields of elastic forces in the crystal lattice and lower the energy barrier for the 
formation of a critical martensite nucleus. 

The experimental results indicate the existence of magnetic heterogeneity in the state of y-phase, whose degree affects 
the running of the martensitic transformation. The imposition of a constant magnetic field under quench cooling enhances 
the existing magnetic inhomogeneity in austenite. Near point M,, the phenomenon of lattice instability in combination 
with magnetostrictive effects from the field initiates the appearance of stress martensite crystals. Thus, the effect of 
intensification of martensitic transformation under hardening in a magnetic field with a strength of 1-2 MA/m can be 
explained, which is of great importance for the practice of heat treatment of steel. 

The research results showed the following. Exposure to a magnetic field under quench cooling made it possible to 
achieve greater completeness of the martensitic reaction, reducing the amount of residual austenite, whereas the early 
formation of stress martensite provided its longer presence in the area of elevated temperatures, which contributed to the 
flow of tempering processes directly under hardening. 
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Abstract 

Introduction. The problem of improving the manufacturing of high-precision products is currently becoming a key one, 
since the requirements for them are constantly being tightened. Maintaining assembly quality and accuracy is an important 
aspect of manufacturing precision products. Standard approaches to this process do not always have sufficient versatility. 
Existing studies that aim to develop universal approaches, such as end-to-end production design or the application of 
parallel engineering principles, also have a number of shortcomings. These include the fact that the given approaches do 
not fully take into account information about the technology capabilities of a particular production when making design 
decisions, and do not consider the relationship between the manufacturing preparation of machining and mechanical 
assembly industries. That is why studies aimed at developing such universal approaches have high applicability. To solve 
these problems, the authors conceptually developed a set of formalized design procedures for a system of accounting 
requirements for the assembly of high-precision products in the design of machining processes. However, to effectively 
identify the numerous requirements for the assembly of high-precision products (output parameters) and select those that 
cannot be provided by the method of complete interchangeability (critical elements), additional research is needed. The 
research objective is to develop principles for identifying the output parameters of high-precision products and detecting 
critical elements. To achieve this goal, it is required to solve the following problems: to formulate principles for 
constructing generalized surface graphs of high-precision products; to develop standards for classifying output parameters 
and identifying critical ones. 

Materials and Methods. To conduct the research, a high-precision assembly unit was selected — “Stator Package 2”. 
The research was carried out under real conditions of the existing multiproduct manufacture. For this assembly, a 
generalized surface graph was constructed, including information about the nature and sequence of surfaces, requirements 
for the assembly, dimensional tolerances and tolerances of shape and location, with its subsequent analysis. 

Results. This paper presents the results of research on improving the enlarged block of design procedures for analyzing 
requirements for the assembly of high-precision products of the designed system. The paper established the relationship 
between the accuracy of dimensional tolerances and the tolerances of the shape and location of the surfaces of the product 
to which these dimensions belonged. Based on the relationship obtained, an order was determined for the unambiguous 
identification of critical elements. 

Discussion and Conclusion. The application of this technique makes it possible to increase the reliability of the source 
information obtained during the implementation of an enlarged block of design procedures, as well as the validity and 
efficiency of identifying rational manufacturing technologies at subsequent stages of the system implementation, while 
providing the specified quality, accuracy of products, and reducing the complexity and cost of their manufacture. 


Keywords: mechanical assembly production, process design, generalized surface graph, assembly, high-precision 
product, design dimensional analysis 
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AHHOTalna 

Beedenue. WipoOnema coBepuieHCTBOBaHHA TIPOM3BOCTBa BbICOKOTOUHBIX H3eIMi B HacTOAUee BpeMA CTaHOBHTCA 
KJIIOYeBOH, MOCKOJIbKY TpeOoBaHHaA K HAM IOCTOAHHO yxecTouaroTca. OOectleveHue KayecTBa HM TOUHOCTH COOpKH 
ABJIACTCA BAXKHBIM ACIICKTOM IIpOH3BOACTBa BbICOKOTOUHOHM MposykKuMH. CranyapTHble MOAXObI K STOMY TIpoweccy He 
Bcerga OONaqaloT WOCTaTOUHOM YyHUBepcabHOCcTHIO. CyulecTByIOMHe HCCeqOBaHHA, WesIbIO KOTOPBIX ABJIACTCA 
pa3pa0oTKa yHHBepcalIbHbIX MOAXO0B, TaKHX KaK CKBO3HO€ TeXHOJIOTMYeCKOe MpOeKTHpOBaHHe HIM MIpHMeHeHHe 
MIPHHUMMOB WapasiesIbHOU HWKeHEpHOH paspaooTKH, TakwxKe OOaaroT pAOM HeOCTATKOB. K ux YHCILY MOXKHO OTHECTH 
TO, 4TO B HHX He B MONHOM Mepe OOecreyeH yueT HH(POpMallMH O TeXHOJIOFHYeCKHX BO3MOXKHOCTAX KOHKpeTHOrO 
TIpOW3BOACTBa IIpH MPHHATHH KOHCTPyKTOPCKHX pelleHHii H He Y4HTbIBACTCA CBA3b MCKAY TeXHOOrM4eCKOH 
TMOATOTOBKOH MexaHooOpabaTpIBalolero HM MexaHOcOopoyHoro mmpou3BoycTB. MmMenHo nosToMy UcciIeqoBaHHA, 
HallpaBJIeHHble Ha pa3paOOTKy YHMBepCaIbHBIX TOAXO0B, OOaqaloT BbICOKOM akTyasbHocTbIO. J[It pelwleHna 
oOOo3HayeHHBIX TpoOseM aBTOpaMH OBIT KOHIelTyaIbHO pa3paboTaH KOMIVIeKc OpMaJIH30BaHHbIX IPOeCKTHBIX 
TIpolleqyp CHCTeMbI yueTa TpeOoBaHHit K COOpKe BbICOKOTOYHBIX H3AeHM IpH MpOeCKTHpOBaHHH TeXHOJIOrM4eCKHX 
TIpOWeccoB MexaHHyecKoH oOOpadoTKH. OpHaKo IA 3P(PeKTMBHOTO BbIABIICHHA MHOXKeCTBa TpeOoBaHHli K cOopKe 
BbICOKOTOUHBIX M3 eM (BLIXOAHBIX NapaMeTpos) H BbIOOpa M3 HX TeX, KOTOPbIe He MOTYT OBITS OOecrIe4eHbI MCTOOM 
TIOHOW B3aMMO3AMCHACMOCTH (KPHTHUHBIX 3JIEMCHTOB), TpeOyeTcA MpoBeeHHe JOMOHHTeIbHBIX MCCIeqOBaHH. 
Llemb1o JaHHOM padorTsl ABIAeTCA pa3spaOOTKa MIPHHUMMOB BLIABIICHHA BbIXOJHBIX MapaMeTPOB BBICOKOTOUYHBIX H3(euIHi 
WM OMpeeueHHA KPHTHUHBIX 9eMeHTOB. JIA WOCTWMKeHUA MOCTaBMIeHHOM We HeOOXOAHMMO pelIHTb ciezyrouHe 
3aaqH: c(POpMHpoBaTb MPHHUMIMIbI MOcTpoeHuA OOOOUIeHHBIX TpaoB MOBepXHOCTei BbICOKOTOYHBIX v3eHi, 
CPopMupoBaTh CTAaHAapThI KaccHPuKalM BLIXOJHBIX WapaMeTPoOB H BbIABJICHHA U3 HAX KPHTHUHBIX. 

Mamepuanoi u memoooi. [na nposneyenua UccieqOBaHHA Oblia BbIOpaHa BbICOKOTOUHAaA COopouHasd eqMHULa — «lake 
ctaTopa 2». ViccneqoBaHua poBOAWIMCb B PpeaIbHBIX YCJIOBHAX leHCTBYIOWerO MHOrOHOMeHKaTypHOro 
mpou3BoyctTBa. Jia WaHHol cOopKu Ob WocTpoeH oOoOmleHHEIM rpad) MoBepxHocTel, BKIHOUAFOIIMM B cebs 
MHopMalluio O XapakTepe HM MOCIeAOBAaTebHOCTH MOBEPXHOCTeH, MpeAbABIAeMbIe TpeOoBaHHA K COoOpKe, AOMycKH 
Pa3MepoB H AOMycKH POpMBI HM pacHoO*KeHHA C MOCMCAYIOWIUM ero aHasIH30M. 

Pesyibomamel ucciedoeanua. B JjaHHOM CTaTbe TpesCTaBsIeHbl MTOM UccieqoBaHuii M0 COBepIIeHCTBOBaHHIO 
YKpylHeHHOro Os10Ka MpOeKTHBIX Mpoleyyp aHasiM3a TpeOoBaHHi K COOpKe BbICOKOTOUHBIX H3eHM mpoekTupyemon 
cuctembl. B padote ycTaHoOBsIeHa B3aMMOCBA3b M@%KY TOYHOCTHIO JOMYCKOB Ha pa3Mepbl HW AOMycKoB dopMbI u 
pacrouIoxKeHHA NOBEPpXHOCTet H3AeMA, KOTOPbIM IIpHHasiexKaT JaHHble pasmeppl. Ha ocHose MosyYeHHOH B3aHMOCBA3H 
OBL ONpeseeH MOPATOK DJIA OHOZSHAYHOTLO BBIABJICHHA KPHTM4HBIX 3JICEMCHTOB. 

O@cyorcoenue u 3axuo“enue. IIpumMeHeHnve aHHOM MeCTOAMKH MO3BOJIACT MOBLICHTb JOCTOBEPHOCTb HCXOHOM HHPopMallMH, 
MOJYYCHHOH TIpH pearM3allWu yKpyNHeHHoro O0Ka MpOeKTHBIX Mpolesxyp, a TakKe OOOCHOBaHHOCTS U 9(ptpeKTHBHOCT 
BBIABICHHA PallMOHAJIbHbIX TEXHOJOTHM W3TOTOBICHHA Ha MOCIeYIOWMX 9Talax peasu3allMH CHCTeMbI TIpH oOecneyeHHH 
3aJJaHHbIX KavecTBa, TOUHOCTH H3eIHi UH CHWKEHHM Tpy{OCMKOCTH H CeOeCTOMMOCTH HX H3rOTOBJICHHA. 


KsroyeBbie c10Ba: MexaHOCOOpOuHbie MpOH3BOACTBa, TeEXHOJIOrMYecKad MOATOTOBKa IpOH3BOACTBa, OOOOMeHHBIT 
rpad MoBepxHocTeli, cCOopKa, BbICOKOTOUHOE U3eIMe, KOHCTPyKTOPCKHH pa3sMepHbIi aHas3 


BaarojxapHocru. ABTODBI BbIparKaroT OarolapHOcT peak 3a BHHUMAaTeCJIBHOe OTHOIMICHHE K COJep2KaHHtO CTaTbH U 
BbICKa3aHHble UpeIOKCHHA, KOTOPbIC MOSBOJIMJIM MOBbICHTb Ce KAaYeCTBO. 
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Introduction. The production of modern devices and machines requires a clear organization of the assembly process 
with careful manufacturing preparation [1]. This requirement is particularly urgent in the manufacture of high-precision 
products in precision machine tool industry, aerospace industry, military-industrial complex, bearing industry, etc. [2]. In 
this paper, high-precision products mean such products, whose assembly process is complex due to difficulties in meeting 
specified requirements (critical assembly requirements, or critical elements), such as dimensional, dynamic and other 
characteristics, which may result in the need to manufacture a large number of spare parts and kits. This, in turn, causes 
an increase in the volume of work in progress and requires more time and resources to complete the assembly production 
process (PP) [3]. 

Modern technology makes it possible to improve the quality and accuracy of production, reduce material costs, and 
introduce resource-saving technologies. However, today there is no unified concept of an integrated approach to solving 
the above-mentioned problems. When assembling products, classical methods of maintaining accuracy are often used, 
which in most cases are not universal. They depend on the type of production and are often economically unjustified. In 
such cases, it is necessary to develop unique assembly approaches that take into account the specifications of a particular 
product and provide the required accuracy without unreasonable costs. In addition, it should be borne in mind that the 
traditional PP design of mechanical processing and assembly is subjective. This is primarily due to time constraints and 
physical capabilities, which make it impossible for a person to compete with modern computing technology [4]. To solve 
these problems, it is needed to develop unique approaches that combine various aspects of production and assembly of 
products, as well as taking into account economic factors. Only a comprehensive solution can provide achieving optimal 
results and effectively cope with the tasks set. 

Currently, a system of automated process planning (SAPP) has been conceptually developed [5]. Nevertheless, the 
weakness of this automated system is the lack of correlation between the multiple-alternate design of the manufacturing 
process for parts included in the assembly unit and the requirements for the technology and accuracy of assembly of a 
high-precision assembly unit, based on which it is possible to efficiently solve the problem of reducing the volume of 
work in progress and the cost of assembly. At the same time, the correlation between design and manufacturing 
preparation of production has not yet been fully worked out. As a result, research in this area is urgent for the modern 
development of mechanical assembly production systems. The objective of the article is to develop principles for 
identifying the output parameters of high-precision products and determining critical elements. To do this, it is needed to 
formulate the principles of constructing generalized graphs of surfaces of high-precision products, determine the standards 
for classifying output parameters and identifying critical ones from them. 

Materials and Methods. To solve the above problems, a system with account for the requirements for the 
assembly of high-precision products in the design of manufacturing processes of mechanical processing was 
developed (SRPPMP) [6]. However, to provide a more comprehensive integration of this system into the SAPP 
structure, it is required to actively look for options of upgrading the development strategies used. This will allow us 
to move on to the analysis of the manufacturability of products and provide a more effective implementation of the 
system into the production structure. 

As part of the research to improve the methodology of formalization of the block under consideration, strategies for 
performing the procedures of this block were formed, and the one that provides identifying the maximum number of 
critical elements was selected. This strategy is a combined one and matches up the advantages of the other approaches. 
In addition, for a more accurate calculation of dimensional chains, a special size indexing system has been introduced, 
which is most effective. The specified indexing is considered in detail in [7] and is not fully cited in this paper. 

According to this strategy, detailed segmentation of high-precision products is performed, which provides determining 
the required parameters for the formation of a set of assembly requirements (Mrc), carrying out a design dimensional 
analysis and, based on its results, identifying critical elements. To implement the described procedures, first, it is needed 
to identify those elements from the set Mr.c that satisfy the following condition: 

on ae 
TBA) <2 TAK (l) 
where b — assembly requirements (closing links); A — component links. 

The above condition makes it possible to identify those elements that cannot be provided by the method of complete 
interchangeability, i.e., critical elements. The specified elements are included in the corresponding subset My. The 
remaining elements are provided by the method of complete interchangeability. They form subset Mrz. and are no longer 
considered. 
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M; =—~. (2) 
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Next, it is necessary to determine methods for maintaining accuracy for each critical element and distribute them 
into appropriate subsets: Mu.3. — critical elements provided by the method of incomplete interchangeability; 
Mrz. — critical elements provided by the method of selective assembly; Mu.n., Mp.z. — critical elements provided 
by compensation. 

Research Results. For the distribution of critical elements into groups, a technique based on the construction of a 
generalized graph of surfaces of high-precision assembly units and nodes has been developed, and distribution rules have 
been formulated. 

The classification of modules of working and connecting surfaces of parts and assembly units (PAU), proposed by 
Professor Bazrov B.M. [8], is shown in Figure 1. 
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Fig. 1. Modules of working and connecting surfaces of PAU [8] 


Consider the principles of constructing a generalized graph of surfaces using the example of the assembly unit (AU) 
“Stator Package 2”, which is part of the high-precision ““Gyromotor’” product. Figure 2 shows an axonometric projection 
of the AU “Stator Package 2” and a schematic representation of this AU, indicating the basic assembly requirements, 
shape and location tolerances imposed on it, as well as the PAU of which it consists. 
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Figure 3 shows a sketch of the AU “Stator Package 2”, on which the surfaces of the composite PAU are numbered, 
and a code is indicated for each surface in accordance with the requirements of the SAPP [9] (Table 1). 

Based on the information obtained about the nature and sequence of the PAU surfaces, a graph of the AU surfaces 
“Stator Package 2” is constructed. In this graph, the surfaces of the PAU included in the considered AU are indicated. For 
the surfaces of the PAU, parallel to the axis of rotation of the AU (marked “01” in the graph), the connection with this 
axis is conditionally shown. Additionally, the PAU surfaces that are in direct contact in the assembly unit are marked 
(e.g., surfaces 1.7 and 3.3 in Figure 4a). Further, the assembly requirements (dimensional chains) posed on this AU are 
presented on the specified graph, indicating the closing links and the directions of the constituent links. It is important to 
note that the components and closing links going to or from the group of contacting surfaces are assembly dimensions. 
Also, the graph indicates the basic shape and location tolerances required for the AU assembly (Fig. 4b). The assembly 
requirements and basic tolerances of shape and location are taken from Figure 25. As a result, a generalized graph of AU 
surfaces “Stator Package 2” is formed. 

Thus, the presented graph provides obtaining the required data on the dimensional relationships between individual 
surfaces (groups of surfaces) of PAU, as well as the surfaces of the mating PAU in contact with each other. The obtained 
data form the initial information for the further stages of the implementation of the SAPP procedures. 
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Table | 
No. of PAU 
° Cod T 
surfaces aces as 
Part “Bushing” 
1.1 0311 Outer end 
1.2 0221 External chamfer 
1.3 0111 External cylindrical 
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1.4 0411 Curved external 
1.5 0211 Outer cone 

1.6 0111 Outer cylindrical 
1.7 0311 Outer end 

1.8 0411 Outer curve 

1.9 0111 External cylindrical 
1.10 0211 Outer cone 

1.11 0311 Outer end 

1.12 0111 External cylindrical 
1.13 0311 Outer end 

1.14 9211 Internal cone 
1.15 9111 Internal cylindrical 
1.16 9211 Internal cone 
1.17 9111 Internal cylindrical 
1.18 9114 Cylindrical hole at an angle to the axis of rotation 
1.19 9411 Curvilinear internal 
1.20 9311 Inner end 

1.21 9411 Curvilinear internal 
1.22 9111 Internal cylindrical 
1.23 9311 Inner end 

1.24 9111 Internal cylindrical 
1.25 9411 Curvilinear internal 
1.26 9111 Internal cylindrical 

Part “Axle” 

2.1 0311 Outer end 

2.2 0211 Outer cone 

2.3 0111 External cylindrical 
2.4 0111 External cylindrical 
2.5 0411 Outer curve 

2.6 0221 External chamfer 
2.7 0111 External cylindrical 
2.8 0311 Outer end 

2.9 0411 Outer curve 
2.10 0111 External cylindrical 
2.11 0211 Outer cone 
2.12 0111 External cylindrical 
2.13 0221 External chamfer 
2.14 0311 Outer end 
2.15 0111 External cylindrical 
2.16 0111 External cylindrical 
2.17 0211 Outer cone 
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2.18 0311 Outer end 
2.19 9211 Internal cone 
2.20 9111 Internal cylindrical 
2.21 9111 Internal cylindrical 
2.22 9411 Curvilinear internal 
2.23 9111 Internal cylindrical 
2.24 9311 Inner end 
2.25 9111 Internal cylindrical 
2.26 9411 Curvilinear internal 
2.27 9311 Inner end 
2.28 9411 Curvilinear internal 
2.29 9114 Cylindrical hole at an angle to the axis of rotation 
2.30 9211 Internal cone 
Part “Spacer” 

4.1 0311 Outer end 

4.2 0111 External cylindrical 

4.3 0311 Outer end 

4.4 9211 Internal cone 
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Fig. 4. Generalized graph of AU surfaces “Stator Package 2” 


The next step is to establish the relationship between the assembly requirements and AU surfaces. In the 
presented example, the assembly requirements are evaluated, interconnected through the surface of the part 
“Bushing” 1.13 “Outer end” (Fig. 5). 

For the qualitative distribution of critical elements into groups and the detection of those that cannot be provided by 
the method of complete interchangeability, it is necessary to conduct a correlation analysis of the accuracy of dimensional 
tolerances between the surfaces of the PAU and the specified tolerances of the shape and location of these surfaces. 

The performed studies included the analysis of the interfaces of working or connecting surfaces predominantly found 
in the PAU of high-precision products and their correlation with shape and location tolerances (total shape and location 
tolerances), which are most characteristic of these combinations. 

Based on the conducted research, an order was determined for the unambiguous identification of critical elements that 
cannot be provided by the method of complete interchangeability. As a result, specific dimensional chains and design 
parameters of the unique parts included in them were identified: 

1. The critical covering and covered dimensions of the mated PAU must be at least 9 accuracy standards. 

2. Linear dimensions must be at least class f (according to the requirements of GOST 30893. 1-2002!). 

3. The tolerances of the shape and location of the mating surfaces must have high relative geometric accuracy 
(C — 12% of the size tolerance for cylindricity, roundness, and longitudinal profile, 25% — of the size tolerance 
for the rest) [10]. 

4. The tolerance of the angular dimensions of the cones should be no coarser than AT 10 (the average value from the 
proposed range in GOST 8908-81? was selected). 

The performed approbation of the developed graph construction techniques and the formed rules on the example of 
the analysis of the requirements for the AU “Stator Package 2” under the conditions of the current multiproduct 


' GOST 30893.1-2002 (ISO 2768-1-89). General Tolerances. Limit Deviations for Linear and Angular Dimensions without Tolerance Indications. 
URL: https://files.stroyinf.ru/Data2/1/42948 16/42948 16276.pdf (accessed: 29.11.2023) (In Russ.). 

2 GOST 8908-81. Standard Angles and Angle Tolerances. Basic Norms of Interchange Ability. URL: https://docs.cntd.ru/document/ 1200011833 
(accessed: 29.11.2023) (In Russ.). 
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manufacture showed their operability. Identified critical assembly requirements that do not meet the conditions of 
complete interchangeability include: 

— permissible value of the deformable (rolled) part of the “Bushing” workpiece obtained during the assembly of the 
AU “Stator Package 2”; 

— permissible value of the chamfer protrusion of the “Axle” part in the AU “Stator Package 2”. 
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of the Axle” part in AU “Stator Package 2” 


Fig. 5. Requirements for assembly of AU “Stator Package 2” with indication 
of surfaces between which the corresponding chain links are located 


As a result of the implementation of this enlarged block of SRPPMP, initial information is formulated, including 
numerous assembly requirements (closing links of dimensional chains) that cannot be provided by the method of complete 
interchangeability. The dimensions of the PAU are the constituent links of the dimensional chains of the specified set and 
the set of SRPPMP to which these closing and constituent links belong. In the example under consideration, such a set 
will consist of “Bushing”, “Axle”, “Stator Package 1”. The dimensions of the part “Spacer” are not involved in the 
calculation of critical assembly requirements. At the next stages of the implementation of the SRPPMP for the received 
set of PAU in the SAPP, all possible options for PP are generated, and their in-process dimensional analysis is carried 
out. Next, the most rational options are selected from the resulting set of PP using the multicriteria optimization method. 
To carry out the selection procedure, a system of criteria was developed. The criterion for obtaining the maximum number 
of assembly kits is the basic from them. The technique of generation procedures, in-process dimensional analysis, and 
selection is described in detail in [6]. It should also be noted that the improvement of the methodological apparatus of the 
SRPPMP will allow us to proceed to the assessment of the production manufacturability of the product [11]. 

Discussion and Conclusion. The paper proposes the principles of constructing generalized graphs of high-precision 
AU surfaces and nodes, develops standards for classifying high-precision output parameters into appropriate groups and 
identifying critical elements of high-precision products from them. The described technique creates conditions for the 
complete formalization of the considered design procedure and the clarity of the graphical results of dimensional analysis. 
This technique makes it possible to increase the reliability of the initial information obtained during the implementation 
of an enlarged block of design procedures, as well as the validity and efficiency of identifying rational manufacturing 
methods for PAU at subsequent stages of the SRPPMP implementation, while providing the specified properties, accuracy 
of products, and reducing the complexity and cost of their manufacture. 

The improvement of the scientific principles of the implementation of this enlarged block of design procedures of the 
SRPPMP increases the level of digital transformation of design and manufacturing preparation, and, along with the issues 
of automation of design, production and assembly, creates prerequisites for the implementation of a systematic approach 
to assessing the manufacturability of the products. 

The next stage in the development of the proposed study is planned to supplement the generalized graph with 
information about the design and processing bases of the PAU, and the nature of combinations of mating connecting and 
working surfaces according to the classification of Professor B.M. Bazrov. The introduction of the specified information 
into the generalized graph will make it possible to finally form the methodological support for this enlarged block of 
design procedures. 
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Abstract 


Introduction. Increasing the durability of gas turbine engine (GTE) blades is achieved through the use of special 
protective coatings on their surface. For the development of such coatings, the basic source information is the geometric 
profile of the blade section. To transfer a given blade cross-section profile to the appropriate CAD/CAM system or 
engineering analysis package, parametric modeling methods are used to automate this operation. However, the known 
approaches to creating a parametric model of a blade profile are not without a number of disadvantages, and a generally 
accepted method for creating it does not currently exist. The research was aimed at creating a technique for approximating 
the profile of gas turbine engine blades, convenient for use in the subsequent analysis of the operating conditions of 
special coatings on the surface of the blades. 

Materials and Methods. When constructing parametric models of the profile of gas turbine engine blades, a method based 
on the orthogonal Legendre polynomials was used. This made it possible to provide high accuracy of approximation and 
construction of a continuous mapping for the parameters of the blade profile approximation. A Python application was 
created for automated processing of source profiles. It provided the calculation of the coefficients of approximating 
polynomials for the contour lines of the blade, visualization of the calculation results, and creation of a dxf file based on 
the points of approximating functions to transfer it to the CAD system. Next, geometric models of blades were used to 
solve the problem of a stationary aerodynamic flow around a blade. The results of solving this problem were used to study 
the effect of the blade profile on its cooling in an aerodynamic flow. 

Results. As an example, three options of blade profiles belonging to different types of GTE were considered. It was shown 
that for all three studied profiles, the proposed technique provided obtaining parametric models that maintained high 
accuracy in constructing approximating lines, which was confirmed by the values of the determination coefficients close 
to unity. To illustrate the possibility of using the obtained models, examples of solving the gas dynamic problem with a 
potential flow around a blade in a stationary aerodynamic flow were given. The distributions of pressure and temperature 
on the surface of the blade were calculated using the finite element method. 

Discussion and Conclusion. The calculation results show that the proposed technique of approximating the profile of the 
GTE blade, based on the use of orthogonal polynomials, is a convenient tool to automate the creation of a geometric model 
of the blade and compare different types and profiles of blades, solving the corresponding gas dynamic problems. At the 
same time, for a given blade profile and GTE operating conditions, it is possible to obtain the distribution of temperatures 
and forces acting on the surface of the blade, which is required for predicting the durability of special coatings. 


Keywords: blades of gas turbine engines, section profile, approximation, orthogonal polynomials 
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AHHOTalna 

Beedenue. Wospnuenue JONTOBeYHOCTH JOMaTOK ra30TypOHHHBIX ABuratTene (IT) qocTuraetca 3a cueT IpHMeHeHHA 
CHeUMaIbHBIX 3all[MTHbIX IOKPbITHH MX NoBepxHocTu. Jt pa3paOoTKH TaKHX MOKPbITHM OCHOBHOM HCXoOZHO 
MHopMallvet ABJIAeTCA TeOMeTPHYeCKHH Upodurs cevenua NoMaTKH. Ja Nepeqaun 3ayaHHoro Mpodusa cevenna 
louaTKH B cooTBeTcTByrollyIo CAD/CAM-cuctemy wim lakeT WHXKeHepHOrO aHasIM3a MIPHMCHAIOTCA MeTOIbI 
HapaMeTpH4ecKOrO MOJeJIMpOBaHHA, MO3BOJIAIOUIMe ABTOMATH3HpOBaTb JaHHy!O Orepaluio. OHako M3BeCTHBIC 
MOAXObI K CO31aHHIO MapaMeTpH4ecKON MOE NpOPUA WOMaTKH He JIMIMCHbI pia HeAOCTATKOB, H OOMIeNpHHATON 
MeTOAUKM ee CO3aHHa B HacTosuee BpemMa He cyujecTByeT. Llemb1o WaHHOM paOoTb! ABJIAeTCA CO3qaHHe MeCTOAMKH 
almpokcuMalun mpoduna sonatoK CTI, ynoOuol Wa Mcnob30BaHHA TPH MOCMeAyIOUeM aHasIM3e YCOBHM paOorTsl 
CIeMMasIbHBIX MOKPbITHH MOBepXHOCTH JIOMAaTOK. 

Mamepuanoi u memoooi. pu noctpoenuu MapaMeTpHyeckux MOJeseH MpodpusaA JOMAaTOK ras0TypOMHHBIX JBHraTesen 
aBTOPaMH MCHOIb30BaH MCTON, OCHOBaHHbIM Ha IPHMCHeCHHM OPTOTOHAJIbHBIX MOJMHOMOB JlexaHypa. ITO WO3BOIMIIO 
OOecIe4YHTb BbICOKYIO TOYHOCTS alllIpOKCHMallMH WM MOcTpoeHHe HelpepbiBHOrO oTOOpaxKeHHA JIA MapaMeTpoB 
allmpokcuMalun npopusa TomaTKH. J[1t aBTOMaTH3HpOBaHHOM OOpaooTKH HCXOAHbIX Ipoduse co3aHo IpH0xKeHHe 
Ha s3bIke Python, mosBosAIOlWee BEIMHCIATS KOOPPUUMEHTH! alMpOKCHMUPyIOWHX MOMMHOMOB JIA JIMHUM KOHTypa 
JIOMAaTKU, BH3YaIM3MpOBAaTb pe3yJIbTaTbI pacueta H CO31aBaTb 10 TOUKAM alllIpoKCHMUpylowMx dpyHKunl dxf-datin wa 
Tlepegqaun ero B CAD-cnctemy. JJaee reoMeTpHueckHe MOJCJIM JIOMATOK HCHONL3OBAaIIN JIA pewleHHaA 3aaun OOTeKAaHHA 
JOMWAaTKH CTalMOHapHbIM adposMHAaMHYeCKUM TOTOKOM. Pe3yJIbTaTbI PeWIeCHHA 3ZTOM 3aadH UCTIONb30BaHbI IIpv 
MCCIe{OBAHHH BIMAHHA IIpOpUsA JOMaTKH Ha ee OCTLIBAHHE B AaIPOMHAMMYECKOM HOTOKE. 

Pe3yibmamot uccnedosanua. B kadecTBe IpHMepa paccMOTpeHbI Tp BapHaHTa Mpopuset OMAaTOK, OTHOCALNUXCA K 
pasHpim TunamM ITJ]. Iloxa3aHo, 4ro fd BCeX TpeX H3Y4eHHBIX Mpodusei mpeqOKeHHAA MeTOAMKa MO3BOIAeT 
TOy4aTb MapaMeTpHyeckHe MOJeIH, OOeCHeYMBaIOMIMe BbICOKYIO TOUHOCTS MOCTPOeHHA alINIpOKCHMUpYIOWINX JIMA, 
4TO NOATBepxRyaeTcA OH3KHMH K eAMHUWe 3HaxyeHHAMM KOSdduuMeHTOB WeTepMuHauuH. JIA wisrOcTpayHn 
BO3MOXKHOCTH HCIHOJIb30BaHHA MOJYYCHHBIX MOJeNeH MpHBeAeHbI IpHMepb! peuleHHA 3aa4H ra30BOl AHHaMUKU TIpu 
TMOTCHIMaIbHOM OOTeKaHHH JIOMATKH B CTalMOHapHOM a3poyqMHamMMyeckom noTOKe. MeToOM KOHeYHBIX 3JIEMeCHTOB 
paccunTaHbl paciipeyeeHuaA JaBIeHHi M TeMiepatyp Ha MOBepXHOCTH JIOMAaTKH. 

O6cyotcoenue u 3aksu0o4enue. PesyibTaTbl BEIMMCICHHH NOKa3asIH, YTO TIpewiaraemMasd MeTOAMKa allmpoKCHMallMH Mpopwa 
nonatku IT], ocHopaHHat Ha HCMOb30BaHHH OPTOTOH@JIbBHbIX MHOFOYWICHOB, ABJIAeTCAH yOOHbIM HHCTpyMeHToM, 
TIO3BOJIAIOLIMM ABTOMATH3HpOBaTb CO31aHe reOMeTPHYECKOM MOJ{eJIH JIOMATKH H IpOBOAUTb cpaBHeHHe pa3JIM4YHbIX THIOB 
MW Mpoduset JoMaTOK, peliiat COOTBeETCTBYIOMIHe 3aayH Ta30Bol AHHaMUKU. IIpu 3TOM AWA 3ayjaHHOro UpopusA OMAaTOK 
MW ycnoBuli padots! IT] MoxkHo noslyuuTb paciipeyeseHue Temiepatyp UM ycuIMi, JevicTByIOWIMX Ha MOBeEpXHOCTH JIoMaTKH, 
4TO HeEOOXOAHMO JVI NPOTHOSHpOBaHHA JOJITOBCYHOCTH CIeLMaIbHBIX MOKPBITHH. 


Ksnrouesble CJI0Ba: jOnaTKH Tra30TYPpOHHHBIX Buratelen, pou ceqeHHaA, allupOKCHMalHaA, OpTOrOHaJIbHble 
MHOrouWJICHbI 


BaarogxapHocru. ABTODBI BbIpaxKaroT OaroapHOcTB peak U peweH3eHTaM 3a BHHMAaTeCJIBHOeC OTHOIMICHHE K CTAaTbe 
HW YKa3aHHble 3aMeuaHhA, YCTpaHeCHHe KOTOPBIX MO3BOJIHJIO MOBbICHTb Ce KadeCTBO. 


Aaa WaTnupoBanna. Conosses M.E., lynesa 1O.H., bangaes C.JI., bangaes JI.X. Anmpokcumana npopusa slomaToK 
Ta30TypOHHHBIX jBuraTelelti. Advanced Engineering Research  (Rostov-on-Don). 2024;24(1):78—-87. 
https://doi.org/10.23947/2687-1653-2024-24- 1-78-87 


Introduction. Turbine blades are the most highly loaded parts of a gas turbine engine (GTE), the durability of which 
mainly determines the engine overhaul time [1]. One of the ways to increase the durability of turbine blades is to apply 
special protective coatings on their surfaces [2]. The coating is a complex composite structure. To select materials for it 
and determine its optimal geometric characteristics, it is necessary to know the operating conditions of the part: the 
distribution of temperatures, pressures, and shear stresses on the surface of the blade. 
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The highest requirements are posed on the shape of the blades, as well as on the selection of materials, their 
manufacturing technique, and special coatings [3]. When designing special coatings, the basic initial information is the 
geometric profile of the surface of the part, since the parameters of the coating manufacturing process are set on its basis. 

Among the characteristic elements of the profile, the following can be distinguished: chord — distance between the 
most distant points with a horizontal arrangement of the blade, like an airplane wing; suction side — upper part of the 
profile; pressure side — lower part of the profile; leading and trailing edges [4]. 

The shape of the blade cross-section profile is selected both on the basis of experimental studies of pilot units and 
actual engines [5], as well as on the results of numerical experiments on models [6]. In the course of numerical simulation, 
we study the task of creating a geometric model of the blade based on a given cross-section profile and transferring it to 
the appropriate CAD/CAM system [7], or an engineering analysis package [8], in which the parameters of the physical 
properties of the materials of the parts are set, the boundary conditions of the problem are determined, and calculations 
are performed. 

Traditionally, when creating a blade drawing, its profile is described by a set of circular arcs, which is then transferred 
point-by-point to the module for developing a geometric model and generating a grid of finite elements. The disadvantage 
of this approach is that when building each new model, you have to manually create its profile in the graphical editor. In 
this regard, methods of parametric modeling of the blade cross-section profile are proposed to automate this operation. 
One of the most well-known models is the nine-parameter RATD (Rapid Axial Turbine Design) model described in [9]. 
At the same time, the practice of using RATD has revealed certain disadvantages inherent in this model, specifically, the 
inconvenience of using it when optimizing the geometry of the profile, as well as insufficient accuracy. On this basis, the 
authors [10] proposed a modified version of this model, which included a larger number of parameters. Methods of 
parametric modeling of blade cross-section profiles are also actively used in domestic practice. Thus, in [11, 12], 
algorithms for the automated construction of the blade profile using approximation by curves from a set of parabolas and 
Bezier curves of the second order were proposed. In [13], a method for designing a grid of GTE profiles based on these 
algorithms was described. The use of parametric models of blade profiles in solving optimization problems was 
considered by the authors in [14, 15]. 

At the same time, it should be noted that the existing methods of parameterization of the blade profile are not without 
a number of disadvantages. The main disadvantage of traditional blade profile parameterization schemes based on the use 
of second-order curves is that different functions have to be used for different parts of the workpiece to maintain the 
accuracy of the description. At the same time, these functions do not form an orthogonal system, and the values of the 
parameters determined by the approximation of the existing profile by the least squares method turn out to be correlated. 
This leads to the fact that it is not possible to construct a continuous mapping for a set of parameters when the system 
approximates the cross sections of the three-dimensional surface of the blade of double curvature. 

It should also be noted that there is currently no universal technique for parametric modeling of GTE blade profiles, 
and the selection of a specific technique depends both on the goals of modeling and on the characteristics of the specific 
type of engine being designed. However, the forms of input information can vary significantly. It can be a geometric 
model in the form of a file in the format of one of the CAD systems, or just a set of points from a database of profiles 
of the type [16]. 

Thus, an additional requirement for the parametric modeling technique of the blade profile, due to the purpose of the 
model, is versatility with respect to the format of the source data. 

This research was aimed at the creation of the technique for approximating the profile of the GTE blades, free from 
the above disadvantages and convenient for use in the subsequent analysis of the working conditions of special coatings 
on the surface of the blades. This was demonstrated by the example of solving the problem of gas dynamics for three 
different blade models, whose profile was parameterized through the proposed technique. 

Materials and Methods. The authors used an option of parametric models of the blade profile based on the orthogonal 
Legendre polynomials [17]. A Python application was created for automated processing of source profiles. This made it 
possible to calculate the coefficients of approximating polynomials for the contour lines of the blade, visualize the 
calculation results in the form of graphs of initial points and approximation curves for the suction and pressure sides of 
the blade, save the initial points of the profile, an array of coefficients of approximating polynomials and the coefficient 
of determination of the model in the database, and create a dxf file based on the points of approximating functions to 
transfer it to the CAD system and an engineering analysis package. 

The constructed two-dimensional geometric models of the blades were used to solve the problem of a stationary 
aerodynamic flow around the blade in the approximation of potential flow. The velocity distributions calculated as a result 
of solving this problem in the flow were used to study the effect of the blade profile on its cooling in the aerodynamic flow. 

Calculation method. The orthogonal polynomial systems used as the basis of the proposed method provide linear 
independence of the approximation coefficients and are devoid of a disadvantage that causes difficulties in constructing 
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models that are continuous in parameters. The degree of polynomials is selected to be high enough — of the ninth order, 
so that the set of two polynomials for two parts of the profile (suction and pressure sides) is twenty parameters, which is 
sufficient for approximating very complex profiles. The type of the approximating function and the method of calculating 
estimates of its coefficients are described below. 

We denote the approximating line of the blade profile using function (x). The expression of this line in the case of a 
number of Legendre polynomials of the ninth degree is represented as: 


y(x)=)o ails(x), (1) 


where a; — coefficients that are parameters of the model; L(x) — Legendre polynomials calculated by formulas: 


L, =(3x?-1)/2; 
L =(5x?-3x) /2; 
Ly = (35x* -30x? +3)/8; 
is = (63x3 - 70x? +15x)/8; (2) 
Lo = 231/16x° —315/16x* +105/16x? — 5/16; 
Ly = 429/16x7 —693/16x° +315/16x? —35/16x; 
Lg = 6435/128x* —3003/32x° +3465/64x* —315/32x7+35/128; 


12155 6435 9009 1155, 315 
128x° 32x? 64x° 32x? 128x- 


Let the coordinates of the initial blade profile be given by a set of points yf = y°(x;),i=1,...,N. Substituting 


values x; in functions (2), we obtain matrix X of size N x 10. The coefficients of the approximating polynomial for the 
given blade profiles were found by the least squares method from the minimum condition of the sum of squared deviations 
of the given values yf and calculated through the regression equation (1): 


A N c 2 i 
a=) (9 y(ai, x:)) > min. (3) 
The estimates of coefficients a; were calculated using the well-known regression analysis formula: 
-1 
a=(X"X) XTy, (4) 
where a — notation for the vector of coefficient estimates; y — for the vector of points of given profile y/. 


The approximation accuracy was estimated by the determinacy coefficient of the model (the coefficient of 
determination R?), which, with a good approximation, should be close to unity: 


way Dlv)-yF) 


, (5) 


where y° — average value of y/. 

In design practice, it is customary to describe the leading and trailing edges of the contour of the blades with the radii 
of circles. In this regard, the coordinates of the points of the approximated profiles between the arcs of the circles of the 
leading and trailing edges were considered in this paper. For the purposes of approximation, the contour of the blade in 
question was assumed to be horizontal, so that the blade chord line coincided with the abscissa axis of the Cartesian 
coordinate system. To preserve the orthogonality condition of the models, size normalization is mandatory. In this regard, 
all dimensions were normalized by the length of the blade chord, so that the abscissas of the profile points lay in the range 
of 0.1. This provided the orthogonality of functions (2) with coefficients calculated by formulas (4). 

Research Results. Table 1 shows the coefficients of approximating polynomials for blade profiles belonging to 
different types of GTE to study the heat fatigue of a thermal-protective coating consisting of two metal and ceramic layers: 
profile of the compressor blade C8626 [13], profile of the NASA blade of a highly efficient high-pressure GTE [18], and 
profile of the simulator blade used in [19]. 
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Coefficients of approximating polynomials for blade profiles from different sources 


Table 1 


Coefficients [13] [18] 

suction side | pressure side suction side pressure side suction side pressure side 
do 142.4250 —123.0993 —1275.6471 717.8336 —1043.6707 527.0915 
a —374.3446 324.9875 3328.6126 —1868.9670 2719.5055 —1369.9895 
ar 477.9332 417.7186 —4181.1239 2338.1048 —3406.1466 1707.9198 
a3 —444 4554 392.7204 3794.0526 —2108.4202 3077.0694 —1531.8025 
as 324.6819 292.0489 —2684.7257 1479.0172 —2164.2992 1066.1471 
ds —189.2284 174.3673 1502.4170 —818.4219 1201.7203 —583.7501 
a6 86.8260 —82.4103 —654.4877 351.5869 —518.2516 247.3087 
aq —30.0775 29.5457 211.8930 —111.9355 165.6344 —77.2958 
dg 7.1533 —7,3024 —46.0284 23.8450 —35.3931 16.0493 
ay —0.9024 0.9652 5.1021 2.5861 3.8411 —1.6752 
R2 0.9999 0.9993 0.9896 0.9990 0.9982 0.9994 


Figure | shows the starting points of these profiles and the corresponding approximating lines. 
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Fig. 1. Initial points and approximating lines of blade profiles presented in Table 1, 
sources of blade profiles: a — [13]; b — [18]; c — [19] 
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As mentioned above, the aim of parameterization in this paper is to approximate the profile by a continuous function 
followed by modeling the distribution of external parameters over the surface of the blade. To illustrate the possibility of 
solving this problem, we consider a simple two-dimensional model of a blade in a stationary aerodynamic flow, in the 
approximation of potential flow. In this case, the distribution of the velocities of flow v(x, y) is described by the equation: 


ve , (6) 


where function y is found from the solution to the Laplace equation: 
Ay = 0, (7) 
which in this case corresponds to the incompressibility condition (V v= 0) and vortex-free flow (V Xv= 0). 
The boundary condition for function y will be equal-zero velocity in the normal direction of the blade surface, which 
means that w is constant on this surface. 
We denote the profile line of the blade surface by S and search for a solution to equation (7) in region Q, bounded by 
the outer part with respect to S and the inner part with respect to circle C of a sufficiently large diameter, compared to the 


length of the blade chord, so that the flow is homogeneous on surface C. That is, the conditions are set at the boundary of 
the domain 0OQ=SUC: 


WV] 5 =O. W[6 = Vex, (8) 
where Vv... — uniform flow velocity vector. 


After calculating the velocities, the pressure distribution can be found using the Bernoulli equation: 


2 
v 


2 
where p — gas density; B — Bernoulli's constant, which can be set equal to ps + ae where p. — pressure at an 


infinite distance from the blade. 
The variational formulation of problem (7) looks like this: 


[vvvw=o. Vwe H)(Q). (10) 
Q 


The solution to this problem was carried out by the finite element method. For comparison, the calculation was 
performed for three variants of blade profiles, shown in Table 1, with the same chord length I. 

The obtained results were further used in solving the problem of cooling the blade in an aerodynamic flow. Such a task 
makes sense for thermocyclic testing of coatings on simulator blades. In this case, the problem statement looked as follows: 


0,7 -V-(KVT)+v-VT =0, (11) 
where 7 — temperature; t — time; « — thermal diffusivity coefficient. 
Initial temperature distribution: 

T(t =0, x)=Tp(x) (12) 
was assumed to be as follows: inside the blade, a constant temperature was 400°C, ambient temperature was 0°C. The 
boundary conditions corresponded to the absence of heat flow on the outer contour C: 

ar 
on 


= 0, tle cage (13) 


Cywen>O0 
The distribution of flow velocities in problem (11) was calculated when solving problem (10), while the velocity on 


contour C was assumed to be v=10[//s] (length dimension in units of chord length). The thermal conductivity 


coefficients of the blade and gas were assumed to be 0.1 and 0.01 [! hs s|, respectively. 


The program code for solving problems (10), (11) was written in the input language of the universal finite element 
package FreeFem++ [20]. 
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Figure 2 shows the pressure and temperature distributions calculated as a result of solving problems (11), (12), when 
the blade cools down for 25 seconds in the aerodynamic flow for the three studied blade profiles shown in Figure 1. 
Table 2 shows the calculated values of the minimum relative pressures near the blade surfaces and the maximum 
temperatures of the blades. 


°) 


Fig.2. Calculated distributions of pressure (p, left) and temperature (T, right) when modeling the flow and cooling in the flow of 
the blades of the studied profiles (profile designations — as in Fig. 1). Flow vector direction angle with blade chord direction was 10° 


Table 2 


Minimum relative pressure near the surface of the blade 
and maximum temperature of the blade after 25 s cooling 
in the aerodynamic flow 


Parameters Blade profile 
fs [18] [19] 


Tris (t= = 25), 156.9 p52 | 2 | 2834 4 
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Discussion and Conclusion. Judging by the data presented in Table 1, the coefficients of determination R? for all 
studied blade profiles were close to unity, despite the fact that the shapes of the profiles differed significantly. This 
confirms the good quality of the approximation and is also illustrated by the data in Figure 1. Thus, the proposed method 
makes it possible to construct approximating blade profile lines for various types of GTE with high accuracy. 

The paper also illustrates the possibility of using approximating functions under calculating the distribution of 
pressures and temperatures on the surface of the blade on the basis of solution to the problem of gas dynamics. Further, 
this distribution can be used to optimize the structure of special blade coatings, such as heat-protective coatings. The data 
shown in Figure 2 and in Table 2, confirm that this task can also be successfully solved. The patterns of pressure and 
temperature distributions for different blade profiles have common features, but the detailed picture differs. It can be seen, 
specifically, that the blade of the profile [19], which has a larger cross-sectional area compared to the others, is 
characterized by a higher value of relative pressure in the area of the outlet edge and a lower cooling rate. This is also 
illustrated by the data given in Table 2. Since, from the point of view of the strength of the adhesive bond of the coating 
with the substrate, the most dangerous forces are directed normally away from the surface of the blade, the relative 
minimum pressure near the surface is relevant. Its values are presented in Table 2. With account to the calculation 
obtained, it is possible to reasonably approach the design of a special coating for the surface of a blade of a specific 
profile: adjust the materials and thickness of the coating depending on the degree of loading of individual surface areas. 

Note that these calculation results are presented only as an illustration of the possibilities of the proposed profile 
approximation technique. Therefore, the model of gas dynamics used in the work is quite simplified. In a specific task, it is 
necessary to use a more rigorous model with parameters corresponding to the specified type and mode of the GTE operation. 

Thus, the proposed technique for approximating the blade profile of the GTE, based on the use of orthogonal 
polynomials, is a convenient tool that provides automating the development of a geometric model of the blade and 
comparing different types and profiles of blades, solving the corresponding problems of gas dynamics. At the same time, 
for a given blade profile and working conditions of a GTE, it is possible to obtain a distribution of temperatures and forces 
acting on the blade surface, which is required to predict the durability of special coatings [21]. Based on this calculation, 
it is possible to optimize the coating technology by varying its thickness and/or composition, increasing the coefficient of 
durability in the most dangerous areas. 
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Abstract 

Introduction. In recent years, the development of high-voltage power systems has received a boost due to the need for 
infrastructural support for priority development areas. Universal models and algorithms are required to implement 
processes in power components and identify their optimal parameters. However, there are no such solutions. Accordingly, 
there are no ready-made subsystems with control and optimization algorithms adequate to the tasks under consideration. 
The objective of the presented research is to develop an optimization subsystem for the design of control and measurement 
components of power distribution systems. 

Materials and Methods. Methods for constructing automated design systems, optimization, system analysis, 
mathematical modeling, and adaptive control were used. When selecting methods, we proceeded from the fact that the 
components of power distribution systems consisted of a finite number of elements. The synthesis of a power system 
includes tens or hundreds of sequential operations. This was taken into account in the developed models and algorithms. 
Results. The possibilities of managing and monitoring manufacturing processes (MP) for the production of components 
of low-voltage power distribution systems were shown in terms of checking the operability and correct functioning of 
processing equipment. A modular structure was created to allow the integration of CAD output files into the 
manufacturing processes of energy distribution system components. A functional diagram of a subsystem for control and 
monitoring of the manufacturing processes of the production of components of power distribution systems was developed. 
The proposed schematic diagram of production control showed how the data collection subsystem, management system, 
and operating mechanisms were involved in the control of operations. The multi-level optimization module model created 
within the framework of this research sequentially optimized the service intensity of the i-th block, the input flow 
separation coefficients, and the priorities of the original data flows that form the input flow of the i-th block. 

Discussion and Conclusion. The combined application of modeling, system analysis, and optimization methods 
maintains control of the accuracy of the generated power components. The algorithm for controlling electrical loads opens 
up opportunities for creating a mathematical model of a power supply system that combines management, control, and 
monitoring, which ultimately leads to an improvement in the quality of electric power. The solution can be in demand in 
the development of power systems of priority development areas. 


Keywords: improvement of electric power quality, power distribution system, data flow in the optimization module, 
multi-level optimization model 
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AHHOTalna 

Beedenue. B nocieqHue Tob! pasBHTve BBICOKOBOJIBTHBIX IHEpreTHYeCKHX CHCTeCM TOYO HOBBIM MMITyIIbC B CBA3H C 
HeoOxOAHMOcTHIO MHdpacTpyKTypHoro obecredeHuaA TeppuToOpuli onepexaromlero pa3BuTHa. HyxHbI yHHBepcasIbHEle 
MOJeM MW asrOpHTMbI JIA peasH3all“u MpoleccoB B 3HepreTw4ecKHX KOMMOHEHTAaX UM BbIABIICHHA MX ONTHMaJIbHBbIX 
TlapaMeTpos. OqHaKo Take pellieHHA OTCYTCTBYIOT. COOTBETCTBeHHO, HeT FOTOBBIX IOJICHCTeM C asIrOopHTMaMM yiipaBsieHHa 
MW ONTHMM3allMH, ajleKBaTHbIMH paccMaTpHBaeMbIM 3aqjauaM. Ile MpeycTaBleHHoro UccieqoBaHua — pa3pa0oTKa 
TIOACHCTeMBI ONTHMH3alIMH pH MpoeKTHpOBaHHH KOHTPOJIbHO-H3MePHTeIbHBIX KOMIMOHCHT pacripes(evIMTesIbHBIX 
93HepreTuyeckux CHCTeM. 

Mamepuanot u memodot. Victionb3yrorca’ MeTObI MOCTPOCHHA ABTOMATH3HPOBaHHBIX CHCTeM IpoeKTHpoBaHHUa, 
ONTHMU3allHM, CHCTeMHOrO aHasIM3a, MaTeMaTH4eCKOrO MOJeIMpOBaHHA MW alallTuBHOrO yipaBieHua. IIpu BErdope 
MeTOJOB HCXOJMIIM M3 TOPO, YTO KOMIOHEHTHI paciipeeIMTeIbHBIX 3ICKTPHYCCKUX CHCTEM COCTOAT M3 KOHE4HOTO YNCIIAa 
3JIEMCHTOB. CHHTe3 JHEpreTHYeCKOM CHCTeMBI BKIFOUACT ECATKH WIM COTHM MOCeOBAaTeCIbHBIX Ollepalluil. ITO yaTeHO 
B pa3pa0oTaHHBIX MOJeIAX HW asIrOpHTMax. 

Pe3ynemamot ucciedoeanua. TloKa3aHbI BO3MOXKHOCTH YUpaBICHHA HM KOHTPOJIA TeXHOJOrMYecKHX Mpoueccos (TIT) 
IIpOH3BOACTBAa KOMIMOHCHT HU3KOBOJIBTHBIX pacIIpeyeIMTebHBIX 3IHeEPreTHYCCKHX CHCTeM B IlaHe MpOBepKu 
paOotocnocoOHocTH HM KOppeKTHOCTH (yHKIMOHUpOBaHHA TexHOOrMYecKoro ObopynqoBaHHuA. Co3faHa MOsyIbHad 
CTpyKTypa, Mo3BONAFONIad WHTerpHpoBaTb BbIXOTHbIe dain: CAIIP B mpolecchr mpou3BoycTBa. Pa3padoTaHa 
(PyYHKWMOHasIbHad CXeEMa MOJCHCTeMBI yiIpaBJICHHA H KOHTPOIA TEXHOJOMYeCKHX MpOMeCccOB NPOH3BOACTBa KOMIIOHEHT 
pactipeweMTeIbHBIX 39HepreTuYecKux cucTem. IIpenoKeHHad NPHHIUMMWalbHat CX€Ma KOHTPOJIA TpOM3BOJCTBAa 
TOKa3bIBaeT, KAKUM O00pa30M B KOHTpoOle Ollepalluii 3aeHicTBOBaHbI NOACHcTemMa cOopa JaHHBbIXx, CHCTeMa ypaBJIeHHa 
yupaBJIAtOllwe MeXaHW3MbI. Co3qaHHad B paMKaxX JaHHOM paOoTb] MHOTOYpOBHeBad MOJeb MOAYIA ONTHMU3allMH 
TOCII€IOBaTeIbHO ONTHMU3UpyeT MHTCHCHBHOCTS OOcIyKUBAaHNA i-ro O0Ka, KOIPPMIUMCHTHI pa3qeueHHA BXOHOTO 
TOTOKAa H IIPHOPUTeTHI MCXOAHBIX MOTOKOB JaHHBIX, OOpa3yrOIHX BXOAHOM MOTOK i-ro OI0Ka. 

Ooécyocdenue u 3akiiouenue. KomiiekcHoe TpwMeHeHHe MeTOJOB MOseIMpoBaHHaA, CHCTeMHOrO § aHasin3a, 
ONTHMU3alIMM OOeCIeUNBacT KOHTPONb TOUHOCTH POPpMUpPyeMBbIX IHeEpreTHYeCKHX KOMIOHEHT. AJTOpHTM yipaBeHHa 
3ICKTpHYCCKHMU Harpy3KaMM OTKPbIBACT BO3MO%KHOCTH IA CO3qaHWa MaTeMaTH4eCKOM MOsesIM CHCTeMBI 
3HeprocHaOxKeHHA, KOTOpad OOBeAMHACT ypaBsIeHHe, KOHTPOJIb, MOHHTOPHHT, 4YTO B KOHEYHOM c4eTe BeeT K 
yiIyuleHHro KauecTBa 3IeKTposHepruu. PemieHue MoxeT OBIT BOCTpeOOBaHO Ip pa3BHTHU IHepreTw4ecKUX CHCTeM 
TeppHTOPUH onepexkarollero pa3sBuTHA. 


Ksnouesbie cJI0Ba: yTYUmeHHe KayecTBa JJICKTPOSHeprHun, pacupeesIMTeIbHaA SHepreTH4eCcKaA CHCTeCMa, MOTOK 
JWAHHBIX B MOJTyJIe OUTHMH3aUWHH, MHOTOYPOBHEeBad OLTHMV3aWHMOHHAA MOJCJIb 


baar ‘OapHocTH. ABTOpBI Omar OapHbI BopoHexckoMy WMHCTUTYTY BbICOKHX TeXHOJIOrMH 3a TIOMOMb Ip HOATOTOBKE CTaTbH. 


Jia warupopannsa. Kmmenko FO.A., JIbposwy 5.E., [peodpaxencxuii AI]. Tpoexruposanve KOHTpOsIBHO-H3MepHTeJIBHBIX 
KOMHOHEHT paclipeyeuIMTeyIbHBIX  3HepreTHyecKux cHcTemM. Advanced Engineering Research  (Rostov-on-Don). 
2024;24(1):88-97. https://doi.org/10.23947/2687-1653-2024-24- 1-88-97 


Introduction. The creation and development of power distribution systems with the required parameters is an 
important science and technology task. Active theoretical and applied developments are underway in the high-voltage 
energy area. In recent years, priority development areas have been formed and developing in Russia. These are production 
clusters that need to be provided with high-quality infrastructure, including energy infrastructure. In most cases, we are 
talking about small settlements that initially lag behind in economic, social and infrastructural terms. The solution of 
power problems under such conditions should be facilitated by the use of universal models and algorithms that will, in 
particular, identify the optimal parameters of the components of power processes. Such solutions are not presented in the 
literature. Accordingly, there are no management and optimization subsystems developed on their basis. 

The high cost of modern power equipment should be noted. Its effectiveness assumes that the setting provides optimal 
parameters for the operation of the components. When designing monitoring-and-measuring components of power 
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distribution systems, it is important to specify control and tolerance points regarding the quality parameters of the 
corresponding process steps [1]. During production, processing sequence should be determined. 

Let us consider the last situation, 1.e., production, in greater detail. A large number of process steps can create 
difficulties associated with control and management [2]. Therefore, when releasing components of power equipment, it 
is important to specify the requirements for the subsystem of control and optimization of manufacturing processes (MP). 
Such a subsystem is used in the manufacture of prototypes. Another recommended approach is the examination of power 
modules. It provides determining adequate output parameters. It is also important to take into account the influence of 
external actions [3]. 

The selection of algorithms should be based on the methods of mathematical statistics. This makes it possible to 
control and optimize the quality parameters of processes and use appropriate management procedures [4]. Versatility 
should be noted as an advantage of the algorithms. By the provision of it, within the framework of a systematic approach, 
it is possible to plan and implement research in this direction [5]. It is known that the start of manufacturing new power 
equipment is not always provided with the required statistical data. In this case, it makes sense to apply adaptive 
management methods. Self-optimization procedures are convenient to use when there are changes in the requirements for 
MP, as well as when external conditions vary [6]. The objective of the presented research is to create an optimization 
subsystem in the design of monitoring-and-measuring components of power distribution systems. 

Materials and Methods. Methods of creating automated systems for design, optimization, system analysis, 
mathematical modeling, and adaptive management are used. When selecting methods, it was taken into account that the 
components of power distribution systems consist of a finite number of elements. In more detail, the task can be presented 
as follows. It is planned to produce equipment components for power distribution systems. It is required to create a 
subsystem that will manage and optimize engineering MP. The analysis provides selecting and using the adaptive control 
method of such a subsystem [7]. 

The synthesis of the entire power system involves dozens or hundreds of sequential operations, and this is taken into 
account in modeling and algorithmization. 

In the production of components of power complexes, it is proposed to use simulation modeling based on the MP 
optimization module [8]. At the same time, the quality of the components is considered in relation to the parameters of 
the created power distribution system. Uncontrolled MP parameters can also be taken into account. 

Figure 1 shows the integration of computer-aided design (CAD) output files into the production of components of 
power distribution systems in the form of modules. 


Operator Interface CAD 


Module containing standards 


; Synthesis of power components 


: ombining power components 
Mocule raed ; Implementation of modeling processes 
to equipment selection : 


Module for optimization 


Generating output files for production 


Database 


of manufacturing processes 


Fig. 1. Modular structure of integration of CAD output files and production of components 
of power distribution systems 
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The control module provides transmitting the process data to an automated control system (ACS) [9]. During the MP, 
operators receive data on the reliability of the power components being created. At the same time, the results of the 
physicotechnical expertise are taken into account, and the optimization module is used. The impact of various factors on 
the efficiency of power equipment is monitored. The data on previously manufactured components are analyzed. In further 
developments, the proposed module will make it possible to save and take into account several types of hardware settings. 
They can be used in production to determine tolerance ranges for each process step [10]. 

The ACS should include four modules. 

1. Subsystem of interoperative control support. 

2. Transport operations management module. 

3. Subsystem for managing MP modes. 

4. Subsystem of process steps management. 

Such ACS translates information from designers to production. In practice, various quality management systems 
are used for the components of power distribution systems being created, including those with production control 
subsystems [11]. 

In this paper, an automated subsystem for control and optimization of monitoring-and-measuring components of 
power complexes is proposed. It provides different types of impacts on MP for any stage of production. The previously 
obtained data allow the subsystem to reduce the number of defective elements, i.e., to improve the quality of products. 

Figure 2 illustrates the structure of this subsystem. It was formed on the basis of the requirements for the efficiency 
of production processes. In addition, different types of components being formed were taken into account. Previously 
used values of MP parameters can be applied for analysis. The operator has the opportunity to change them. 


Programs Programs Programs 
for dialog interaction for interface to process messages 


| {4 


Input and output 
module 


Main program 


Automated task setting module 


Target ‘ : Libraries 
Function Compiler es eee acest oa of basic algorithms and models 
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Adaptive for non-standard conditions for passing parameters 


control algorithms 


Module including application models 


Block for adjusting 


mathematical models Block for production 


of power components 


Database 


Precision Setting Power Component Design Data 


Transfer procedures Block Processing Unit 


Block for processing results 
of manufacturing processes 


Control Module 
Optimization Module 


Fig. 2. Functional diagram of subsystem for management and control of MP 
components of power distribution systems 
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The main program is considered as the basis of the automated task module, and for the processes of setting up and 
transferring MP parameters. An analysis is carried out to identify the need to adjust equipment settings. A special role is 
played by programs for processing messages and supporting dialog interaction. They make it possible to adjust the MP 
manually. The operator has the ability to download current settings from the database. They are set manually or through 
an input and output module. 

The optimization module generates a list of top-level management operations for a specific MP. Then, the 
operations are analyzed, the optimization accuracy and production parameters are set. According to the information 
received, the application models are transferred to the module. At the same time, analytical and statistical models 
of basic MP operations are used. For each such operation, the input and output parameters are determined in the 
corresponding part of the system. 

The control module develops and refines the mathematical model of production processes [12]. It is then used to 
determine the tolerance for quality parameters within the framework of individual process steps. The model is developed 
in the module of applied models. The automated task module is designed for procedures to maximize the target function. 

The received settings for equipment, management and process parameters are recorded in the database. Then, 
during the production of power components, the process is analyzed, and the current settings of the equipment 
elements are adjusted. 

The study of MP management and control is required to verify the operability of processing equipment, the correctness 
of its operation and the organization of production. 

Research Results. The monitoring-and-measuring components of power distribution systems consist of a finite 
number of elements. The component production involves dozens or hundreds of sequential operations and relies on 
solving the problem of forming components with maximum accuracy. At the same time, uncontrolled parameters 
should be taken into account for any step. For this purpose, correlation analysis, adaptive management, and 
optimization are used. 

Particular models are combined into a general model for creating power components, which can be adjusted (e.g., 
with account of experimental testing of equipment or requirements for specific components). Some algorithms are created 
on the basis of the results of monitoring the formation of components and, as a result, provide improving the quality with 
changes in operating modes. The main indicators of the efficiency and stability of production can be considered the values 
of the output indicator of valid power components. 

The following are three key features of managing the processes of component formation. 

1. It is important to maintain the specifications at the required levels for input and output components for 
different batches. 

2. Algorithms affect the shutdown of each processing step. 

3. Constant wear of equipment requires regular adjustment of process parameters. 

We describe the model based on the MP of creating power components [13]. 

Suppose that the (i-1)-th operation of MP is being considered: 

uj = F(ui,vi). (1) 


Here, u; — quality parameter of the power components of the current operation; v; — production option. 

It is important to take into account that in practice not the quality parameters of the power components are critical, but 
the design parameters that depend on them (e.g., highspeed response). Then: 

gi =F (uia,k). (2) 

Here, gi — controlled parameters of the current operation; k; — design parameters. 

The production of power components is described as a trajectory with a change of state. According to the final state, 
the adjustment of the MP is carried out. Within the framework of trajectory management, the corresponding tasks are 
solved. The optimal solution is selected from a variety of solutions [14]. At the same time, the specified characteristics of 
the generated power components are taken into account. In this case, the control actions should provide the best match of 
the required and output characteristics of the components being formed K;’. Here, n — number of MP operations. 

The basic industrial technology determines the initial data, the processing sequence, and sets limits on control variables. 

When selecting the number of MP operations (n—1), the target function is as follows: 

F, = Rick yelude |. (3) 


Here, c)...c;, — selectable process control variables. 
We introduce parameter @;, which characterizes the corresponding goals associated with the i-th operation. For (7-1) 
operations in MP, the following expression is valid: 


fn (Kn ) _ min(@n-1 [Kn seas | + Sn-2 (Kia )) =min Ch-1 (KnsOna ). (4) 


Cn-1 
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If there is one operation in the MP, then: 


fi(K2) =min(K2,h). (5) 
cl 
Significantly, 
fo (K 1) =0. (6) 
Under such conditions for MP, the target function is defined as follows: 
Ci (K2,c1) = @(K2,c1). (7) 


The schematic diagram in Figure 3 shows how the power components are generated. The MP data acquisition 
subsystem (sensor system and instrumentation) transmits the values of input and output parameters for each processing 
step 7 to the control system. This happens before and after each step n (1 = 0,1,2, ..., m). The values are recorded in the 
data collection subsystem. 


i Control z = Data 


Control system 


mechanisms acquisition subsystem 


Fig. 3. Schematic diagram of control of production of power components 


In total, m steps are considered. The work ends after performing m+1 steps. 

The use of the optimization model makes it possible to form technologies for converting the input data flow (xsx). The 
main channels for processing input data are: 

— process department (PD); 

— production department (PRD), n = 1; 

— N of structural subdivisions (SS"). 

Figure 4 shows the block diagram of the transformation of the input data flow. 
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Fig. 4. Block diagram of data flow distribution in the optimization module 
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Here, Ygx — input data flow to be processed in PD; Xgrixy,X grxo»Xguxs — data flows at the output of PD; 
XB — input data flow for processing in PRD; X puss X prixo — data flow at the output of PRD; x§% — input data 
flow for processing in the n-th (n = 1, N) SS"; XGhix1,X Suixo — data flows at the output of SS”. At the same time, the 
flow at the output of the optimization module Xgpx =X ae 

The scheme efficiency is determined by three factors. 

1. For all types of input data flows for any block in the system, the intensity and maintenance mechanisms are considered. 

2. Input flow for each block is divided into several output data flows. 

3. Depending on the level and priority, it is needed to observe the data queue in each block. 

It is possible to analyze these processes using a multilevel optimization model (Fig. 5). 


AsxiRiy 


Optimization of service intensity of the i-th block 


Optimization of separation coefficients of the i-th input flow ZL 


eye 
Sid BbIXli 


Optimization of priorities of initial data flows 
forming input flow of the i-th block 


Selection of power component resources 
and formation of data flow patterns 


Fig. 5. Multilevel optimization model of the optimization module 


Here, Asxi, i=1,/ — intensity of the input data flow associated with Xpgx; i-th block; i =1,7 — block numbers for 
power components; R;; — guaranteed level for the j-th (_j =1,/ ) resource in the i-th block; j = 1,7 — numbers related 
to resource provision; p;, i = 1, — intensity of data management associated with the i-th block; 14; — optimal intensity 


value; vi; — data separation coefficient of the i-th flow to the Ist output; / =1,Z — designation of numbers in output 


flows; vj; — optimal value in the separation coefficient; & =1,K — values of the numbers of source data flows for the 


input flow in the i-th block; k* =1,K~ — numbers of the source flows that will show the optimal priorities; 
tea) — time of interruption of data k; management by data kj with mixed priorities; ¢;; ;;, — optimal time for interruption 


of data k; management by data ki. 


Klimenko YuA, et al. Design of Instrumentation and Control Components of Power Distribution Systems 


The optimal pattern of priorities for the lower level is determined. The optimality problem is solved according to the 
Bellman principle, i.e., the optimality principle is proved to the contrary. It is accepted that part of the process is not 
optimal according to the quality criterion. The criteria of the initial and optimal processes are compared. Based on this 
comparison, a conclusion is drawn about the proof of the optimality principle [15]. The selection is made during the 
processing of different types of data according to optimal, absolute, and mixed priorities. The task scheduler can use 
various algorithms to make decisions about the order of tasks. In some cases, optimal priorities are useful, in others, 
absolute ones. Tasks with a higher priority are completed earlier, and, accordingly, the response to them is faster. When 
combining (mixing) priorities for individual applications, non-priority service is possible. 

The optimal mechanism for dividing into several output flows at the middle level is selected by the input data flow 
for each block. 

At the upper level, the following are selected: 

— input flow balanced in intensity; 

— resource provision of the module for optimizing the intensity of management for each block. 

In this case, a gradient approach can be applied. 

For energy components, a three-level parametric optimization is used in the selection of resources and the formation 
of data flow patterns for each block. 

Discussion and Conclusion. A subsystem for the management and control of MP for the production of monitoring- 
and-measuring components of power distribution systems is created. It provides data collection and defect analysis 
depending on the hardware settings. In addition, the solution can obtain requirements for equipment settings to achieve a 
given level of component quality. 

Optimization of the design is required, specifically, when modeling the control of electric loads to improve the quality 
parameters of electrical power in 0.4 kV power distribution networks. The algorithm of the electric load control process 
was studied in an adaptive system of control and management of the power quality [15]. Its operability has been 
confirmed; therefore, the solution can be used in the development of equipment for 0.4 kV distribution networks. In 
addition, this algorithm can be applied to develop a mathematical model of the power supply system, which is a set of 
functions: monitoring, management, control. The use of modeling, system analysis and optimization methods, maintains 
control of the accuracy of the generated power components. Adequate implementation of this approach can improve the 
quality of electric power. 

The results of the presented research are practically applicable, in particular, for solving problems related to the power 
supply of the priority development areas. 
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Abstract 


Introduction. To date, numerous methods of nucleic acid amplification have been proposed, and each method has a 
number of advantages and disadvantages. One of the most popular methods is Loop-mediated isothermal AMPlification 
(LAMP). Unlike thermocyclic reactions, such as PCR (polymerase chain reaction), which require three temperature 
changes and expensive equipment, in LAMP, the entire reaction takes place at one and the same temperature and at the 
maximum rate possible. An important component of LAMP is primers (usually 20—25 nucleotides), which need to be 
matched to a specific part of the nucleotide sequence. It is known that DNA sequence contains four nucleotides: 
A — adenine and T — thymine, G — guanine and C — cytosine. There is a huge variety of permutations of these 
nucleotides, and it is practically impossible to analyze such a large amount of data manually. Therefore, there is a need 
to use modern computer technologies. More than 150 computer programs have been proposed for the design of PCR 
primers, while for LAMP primers there are less than 10 of them, and each of them has a number of drawbacks, e.g., in 
terms of the length of the analyzed site. Therefore, this work is aimed at developing a new domestic computer program 
for the design of specific primers for LAMP. 

Materials and Methods. The primer search algorithm was based on a linear search for a substring in a string, taking into 
account the criteria of primer selection for LAMP. The program complex of LAMP-primer design was implemented in 
Python programming language. The bioPython library was used to work with various DNA and RNA, and the Qt framework 
was used to develop the interface. 

Results. A modification of the direct sampling method using a stencil approach was proposed, taking into account the 
GC composition and annealing temperature of primers depending on their structure. A software package with a friendly 
interface was developed. It took into account the design criteria of primers: certificates of registration of computer programs 
(LAMPrimers iQ No. 2022617417 dated April 20, 2022, LAMPrimers iQ_ loop No. 2023662840 dated June 14, 2023) were 
received. The program is in the public domain at https://github.com/Restily/LAMPrimers-iQ 

Discussion and Conclusion. The developed software packages can be used for research and analysis in molecular biology 
and genetics, to create diagnostic test systems that provide high sensitivity and reliability of detection of specific DNA 
and RNA. The software packages can be used in research institutes and laboratories engaged in the amplification of 
nucleic acids. The results of evaluating the selected sets of primers for the LAMP reaction were tested, and the 
effectiveness of working sets using the LAMPrimers iQ program was experimentally proven by the example of the 
detection of genetic material of the SARS-CoV-2 coronavirus. 
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AHHOTalna 

Beedenue. Ha ceroquamuuii eHub cyujecTByeT MHO%KeCTBO ciHocoboB aMIIMpuKalMH HYKJICHHOBBIX KHCJIOT Uy 
Ka@KOFO ciocoba ecTb pA AOCTOMHCTB MH HEAOCTAaTKOB. OHM V3 HaHOolIee MOMYIAPHBIX CIOCOOOB ABIIACTCA MeTICBAA 
w30TepMuyeckad amimMdpuKaya (Loop-mediated isothermal AMPlification, LAMP). B otmmune oT TepMOWMKIM4eCKHX 
peak, Taxux Kak IIL[P (nomuMepa3Had WelHad peakHA), WJId KOTOpbIX TpeOyeTca CMeHa Tpex TeMIepaTypHBIX 
pe*KUMOB H AoporocTosmjee obopyzoBaHne, B LAMP sca peakiua NpOXOANT Mp OAHOM Temiepatype uC MaKCHMaJIbHOM 
Ha J{aHHbIM] MOMeCHT CKOpocTbi0. BaxkHbIM KOMIOHeHTOM mpoBeyzeHuat LAMP-ammimduKkalMu aABJIAIOTCA MpaliMepbl 
(oO6prHO 20-25 HyKMeoTHOB), KOTOpble HeoOxoqHMo WosOupaTb K oOlpeyeieHHOMy y4acTKy HYKIeOTHIHOM 
TocueqoBatesbHocTH. V3BecTHo, 4uTo nMoceqoBaTenbHocTb JIHK coyepxxuT uerbipe HyKeoTuya: A — aeHuH U 
T — tTumun, T — ryanuu u I] — yutTo3un. BapuaHTos mepecTaHOBOK 3THX HYKJICOTHOB OFTPOMHOE MHO2XKeCTBO, H 
TipOaHaIM3HpOBaTb BPYYHY!IO Takoe OObINOe KOMMYCCTBO Ja@HHbIX MIpaKTHYCCKH HEBO3MO%KHO, MOITOMY BOSHHKaeT 
HeOOXOAMMOCTb B MCIOJb30BaHHH COBPCMCHHBIX KOMIBIOTepHbIX TexHONorHi. Ja gAu3aiina IIL[P-npaiimepos 
Tipefo%xKeHO Oosee 150 KoMMbIOTepHEIX IporpaMM, B TO BpeMsA Kak JI LAMP-mpaiimepos ux Mexee 10, u Kaxxyaa 13 
HHX HMECT pA HEAOCTATKOB, HallpiMep, M0 JJIMHe aHasM3upyemoro yyactka. TlosTomy Wesb1o JaHHOM paoorTsl ABJIAeCTCA 
pa3pa0oTka HOBO OTeYeCTBEHHOM KOMIIBIOTepHON IporpaMMbI AH3aiHa CHeLMPHYHBIX UpaliMepoB HMeHHO Aa LAMP. 
Mamepuanoi u memoooi. B ocHose aulropuTMa MOoMcKa WpaliMeposB JIexKUT JIMHeMHBIM NOMCK MOACTPOKU B CTpOKe C 
yueToM KpuTepHeB 1oA00pa mpalimepos Aa LAMP. IIporpammuniit kommiekc 2u3aitHa LAMP-npatimepos pa3paboTaH 
Ha A3bIKe IporpaMMupoBanns Python. Jia padotsr c pasmmunpimu JIHK u PHK ucnomp3o0Banacb OudsMoteKa bioPython, 
a JIA paspaOoTKu HHTepPelica — PpeiimBopkK Qt. 

Pe3yrvmamot ucciedoeanua. Upeqnoxena MOANPUKalNA MeTOa IIpaAMoro HepeOopa c HcIIOIb30BaHHeM TpadapeTHOrO 
ToAxoa, YAHTHIBarolero GC-cocTaB HW TeMIlepaTypy OTMKHTa lpaliMepoB B 3aBHCHMOCTH OT UX CTpyKTypbl. PaspabotaH 
KOMIIJIeKC TIporpaMM Cc Apy?KeIOOHbIM MHTepdelicoM, YYHTbIBAaIOWHH KpHTepHu AM3aliHa MpaiiMepoB: Mosy4eHbI 
CBUeTeIbCTBa O perucTpalun mporpamM jwia IBM (LAMPrimers iQ Ne 2022617417 or 20 anpena 2022 roga, 
LAMPrimers iQ loop Ne 2023662840 ot 14 urona 2023 rona). [porpamma ecTb B OTKPbITOM WOcTyHe oO aspecy: 
https://github.com/Restily/LAMPrimers-iQ 

O6cystcoenue u 3akio4uenue. Pa3zpa0oTaHHble IporpaMMHBle KOMIICKCbI MOTYT HCHOMb30BaTBCA JIA HCCIeOBaHHH U 
aHaiv3a B OOaCTH MOJIeKyIApHOM OMONOrMH HM YTeHeETHKH, Id CO3faHHA MarHOCTHYeCKHX TeCT-CHCTeM, 
oOOecrIeYMBaIOWIHX BbICOKYIO UYBCTBHTCJIbBHOCTb HM JOCTOBepHOCTs OOHapyxKeHuA creuMpuyeckux JIHK u PHK. 
IIporpaMMHbIe KOMIWICKCbI MOTyT IIPHMeHATBCA B HayYHO-HCCIeAOBAaTeIbCKUX MHCTHTyTax uM sabopaTopusAx, 
3aHMMalOWHXCcd aMIVIMpuKalWei HYKIICHHOBBIX KHCIOT. PesyIbTaTbl OMCHKH WOAOOpaHHbIX HaOopoB UpaliMepos JVIA 
peaxuuu LAMP anpoOupopaHsl, u 3=deKTHBHOCTb padounx HadopoB Cc MOMOLIbIO MporpamMmpr LAMPrimers iQ 
OKa3aHa IKCICPHMeCHTAIbHO Ha IpHMepe OOHapy2xKeHHA TeHeTHUeCKOrO MaTepHasa KopoHaBupyca SARS-CoV-2. 


Kono4uesbie cs10Ba: 2173aiH Tpalimepos, leTeBbIe mpaliMepsl, python, C++, nerepad H30TepMHYeckKad amrumMuKarua, LAMP 
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BaarogapnHoctu. Astop OmarofapuT cBoero HaydHOrO PyKOBOAMTeIA, HOKTOpa (:bu3MKO-MaTeMaTHYeCKHX HayK, 
mpodeccopa V.M. Tyoaiinzynmua u Komer v3 Muctutyta OvoxumMuu u reHeTuKH YOU] PAH 3a koncybTayMn B xoye 
BBINOJIHEHHA paooTHl, a TAaKKe 3a IKCIECPHMCHTAJIBHYIO allpOOalHi0 Pe3YJIbTATOB B JIaAOOpaTOPHbIX YCIOBUAX COBMECCTHO 
C aBTOPOM CTATBH. 


@unancuposanne. Padota BEmlowHeHa pH duHaHcoBoH MoATepxKKe PDODY B pamKax rpauta Ne 20-37-90091 
«Pa3paOoTKa IporpaMMbl Ju3aiina Mpalimepos Aa Loop AMPlification — neTmeBoi u30TepMM4ecKOn aMIIMpukaluu 
Ha OCHOBe TEXHOJIOrMM MalIMHHOro obyyeHH»», https://www.rfbr.ru/rffi/ru/contest/o 2088005 


Jia uuTupopanusa. Axmet3saHosa JI.Y. KommbrotepHaa uporpamMa nogOopa upaiimepos aia LAMP-amnmoukauun. 
Advanced Engineering Research (Rostov-on-Don). 2024;24(1):98-108. https://doi.org/10.23947/2687-1653-2024-24-1-98-108 


Introduction. Nucleic acid amplification is a valuable molecular tool not only in fundamental research, but also in 
applied areas, such as the diagnosis of infectious diseases, hereditary pathologies, establishing kinship, etc. Currently, 
amplification methods are intensively developing, and their application areas are expanding. The most popular and most 
frequently used amplification method is polymerase chain reaction (PCR) [1]. PCR is a reaction that takes place under 
three different temperature conditions: denaturation (95°C), annealing of primers (from 50° to 60°C), elongation (72°C). 
To quickly change these modes, a special device is needed — a DNA thermocycler [2]. At this, temperature variation in 
the amplifier does not occur instantly, but starts only when the desired temperature is reached, and this causes artificial 
containment of the reaction. As a rule, the duration of PCR is 1—1.5 hours. 

The second most popular amplification method is loop-mediated isothermal amplification (LAMP) [3]. A water 
bath or thermostat is sufficient for LAMP, since the reaction takes place at one and the same temperature, and the 
first results can be seen in 15 minutes. 

For both LAMP and any other type of amplification, the key component is primers, which are short sequences 
of nucleic acid. They serve as a starting point for increasing copies of a specific region of DNA. It is the primers 
that determine which DNA sequence will be copied. 

The main difference of LAMP implementation is the number of primers. For a conventional LAMP, at least four 
primers are required (two external, two internal), while for a conventional PCR, two are sufficient (direct, reverse). 

To increase the specificity and accuracy of the reaction, it is important to select the right primers. For the 
automatic selection of primers for PCR, more than 150 different computer programs have been developed that 
provide the selection of primers for any modifications of this reaction [4]. However, there are very few such 
programs for LAMP, no more than ten, and only two of them are available online. These programs also have a 
number of disadvantages, such as restrictions on the length of the analyzed sequence; they do not exclude the 
possibility of formation of homo- and heterodimers of primers, repeats of nucleotides in one primer. And none of 
the programs take into account the close arrangement of primers in one set, which, in turn, reduces the quality of 
primers and the accuracy of reaction results [5]. 

Thus, an urgent task is to develop a new computer program to select (model) high-quality sets of primers for 
LAMP with tougher conditions for choosing primers for nucleotide sequences of any length. 

Materials and Methods. The authors of [3] proposed using two external, F3 (Forward), B3 (Backward), and 
two internal primers, FIP (Forward Inner Primer), BIP (Backward Inner Primer). It was assumed that the internal 
primers had double length (FIP: Flc/F2, BIP: Blc/B2) and were annealed at four regions of the nucleotide sequence. 
Schematically, the location of LAMP primers can be seen in Figure 1. 

External primers are only needed at the initial stage. They are designed to limit the analyzed region of the 
nucleotide sequence and form a single-stranded structure of this region. A pair of internal primers, Fle and Blc, 
start their work already at the second stage, as they are annealed after the formation of new DNA chains. 

Later, the same authors proposed a modified method offering the use of not four, but six primers, annealed 
already at eight regions of the target nucleotide sequence [7]. It was proposed to add two more loop primers (Loop 
B, Loop F), which should react at the third stage after the formation of a dumbbell-like DNA structure and anneal 
between regions F1/F2 and B1/B2, respectively. The use of additional primers implies an increase in sensitivity and 
reliability of the reaction. 

Any amplification reaction has its own sensitivity threshold, and the spread of this indicator is very large due to 
the fact that numerous factors affect the course of both PCR and LAMP. Some papers have noted that LAMP is 
significantly more sensitive than PCR. The authors of [8], e.g., claim that LAMP is 10 times more sensitive than 
PCR. The authors of other studies have found that LAMP is 100 times more sensitive than PCR [9], and in some 
investigations, this indicator reaches 1,000 times [10]. 


Akhmetzianova LU. Computer Program for Primer Design for Loop-Mediated Isothermal Amplification (LAMP) 


Analyzed nucleic acid 


B FP Fl y Ble B2e B3c 


F3c F2c Fle Bl B2 B3 
F3 ____FIP 
Bl 
Fl B2 
F2 
Fic 
+ 
Flic B2 
F2 
Fl 


Fig. 1. Schematic arrangement of the annealing regions of external and internal primers for LAMP [6] 


In addition to sensitivity, any amplification reaction has another equally important indicator — its specificity. 
And here, questions have recently begun to arise about the LAMP method [11], including due to the emergence of 
the so-called primer homo- and heterodimers, which are more difficult to exclude in this reaction than in PCR, 
because of the larger number of primers used and their increased length [12]. 

For successful amplification, it is necessary to select the right primers. When using the LAMP method, the main 
difficulty is in modeling primers with account for all recommended conditions, namely: 

1) length of the primer (18-35 nucleotides for external primers, 30-55 nucleotides for internal ones); 

2) content of guanine (G) and cytosine (C) (GC composition ranging from 40 to 60%); 

3) optimal primer annealing temperature (55—65°C); 

4) close arrangement of primers in one set: the average size of the amplicon (120-220 bp); 

5) exclusion of the formation of dimers of primers; 

6) elimination of nucleotide repeats in one primer (no more than three). 

Table 1 provides brief characteristics and features of the most popular LAMP primer design programs [5]. 
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Table 1 
Brief characteristics of popular computer programs for LAMP primer design 
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The design of primers for LAMP is a very difficult task and requires the development of a special computer 
program with proper functionality, taking into account all recommended conditions, and with the possibility of an 
extended selection of primers and a user-friendly interface. 

The computer program for primer design is developed in the Python programming language. This language has 
the bioPython library, which allows working with nucleotide sequences, as well as the Qt framework for interface 
development. 

Research Results. With account for the structural features of nucleotide sequences and the criteria for selecting 
LAMP primers, a modification of the direct sampling method using a stencil approach has been proposed. It 
considers the GC composition, the annealing temperature of the primers, and reducing the complexity of the 
sampling. 

As is known, the GC composition of primers should be in the range from 40 to 60%. This is one of the important 
criteria for selecting LAMP primers, which depends on the sequence being analyzed, the length of the primer, and 
it partially affects the annealing temperature (7m, °C). 

In this paper, the following formula is used to calculate the annealing temperature of primers: 


Tm =81.5 +16.6-(logio [ Na’ ]) +0.41-(%G+%C)—548/L, (1) 


where [Na*] — molar concentration of sodium ions; (%G+%C) — GC composition in the analyzed sequence, 


expressed as a percentage; L — length of the primer. It was based on a well-known dependence: 


Tm =81.5+ 16.6-(logio [ Na” ])+ 0.41-(%G+%C)-600/length. (2)! 


Formula (1) was selected empirically, the calculated values were compared to the values obtained through the 
convenient OligoAnalyzer? utility, which provided a high-quality selection of primers for all types of nucleotide 
acid amplification. Through determining the length of the primers, GC composition and annealing temperature, all 
possible primers can be found in a nucleotide sequence of any length that will meet the specified criteria. 


' Oligo Calc: Oligonucleotide Properties Calculator. URL: http://biotools.nubic.northwestern.edu/OligoCale.html (accessed: 10.12.2023). 
? OligoAnalyzer™ Tool. URL: https://eu.idtdna.com/pages/tools/oligoanalyzer?returnurl=%2Fcalc%2Fanalyzer (accessed: 10.12.2023). 
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If we imagine the primer as a substring, and the analyzed nucleotide sequence as a longer string, this task can 
be represented as a search of all possible options (direct search), but complicating it through calculating the GC 
composition and annealing temperature of the primers. 

Figure 2 shows a complete block diagram of the direct search algorithm, taking into account: 

— length of the primers (bp); 

— GC-composition, %; 

— annealing temperature of the primers, Tm, °C; 

— homodimers on both DNA strands. 

The total complexity of the modified algorithm in the worst case is O(m-n), where n — length of the primer, 
m — length of the nucleotide sequence. It should be understood that the running time of the algorithm directly 
depends on how often the nucleotide fragments that meet the requirements are found. 

Table 2 shows the search data for all possible primers in nucleotide sequences of different structures over time. 

Further, it is necessary to form sets from all the primers found, taking into account the close distance between 
the primers, the heterodimeric properties in one set, as well as the minimum temperature difference between the 
annealing of the primers. 


The design scheme for forming primers into the LAMP sets is shown in Figure 3. 


Primer search cycle 
Analysis of the first n-1 nucleotides 
GC, flag calculation 


Within search range Returning an array of primers 


Calculation of flags and GC 


for nucleotide of current index 
z Calculation fal 
of parameters 
ii Testing for homodimers _ 
no ; yes 
No homodimers? 
Adding a primer 
Testing for homodimers 
from complementary strand 


es 
: 
Adding a primer 


Calculation of GC and flags 


for the first and last nucleotide 


Fig. 2. Complete flowchart of the algorithm 
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Table 2 


Time spent searching for valid primers in nucleotide 
sequences of different lengths 


Genome name Nucleotide sequence size, bp Search for primers, s 


SARS-COV2 
Escherichia virus T4 168,903 


Mycoplasma 580,076 


Helicobacter pylori 1,624,458 18.11 
Escherichia coli 4,641,652 71.68 (1.2 min) 


Caenorhabditis 100,286,401 1,082.53 (18 min) 


yes 
Loading from a file? 


File download 


Entering a sequence 


Chromosome selection 


Configuration setting 


Search for primers, calculation 
of length, Tm, GC composition 


Sorting primers into sets 


Displaying sets on the 
screen 
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Fig. 3. Block diagram of the formation of primer sets for LAMP 
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At the input, the program reads the nucleotide sequence, either the desired file is loaded, or a fragment is inserted 
through the clipboard. Next, configurations are set, such as primer length, GC composition, primer annealing temperature, 
and temperature difference in one set. Then, all possible parameters satisfying the set configurations are searched, the 
selected primers are sorted into sets and displayed to the user. 

The operation of the computer program: 

1. Uploading a file (simple text format, FASTA format, GenBank), or a fragment of a sequence via the clipboard. 

2. Search for all possible primers: primers are combined according to the following criteria: 

— length of the analyzed region; 

— distances between primers (F3/F2 — 1-10 nucleotides, F2/F 1c — 10—25 nucleotides, Flc/B1c — 0-30 nucleotides); 

— temperature difference of primer annealing (<3); 

— heterodimers. 

If all of the above conditions are met, the set is considered to be working. 

3. Displaying simulated primer sets to the user's screen and/or saving to a file. 

The LAMP program is registered in the Register of Computer Programs under the name LAMPrimers iQ, No. 
2022617417 on April 20, 2022, and LAMPrimers iQ-loop, No. 2023662840 on June 14, 2023. The program code is 
publicly available’. 

The developed software product has a user-friendly and intuitive interface, which can be used directly by end 
users — experimenters engaged in LAMP. 

Discussion and Conclusion. The number of primer sets issued depends on the specified search parameters. The 
stricter the parameters, the fewer sets will be found. In case of hard restrictions, the program may not output a single set. 
The number of primer sets with different selection parameters for the genome of the bacteriophage lambda [13], whose 
length is ~48,500 nucleotides, is shown in Table 3. 

For relatively short nucleotide sequences (up to 2,000 nucleotides), the selection of primers takes less than a second. 
With increasing sequence length, the duration of the primer search grows exponentially. 

Table 3 
Number of the primer sets for the bacteriophage lambda genome, 
depending on the set selection parameters 


Maximum length of amplified region, bp 
GC, % ATm, °C 
300 230 160 

5 213 119 3 
40-60 

2 198 184 3 

5 132 80 0 
45-55 

2 116 64 0 

5 195 185 4 
50-60 

2 181 164 4 

5 160 147 4 
55-65 

2 134 105 0 


Figure 4 shows the effect of the length of the nucleotide sequence on the duration of the selection of primer sets. The 
data for the nucleotide sequence of the bacteriophage lambda on a laptop with the parameters are as follows: Intel(R) 
Core(TM) i7-10750H CPU, 2.60GHz, 6 cores. 16 GB RAM. The parameters of the soft selection of primers are indicated 
(40-60 % GC, ATm = 5, length of the analyzed region is up to 300 bp). 


3 LAMPrimers-IQ. URL: https://github.com/Restily/LAMPrimers-iQ/blob/main/lamp/start_lamp.py (accessed: 10.12.2023). 
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3,000 


2,000 


Set search duration, s 


1,000 


1 2 3 10 20 40 


Nucleotide sequence length x 10°, n 
Fig. 4. Effect of length of nucleotide sequence on duration of selection of primer sets 


It should be noted that the duration of the search for primers depends on the power of the computer. 

To compare and determine the quality of the simulated primer sets, a number of field experiments were conducted to 
detect the RNA of the SARS-CoV-2 coronavirus, whose length was ~ 30,000 nucleotides. To do this, sets of LAMP 
primers were selected for the same region of the nucleotide sequence of the coronavirus using the LAMPrimers iQ 
program and two popular and accessible online utilities from New England Biolabs (NEB LAMP)’ and PrimerExplorer?. 
The designations L, N and P correspond to the sets of primers LAMPrimers iQ, NEB LAMP Primer Design and 
PrimerExplorer; “+”? — samples contained RNA of SARS-CoV-2 coronavirus, “—” — samples did not contain nucleic 
acids. Figure 5 shows the curves of this experiment. 
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Fig. 5. Graph of the comparative experiment (the author’s figure): 
L — primers obtained using LAMPrimers Iq, N— NEB LAMP Primer Design, 
P — PrimerExplorer, samples contained RNA of SARS-CoV-2 coronavirus, 


“oD 


— control samples without matrix (did not contain nucleic acids) 
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4 NEB LAMP Primer Design Tool. URL: https://lamp.neb.com/#!/ (accessed: 24.11.2023). 
5 LAMP primer designing software Primer Explorer. URL: http://primerexplorer.jp/e (accessed: 25.11.2023). 
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The primers obtained through PrimerExplorer, showed the latest rise in amplification curves (P+) compared to 
NEB LAMP (N+) and LAMPrimers iQ (L+). The primers obtained through LAMPrimers iQ, provided a later rise in the 
amplification curves (L+) compared to (N+). However, samples that did not contain virus RNA (P-), showed a later rise 
compared to the set (N—), while (L—) showed no rises even after 50 minutes, thereby providing the highest reliability of 
viral RNA detection. 

The performed experiments showed a higher accuracy and specificity of the primer sets selected through 
LAMPrimers iQ computer program, caused by a decrease in the reaction rate with negative control samples. The 
amplification curves had later rises, or did not have them at all, even after 50 minutes of reaction time. 
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Abstract 

Introduction. Multicriteria optimization, taking into account contradicting criteria, is used to improve production 
efficiency, reduce costs, improve product quality and environmental safety of processes. The literature describes the 
application of multicriteria optimization for production purposes, including the selection of reaction conditions and 
improvement of technological processes. In the presented paper, the object of research is the process of hydrogenation of 
polycyclic aromatic hydrocarbons (PAH) in the production of high-density fuels. To determine the optimal conditions of 
the process, the problem of multicriteria optimization based on the kinetic model is solved. The criteria include 
maximizing the yield of targeted naphthenes and conversion of feedstock. The research objective is to create a program 
implementing the multicriteria optimization non-dominated sorting genetic algorithm-II (NSGA-IJ). Due to this, it is 
possible to calculate the optimal temperature for the PAH hydrogenation process on the basis of the kinetic model. 
Materials and Methods. The NSGA-II genetic algorithm was used to solve the multicriteria optimization problem. 
Modified parental and survival selection within the Pareto front was also used. If it was necessary to divide the front, 
solutions based on the Manhattan distance between them were selected. The program was implemented in Python. 
Results. In the system of ordinary nonlinear differential equations of chemical kinetics, the concentration was designated 
yi, the conditional contact time of the reaction mixture with the catalyst — t. The system was solved for the hydrogenation 
reaction of polycyclic aromatic hydrocarbons. The calculations showed that at t=0 y,(0)=0.025; y2.(0) =0.9; 
yo(0) = 0.067; yo(0) = 0.008; yi(0) = 0, i = 3-5, 7, 8, 10-20; O(0) = 1. The process temperature was considered as a control 
parameter according to two optimality criteria: maximizing the yield of target naphthenes (f;) at the end of the reaction, 
and maximizing the conversion of feedstock (f:). Values f; were in the range of 0.43-0.79; conversion — 0.01-0.03; 
temperature — 200-300 K. The growth of temperature was accompanied by an increase in the yield of target naphthenes 
and a decrease in the conversion of feedstock. Each solution obtained was not an unimprovable one. When modeling the 
process of hydrogenation of PAH, an algorithm was launched with a population size of 100 and a number of generations 
of 100. A program implementing the NSGA-II algorithm was developed. The optimal set of values of the PAH 
hydrogenation reaction temperature was calculated, which made it possible to obtain unimprovable values of the 
optimality criteria — maximizing the yield of target naphthenes and conversion of feedstock. 

Discussion and Conclusion. The NSGA-U algorithm is effective for solving the problem of non-dominance, and deriving 
the optimal solution for all criteria. Future research should be devoted to the selection of optimal algorithm parameters to 
increase the speed of the solution. Based on the obtained theoretical optimal conditions of the PAH hydrogenation 
reaction, it is possible to implement the process in industry. 


Keywords: hydrogenation of polycyclic aromatic hydrocarbons, multicriteria process optimization, nonlinear 
programming problem, Pareto front, NSGA-II method 


© Alexandrova AA, Koledin SN, 2024 


Information technology, computer science and management 


109 


http://vestnik-donstu.ru 


110 


Alexandrova AA, et al. Optimal Temperature Calculation for Multicriteria Optimization of the Hydrogenation ... 


Acknowledgements. The authors would like to thank the reviewers for valuable comments that contributed to the 
improvement of the article. 


For citation. Alexandrova AA, Koledin SN. Optimal Temperature Calculation for Multicriteria Optimization of the Hydrogenation 
of Polycyclic Aromatic Hydrocarbons by NSGA-II Method. Advanced Engineering Research (Rostov-on-Don). 
2024;24(1):109-118. https://doi.org/10.23947/2687-1653-2024-24-1-109-118 


Hayunaa cmamba 
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THAPHpoOBaHnA NOJMUMKIM4eCKHX APOMATH4eCKHX YrieBOLOpoOA0B MeTo.0OM NSGA-IT 
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AHHOTalna 

Beedenue. MuorokputepuaibHy!0 ONTHMV3allHIO C y4YeTOM MIpOTHBOpeyalHx pyr Apyry Kputepves 3ayjevicrBy!OT WIA 
yryuWieHua 9deKTHBHOCTH TIPOv3BOACTBa, COKpall[eHHA 3aTpaT, MOBbIMeHHA Ka4yecTBa MIPOAYKUMA WU IKONOTMYeCKOH 
OesoracHocTH mpoyeccos. B smTepatype onmvcaHoO McHOsb30BaHHe MHOFOKPHTepHabHOM ONTHMM3alWM [IA 
TIPOH3BOACTBeHHBIX [esieH, B TOM 4HCIe Mp BbIOOpe yCNOBMM peakKUHH MW yIyYWIeHHH TexHOsOrM4ecKHXx TIporeccos. B 
TIpey{cTaBseHHow paOoTe oOObeKT UcceqoBaHHuA — 9TO MIpollecc THAPHpOBAaHHA TOJIMUMKIMYeCKHX apoMaTH4eCcKHXx 
yrmeBoyopoyos (IIAY) mp nomyseHHu BbICOKOMIOTHBIX TorIMB. Jd ONpexeueHHA ONTHMAaJIbHBIX ycuOBHi mpolecca 
pellaerca 3aa4a MHOFOKpHTepHasIbHOH ONTHMV3alMH Ha OCHOBe KHHeTHYeCKON MoyerM. Kputepuu: MaKCHMV3alMA BEIXOa 
eIeBbIX Ha:bTeHOB HM KOHBepcHa UcxogHoro cbipps. Lens paboTsr — co3qaHHe MporpaMMBbl, peasiM3yrolle aropuTM 
MHOTOKpHTepHasIbHoH onTHMv3anMH NSGA-II (aura. non-dominated sorting genetic algorithm II). baarogapa 9Tomy Ha 
OCHOBe KHHETHYECKOM MOJeIM MOXKHO pacCuHTaTb ONTHMAJIbHY!O TeMIepatypy ia Ipowecca rugpupoBanua TAY. 
Mamepuanoi u memoooi. na pemienua MHOTOKpUuTepHabHOM 3aqa4dH ONTHMH3alHH MpHMeHsIcA reHeTHueCKHH 
anroputm NSGA-II. Ucnomp3yetca Takxke H3MeCHCHHBbIM OTOOp posHTesel M BbDKMBaHHA B paMKax (ppouta Iapero. IIpu 
HeOOXOAMMOCTH pa3yeseHHA (PpOHTa peleHHA BbIOMpAaIIMCh IO MAHXITTEHCKOMY PaccTOAHHIO MexKy HUMH. IIporpamMa 
peasmm30BaHa Ha a3bIKe Python. 

Peszyibmamei ucciedoeanua. B cucteMe OObIKHOBCHHBIX HeJIMHeMHbIX TupcdepeHuMabHbIxX ypaBHeHHii XHMHUYeCKOH 
KHHETHKH KOHICHTpal{Hi0 OOO3HAYHJIN Yj, YCIIOBHOe BPCMA KOHTAKTAa PeaKIMOHHOH CMeCH C KaTasIHM3aTOpoM — Tt. CucTemy 
peluwin JUL peaklHVM TH {pHpoBaHHA MOMUMKIMYeCKHX apOMATHYeCKHX YIIeBOLOpONOB. PacueTsI NoKa3zasiM, 4To pu t = 0 
yi(0) = 0,025; y2(0) = 0,9; ye(0) = 0,067; yo(0) = 0,008; yO) = 0, i= 3-5,7,8,10—-20; O(0) = 1. B kayectTBe ynpaBaaemoro 
TlapaMetpa paccMaTpHBaIM TeMilepatypy Mpollecca 10 AByM KpHTepHaM ONTHMasIbHOCTH: MaKCUMU3alHA BbIXO[a LeJIeBbIX 
Ha(pTeHos (/|) B KOHIe peak HW MaKCHMH3al{A KOHBeEpCHM MCXOHOPO CbIpHA (f2). 3Hauenua /| ObuIM B rpannax 0,43—0,79; 
KoHBepcun — 0,01—0,03; Temneparyps1 — 200-300 K. Pocr Temnepatypbi CompoBoxaeTcaA YBeIM4eHHeM BbIXOJ{a LeJIeBbIX 
Ha(TeHOB HM CHWKeHHeM KOHBepcHH UcxoqHoro cbIppa. Kaxgoe nomryyeHHoe pelieHHe — Heymyumaemoe. IIpu 


MOJeMpoBaHHM rmporecca ruapwposanua TIAY 3anycrusm anropHTm c pa3mMepom nomymanaH1 — 100, kommuecrBoM 
nokoneHui — 100. PaspaOotana mporpamma, peasmM3yrolad anmroputm NSGA-II. PaccuntaHo onTHMasIbHoe MHO?KeCTBO 
3HaveHH TeMUepaTypbI peak rHapupoBanua TIAY, no3zpossroujee MOyYTh HeyyYWaeMble 3Ha4eHHA KPHTepHeB 
ONTHMaJIBHOCTH — MaKCHMV3al[HH BbIXOJa IeJIeBbIX Ha()TeHOB H KOHBeEpCHH HCXOHOYO CbIpbsA. 

O6cystcdenue u 3aknio“uenue. Anroputm NSGA-II spdextuBen Jia pewieHua 3aa4n HeJOMHHMpOBaHHA UM BbIBOTAa 
ONTHMAaJIbHOrO pelleHHA JIA BCeX KpHTepveB. byzyuMe uccieqOBaHHA ClleyeT MOCBATHTb MOAOOpy ONTHMAJIbHBIX 
TlapaMeTpoOB aITOpHTMa, MO3BOJIMIOWIMX YBeJIMYMTb CKOpOCTb pelieHua. Onupasch Ha Mosy4eHHble TeopeTH4ecKne 
ONTHMAaJIbHbIe YCIOBHA peakHH ruzpupoBanua ITA Y, MoxHO peasM30BaTb IpOWecc B MpOMBILIJIeHHOCTH. 


Kunouesbie c10Ba: THApHpoBaHve MOJMUMKIMYeCKUX aAPOMATH4eECKUX yTIeBONOpOAOB, MHOrOKpuTepvasibHar 
ONTHMH3alHA TEXHOIOPHYeCKOrO Mpouecca, 3aja4a HesIMHeMHOrO MporpaMMupoBaHna, PpoutT Iapeto, metoy NSGA-II 


Baarogxapnocru. ABTODBI BbIpa7KatOoT IIPH3HAaATCJIBHOCTh PeCICH3CHTaM 3a IWeCHHbIe 3aMedaHHs, cnocoOcTBOBaBILIHe 
YJTY4WIeHHtO CTaTbH. 


Aaa wnTwMpopanna. AnexcanypopaA.A., KonequaC.H. Pacuer ontumMaibHol TemlepaTypbl mp 
MHOFOKpHTepHaIbHOH ONTHMU3allMM Mpolecca THApHpoBaHuA MOMMIIMKIM4eCKUX apOMaTHYeCKUX yIIeBOAOpO0B 
mMetToyom NSGA-II. Advanced Engineering Research (Rostov-on-Don). 2024;24(1):109-118. https://doi.org/10.23947/2687- 
1653-2024-24-1-109-118 
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Introduction. Optimization of multistage reactions is used in chemical, oil and gas, food, and other industries. In 
practice, optimization tasks are multicriteria, and the criteria are often contradictory and have an optimum at different 
points. Multicriteria optimization remains urgent, as it provides taking into account several parameters and selecting the 
best solution from a variety of options. 

Within the framework of this study, the object of research is the catalytic hydrogenation reaction of polycyclic 
aromatic hydrocarbons (PAH). They are a class of organic compounds whose molecules contain at least two benzene 
rings [1]. PAH are common in the interstellar medium. They are part of heavy oil fractions, and are formed by laser 
irradiation of carbon materials. The study of these compounds is of interest from the point of view of establishing the 
dependences between their chemical structure and physico-chemical properties. In addition, the data obtained as a result 
of such scientific research can be used to create new organic and hybrid compounds with a strong carbon frame, which 
are applicable in nanoarchitectonics. 

Applied science correlates the presence of PAH with the production purposes. As an example, the presence of PAH 
in raw materials is desirable if it is used to produce coke with a given structure [2]. However, in the production of fuel, 
PAH can negatively affect the performance characteristics of the product, specifically, the density [3]. 

Extremely strict requirements are imposed on the production of high-density jet fuel. At high density, it should have 
a boiling point no higher than the upper limit of the boiling point of the kerosene fraction. Another mandatory criterion 
is the low content of aromatic hydrocarbons. We also note the high cost of well-known technologies for producing high- 
density fuels. 

Considering the above, the problem must be solved according to the Pareto dominance principle to determine the set 
of unimprovable options using a genetic algorithm of non-dominated sorting [4]. 

The research objective is to develop a program that implements NSGA-II (Non-dominated sorting genetic 
algorithm II) multicriteria optimization algorithm and provides calculating the optimal temperature for the PAH 
hydrogenation process based on a kinetic model. 

Materials and Methods 

Mathematical Model. We describe the changes in the concentrations of the components depending on the reaction 
time. To do this, the equations of chemical kinetics are used, which are a system of ordinary nonlinear differential 
equations (SONDE): 


dyi_X 
OPV yyw i= hl, (1) 
dt & 
j=l 
i I n - 
me ‘yieil_ A... A tees dus i a we a 2 
wj =k; [To k_j [Too »kj = Aj-exp RT »k_j = A_j exp RT | (2) 
Here are the initial conditions: y(0) =y at t=0; te[0, t']; y; — reagent concentrations, mole fractions; 


t — conditional contact time of the reaction mixture with the catalyst, kg-min/mol; J — number of stages; J — number 
of substances; vj; — stoichiometric matrix; w; — speed of the j-th stage, 1/min or mol/(kg-min); kj, k_;— speed constants 
of the stages (given), 1/min; a — negative elements of matrix vj; B— positive elements vj; 4;, A_; — pre-exponential 
multipliers, 1/min; £;*, Ej; — activation energies of direct and reverse reactions, kcal/mol; R — gas constant, 
2 cal/(mol-K); 7 — temperature, K; t" — reaction duration, kg-min/mol. 

The model of catalytic hydrogenation of PAH takes into account the dynamics of the molar composition and volume 
of the reaction mixture. Therefore, changes in the concentration of components at each point in time are considered [3]: 


dO _ wordy: es 
7 aes Q(0) = 0°, 


i=l 
I I 
Wj =k; | | eae le ae | | (=). 
i=] Q i=l Q 


To describe a nonstationary reaction that occurs with a change in the volume of the reaction mixture, it is necessary 
to solve a system of nonlinear differential equations at each moment of time. The direct kinetic problem is the solution to 
SONDE (1)-(3). 

In the process of hydrogenation of PAH, naphthenes are obtained from the source aromatic hydrocarbons, which have 
a higher density and can be used as rocket fuel. For this purpose, nickel catalysts are used, and the control or variable 
parameter is the process temperature, which should be within 200-500 K. Optimality criteria — maximizing the output 
of target naphthenes at the end of the reaction and maximizing the conversion of feedstock. 


(3) 
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Research Algorithm. Multicriteria optimization is the selection of the best solution from a variety of alternatives, 
taking into account several criteria. The importance of each of them is determined by the weight (priority). 
Suppose, f(x) — an objective function, and the constraints given in the form of equalities 1\(x).../4m(x) and inequalities 


2m+1(x)...gp(x) are represented by a column vector of components x = [x1,..., Xn]” in n-dimensional Euclidean space. 
We formulate the problem of nonlinear programming [5]. 
Optimize 
f (x) > extr, x € E" (4) 
with m linear or nonlinear constraints in the form of equalities 
h(x) =.) = laut (5) 
and with (p — m) linear or nonlinear constraints in the form of inequalities 
gj(x)20, j=m+l,...,p. (6) 


Decomposition (4)-(6) is the formulation and solution of a linear and quadratic programming problem. Each of them is 
determined by the type of equations (4)-(6). Thus, in the case of quadratic function (4) and linear equations (5), (6) — the 
quadratic programming problem described below. 

We define the extremum of function 


f (x)=a0+e"x+x' Ox > extr (7) 


with constraints: 
a’x>b,x>0. (8) 


In equations (7), (8) Q — nonnegatively defined quadratic symmetric matrix; a, b, c — coefficient matrices. 
When setting multicriteria Pareto optimization problem, (4)-(6) will have the form: 


extrF (x)= F(x*)=F". (9) 
In equation (9), Fx =(fi(X), fo(X),.....fm(X)) —vector-function of optimality criteria fi and fs. Set x* — desired 


solution to the problem in the region of variation parameters. Set F* represents the desired solution to the problem in the 
region of optimality criteria, which is unimprovable in terms of Pareto approximation. Then, x” defines the Pareto set, 
F* — Pareto front. 

Priori and posteriori Pareto approximation algorithms were used to solve (9). One of them is the ideal point method, 
which represents the best solution according to all criteria [6]. To find it, we must first determine the minimum and 
maximum values of each criterion of all the solutions under consideration. Then, for each criterion, the maximum value 
of all the minimum and the minimum of all the maximum values are selected. 

However, this approach has some drawbacks. Firstly, it may be ineffective if the ideal point is outside the range of 
acceptable criteria values. In this case, other methods of solving multicriteria optimization problems are used. In addition, 
the ideal point method does not take into account the relationship between the criteria and may cause the selection of a 
compromise solution that is not optimal for all criteria. Therefore, when using this method, it is needed to additionally 
analyze and verify the optimality of the solutions obtained [7]. 

The lexicographic ordering method is also used in solving multicriteria optimization problems. In this case, the criteria 
are ordered by priority and considered sequentially. If the solutions cannot be sorted by the first criterion, then they are 
sorted by the next criterion, etc. [7]. The advantages of lexicographic ordering are simplicity and transparency. Using this 
method, it is possible to obtain a single optimal solution that is easily interpreted. As for the disadvantages, we note, 
firstly, the inability to take into account compromise solutions that may turn out to be optimal by all criteria. Secondly, 
there is a risk of selecting an unfavorable solution if the first criterion has a lot of weight, but is not the most important 
for this task [8]. 

The presented work provides a solution to the problem of multicriteria optimization of the hydrogenation process of 
polycyclic aromatic hydrocarbons through the well-known NSGA-II method. It is based on a genetic algorithm and uses 
several techniques to solve the problem of non-dominance [9]. The main steps of the algorithm are described below [10]. 

1. Initialization of the population. The initial population is randomly generated. 

2. Population assessment. Each element of the population is evaluated according to several criteria. 

3. Sorting the population. The elements of the population are sorted by non-dominance level. Dominant and non- 
dominant elements are placed in the first level. Elements dominated only by elements of the first level are placed in the 
second level, etc. 

4. Selecting parent elements. To create a new population, parent elements are selected from the first few levels. 

5. Crossover and mutation. Parent elements undergo crossing-over and mutation to create new elements of the population. 

6. Assessment of the new population. New elements are evaluated according to the criteria. 
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7. Sorting the new population. New items are sorted by non-dominance level. 

8. Selection of a new population. Elements for the next generation are selected from the new population. 

9. Repeat steps 4-8 until the stop criterion is reached. 

NSGA-II allows working accurately and efficiently with multicriteria optimization tasks. It effectively solves the 
problem of non-dominance, which provides getting optimal solutions for all criteria. 

NSGA-II is based on a genetic algorithm with parental selection and survival. Individuals are selected along the fronts, 
while the front is divided if not all individuals can survive. Solutions in the split front are selected on the basis of the 
distance between them, which is the Manhattan distance in the criteria space [9]. The endpoints are saved at each 
generation and are assigned a conditionally infinite distance for use in subsequent iterations [11] (Fig. 1). 


Fig. 1. Visualization of the Pareto front and selection of solutions based on distances 


Figure | shows an example of a set of solutions for a multicriteria optimization problem using criteria f, and fr. The 
red dots represent the Pareto front. The calculation of the cluster distance for solution i is shown — this is the average 
length of the side of the cuboid in which solution i is located (marked with a blue frame). 

To enhance the impact on the selection of parents, NSGA-II uses binary tournament selection [9]. Each individual is 
first compared by rank, and then — by the distance between them. 
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Research Results. The program implementing the algorithm for solving the multicriteria optimization problem is 
written in Python. 

In the course of the study, a system of differential equations (1-3) was solved. For the hydrogenation reaction of 
polycyclic aromatic hydrocarbons, it has the form [12]: 


dy, 
— =(-w, -w.+w 
de ( 1 2 1) 
dj 
- = —3W — w2 —2w3 —3w4 —3ws — 36 —W7 vy = RV _ fen -y 
3 
Ws —3Wo — Wio — Wi — Wi2 — Wi4 — WI Q Q 
ap Wp PI 
SS WI QO? 
dt 
2 
A 45 ancy ie ws y2 kig y7 
dt Q Q 
Laem _ kar yr y3 
dt “ 4 
dy 7 
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geo vos = EELYO-D2 hig + Yio 
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d. Q Q 
ail = W3 — Wo — Wo — 2W16 3 
dt a _ ko Yor ki yio 
6 
a 4 
a W4 + W1i6 g Q 
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7 We = Rai vu0 + y2 
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dt k 3 
dy we a ENS 
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aoe yee kis yi2*y2 
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at 
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dt 
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At t=0, v1(0)=0.025; y2(0)=0.9; y6(0)=0.067; yo(0)=0.008; y(0)=0, =3—5,7,8, 10-20; O(0)=1. 
Figure 2 shows the calculated set of solutions satisfying the constraints of the system and the Pareto front, 
where f| — output of the target naphthenes, and f: — conversion of the feedstock. 
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Fig. 2. Pareto front of the hydrogenation process of polycyclic aromatic hydrocarbons 


The controlled parameter is the process temperature [13]. The optimality criteria are maximizing the output of target 
naphthenes at the end of the reaction and maximizing the conversion of feedstock [14]. Calculations using the proposed 
algorithm gave the results of matching the optimality criteria and the values of the temperature parameter, which are 
summarized in Table 1. 
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Table 1 
Compliance of optimality criteria and temperature parameter value 
of the hydrogenation of polycyclic aromatic hydrocarbons 


Output of targeted naphthenes /; Feedstock conversion f2 Temperature T, K 
0.43 0.03 200.00 
0.59 0.03 221.20 
0.71 0.01 250.00 
0.76 0.01 271.10 
0.79 0.01 300.00 


Thus, as the temperature grows, the output of the target naphthenes increases, and the conversion of the feedstock 
decreases. Each of these solutions is unimprovable, and the selection of specific values depends on the decision maker. 

For NSGA-II to work successfully, it is required to select the algorithm parameters correctly. Specifically, the criteria 
of population size, number of generations, probability of crossover and mutation, should be optimized. When launching 
the algorithm on the model of the hydrogenation process of polycyclic aromatic hydrocarbons, the following parameters 
were used: population size — 100, number of generations — 100. An insufficiently large population may cause premature 
convergence of the algorithm to the local optimum. Too large population can slow down the optimization process [15]. 
An excessive number of generations potentially leads to retraining the algorithm, whereas with an insufficient number of 
generations there may not be enough time to achieve optimal solutions. 

Discussion and Solution. A program has been created that implements the NSGA-II multicriteria optimization 
algorithm. Working with the corresponding problem within the framework of this method includes solving a system of 
differential equations, visualizing a set of solutions satisfying the constraints of the system, and constructing a Pareto 
front. In addition, the values of the variable parameters were found to achieve optimization goals. For the PAH 
hydrogenation process, based on the kinetic model, a set of temperature values has been calculated that are optimal for 
obtaining the unimprovable values of two optimality criteria: naphthene output and feedstock conversion. With increasing 
temperature, the reaction rate and the output of naphthenes increases. However, the conversion of raw materials is 
decreasing. In addition, too high temperatures can cause adverse reactions and decomposition of products. 

The data obtained in the framework of the presented research can be useful for optimizing the process of hydrogenation 
of PAH under industrial conditions. It is important to take into account the impact of temperature on the output of 
naphthenes and the conversion of raw materials when developing a production strategy. In addition, other parameters, on 
which the kinetics of the reaction depends, should be considered. These are, e.g., pressure, flow rate of reagents, and the 
role of catalysts. 

Thus, the developed program and the proposed algorithm make it possible to simultaneously analyze several 
criteria for the optimality of the process based on a kinetic model, and generate a set of unimprovable values of 
variable parameters. 
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O6 aemopax: 
Anactacna AjeKcaHsApoBHa AslekcaH{poBa, MarvcTpaHT Kadepbl HHPOpMalMOHHBIX TEXHOJOPHii H NpHKaqHOH 


MaTemaTHKH Y®uMckoro rocyyfapcTBeHHoro HeiTaHoro TexHHuecKoro yHuBepcuteta (450064, P®, r. Yoda, 
yn. KocmMonastos, |), SPIN-Kom: 4026-5240, ORCID, nastenal425@gmail.ru 


Cepreii Hukonaesuy Kosequn, KanqugatT :u3ukKO-MaTeMaTHYeCKHX HayK, JOWeHT Kadespbl UHPoOpMalMOHHBIX 
TeXHOOrHH WM WpukTaqHOW MaTeMaTukKH YduMcKoro rocyapcTBeHHOTO He@TAHOTO TexXHHYeCKOrO YHMBepcHuTeta 
(450064, P®, r. Ycba, yn. Kocmonastop, 1), SPIN-Kom: 4243-6265, ORCID, koledinsrg@gmail.com 


3aAeleHHblu 6K1aO aemopoe: 
A.A. AnekcaHapopa — pa3paboTka nporpaMMHOro oGecrtedeHHA, MOATOTOBKa TeKcTa, (POPMYJIMPOBKa BbIBOOB. 
C.H. Konequn — Hay4Hoe pyKOBOJCTBO, MpeOcTaBsIeHHe HCXOHBIX JJAHHBIX, KOPpeKTHPOBKa BBIBOOB, JOpaboTka TeKcTa. 


Konqbauxm unmepecos: aBropbl 3a BIAIOT OO OTCYTCTBHM KOHAMKTAa MHTepecos. 
Bce aemopoi npowumaau u odo6puau OKOHYaMeNbHbItL BaPUaHM PyKONUCU. 
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